BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,917	C→T	51.7%	intergenic (‑1000/‑168)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,917	0	C	T	51.7%	‑3.0 / 152.2	114	intergenic (‑1000/‑168)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (24/31); new base T (38/20); total (63/51)
Fisher's exact test for biased strand distribution p-value = 2.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

Only 100 of 114 total aligned reads displayed.

GACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTT  >  CP000730/1929779‑1930052

gACTGTGCTTGTATATTTGCCGATCATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCCCATCATAGCACCAGTCAGCAGTGGCTGTGCCAt                                                                                                                                       <  1:134979/141‑1 (MQ=255)

   tatgattGTATATTCACAGATGGTATAATAATAGAATCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TT‑GGAGCATCATTGTCGGTCACAGCATAGCCCCAGTCATCAGTGGCCGTGCCGGTGCg                                                                                                                                   >  2:41985/6‑141 (MQ=255)

   tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCGGTGGCTGTGCTGTTGCg                                                                                                                                   >  1:117086/3‑141 (MQ=255)

    ttgCTTGTATAGTTGCCGATGATATAATAGTAGCGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                                                                   <  2:101255/140‑1 (MQ=255)

    gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                                                                   >  1:416956/1‑141 (MQ=255)

     tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGta                                                                                                                                 >  1:32426/1‑140 (MQ=255)

       cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                                                               <  1:557506/141‑4 (MQ=255)

        ttGTATATTCACAGATGGGATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                              <  1:594485/141‑2 (MQ=255)

          gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAA‑TCATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCCGTGCTGTTGCGtattata                                                                                                                            >  2:374439/1‑138 (MQ=255)

           tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtgttatac                                                                                                                           >  1:187644/11‑137 (MQ=255)

           tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatat                                                                                                                           >  2:605751/11‑137 (MQ=255)

           tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTGGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                                                           >  1:451585/11‑137 (MQ=255)

           tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGACGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TAATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                                                           >  2:605744/11‑137 (MQ=255)

           tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATTC‑TTATCTTCGTTCGCATCATAGCACCAGTCATCAGTGCCTGTGCCATTGCGttttttt                                                                                                                            <  1:407313/141‑1 (MQ=255)

            atattcacagatggtATAATAATAGAGTCGCGTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTACACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                                                                          >  1:491760/10‑136 (MQ=255)

             tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACtt                                                                                                                         >  1:303874/9‑141 (MQ=255)

             aatTTGCCGATGATATAATAGGAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt                                                                                                                          <  2:262162/140‑1 (MQ=255)

              tatcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                                                        >  2:337992/8‑141 (MQ=255)

              attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                                                        >  2:121783/8‑141 (MQ=255)

              atTTGCCGATGATATAATAGTAGAGTCGACTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                                                                                                         <  1:463886/141‑1 (MQ=255)

                  acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGCGGAGGTGCATTAC‑T‑TT‑GCTAG‑TCCTTGTCGCTCACATCAGAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                    >  2:640275/4‑141 (MQ=255)

                   cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                   <  2:591541/139‑1 (MQ=255)

                   cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                   >  1:9688/3‑141 (MQ=255)

                    cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                   <  2:40546/141‑1 (MQ=255)

                       tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                               >  1:183012/4‑141 (MQ=255)

                        ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATTATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                              <  2:20863/139‑1 (MQ=255)

                        ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTCCGTATTGTACTTCTTGGCTTAgg                                                                                                              >  2:2184/3‑140 (MQ=255)

                        ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                              >  2:621136/3‑141 (MQ=255)

                        ggtacaatgttagAGTCGCCTATCTCCTAGGCGGCACTT‑‑GCCGAAGAGAGGAGGTGCATTCC‑T‑TT‑GCTAG‑TCATTTCCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                              <  2:112381/131‑1 (MQ=255)

                         gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGCTCACCTCAAAACACCTGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                             >  1:126114/2‑141 (MQ=255)

                         atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                              >  2:647888/1‑141 (MQ=255)

                           ctaataaGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                           <  1:533308/140‑3 (MQ=255)

                           ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                           >  1:328466/1‑139 (MQ=255)

                           ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                           >  1:145633/1‑139 (MQ=255)

                            taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                          <  1:220128/141‑4 (MQ=255)

                            taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑GCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATAATTATTGGCttggtaaac                                                                                                          >  1:91844/1‑133 (MQ=255)

                             aataatACAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                         <  2:430551/141‑1 (MQ=255)

                              ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTacctcgacc                                                                                                         >  2:429641/1‑133 (MQ=255)

                               taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TAATTTTCGTTCACATCATAGCAGCAGTCATCAAAGGATGTGCTGTTGCGTCTTATACTTATTGGCTTAGTAACCGCa                                                                                                       >  2:281208/4‑141 (MQ=255)

                               aaatagAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                       <  2:581824/138‑1 (MQ=255)

                                  tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                    <  1:47562/141‑3 (MQ=255)

                                  tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                    <  2:299261/141‑3 (MQ=255)

                                  tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                    >  1:293777/1‑139 (MQ=255)

                                  tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                    >  1:2819/1‑139 (MQ=255)

                                    gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                   <  1:459516/141‑1 (MQ=255)

                                    gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                  >  2:282973/1‑137 (MQ=255)

                                      gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAACCGCaacaaata                                                                                                >  1:213076/1‑135 (MQ=255)

                                      gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                >  2:430302/1‑135 (MQ=255)

                                       tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                               >  2:497279/1‑134 (MQ=255)

                                         gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                             >  2:154913/1‑132 (MQ=255)

                                         gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                             <  1:464365/141‑10 (MQ=255)

                                           cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCGGTCGTCATCGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                           >  1:154206/1‑130 (MQ=255)

                                           cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGGGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                           <  2:45839/141‑12 (MQ=255)

                                            tacccctcAGGCGTAAATTTAGACGCAGAGAGGAGGTGGA‑TAA‑G‑GT‑GAGGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAgg                                                                                           <  1:190978/136‑1 (MQ=255)

                                            tATCTCTCAGGCGTCAATT‑‑GACGAAGAGCGGAGTGGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATTAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                          <  2:372869/141‑13 (MQ=255)

                                            tATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCA‑TAACT‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCCCCAGTCATCAGGGGCTGTGCTGTTGCGTATTATACTTATTGGCTGAGTAAACGCaacaaataattaat                                                                                          <  1:603319/141‑13 (MQ=255)

                                             aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAC‑TCATTTTCTTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                         <  1:349945/141‑14 (MQ=255)

                                             aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaatccatactc                                                                                         >  2:328789/1‑128 (MQ=255)

                                             aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                         >  2:101483/1‑128 (MQ=255)

                                             aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                         >  1:377099/1‑128 (MQ=255)

                                              tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGCCTAAGTAGCAGCCATACAACATAGGTg                                                                                         >  2:528254/1‑141 (MQ=255)

                                              tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTTGTAAACGCaacaaataattaatag                                                                                        <  1:426032/141‑15 (MQ=255)

                                              tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattcatag                                                                                        >  2:195512/1‑127 (MQ=255)

                                                 ctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACa                                                                                      >  2:389607/1‑141 (MQ=255)

                                                 cccacgCGTCAATTTATCCGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCAGCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACa                                                                                      <  1:401522/136‑1 (MQ=255)

                                                   cAGGCGTCAATTTAGACGAAGAGAGGAGGTGTA‑GAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                                                                    <  2:33137/141‑1 (MQ=255)

                                                       cggcAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGc                                                                                <  2:463808/138‑1 (MQ=255)

                                                             tttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                                                                          <  1:425465/141‑1 (MQ=255)

                                                             tttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTATGTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                                                                          >  2:368015/1‑141 (MQ=255)

                                                                   cGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCGAGGTAGATGCAATACAAAATAGGTGCCGTTTAGCCGCACCATTaga                                                                    >  1:17030/1‑139 (MQ=255)

                                                                    gCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                                                                   >  2:659479/1‑141 (MQ=255)

                                                                    gCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTAACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTGGCTAAGTAGACGGAATACATAATAGGTGACATATAGCCGCACCAATaaaa                                                                   >  2:321984/1‑141 (MQ=255)

                                                                          aggaggTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGACATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcccc                                                             >  1:628334/1‑141 (MQ=255)

                                                                           ggaggTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                            <  2:543100/141‑1 (MQ=255)

                                                                              ggTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                                                         <  1:589440/141‑1 (MQ=255)

                                                                               gtgtA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACt                                                        >  1:231736/1‑141 (MQ=255)

                                                                                  ta‑taA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAcc                                                     >  1:620603/1‑141 (MQ=255)

                                                                                  ta‑taA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAcc                                                     >  1:281535/1‑141 (MQ=255)

                                                                                     taA‑G‑GT‑GATTC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc                                                   <  1:24990/141‑1 (MQ=255)

                                                                                     taA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc                                                   <  1:494761/141‑1 (MQ=255)

                                                                                         g‑gT‑GAGTC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa                                                <  2:232352/141‑1 (MQ=255)

                                                                                              gATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAACATAGGAGTCATAGAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAg                                             <  1:293784/141‑1 (MQ=255)

                                                                                              gATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGACATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAg                                             <  1:437436/141‑1 (MQ=255)

                                                                                              gATGC‑TCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAg                                             <  1:2760/141‑1 (MQ=255)

                                                                                                 gC‑TTATTTTCGTTGCAATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGa                                          <  1:2023/141‑1 (MQ=255)

                                                                                                 gC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGa                                          >  1:65744/1‑141 (MQ=255)

                                                                                                   gttATTTTCGTTCACATCATAGCATCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGGGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                         <  2:619904/140‑1 (MQ=255)

                                                                                                     tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                       >  2:220031/1‑141 (MQ=255)

                                                                                                       ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAACAATCCCCTCACTCCCGCAAATCGTGAgggca                                     <  2:661663/141‑3 (MQ=255)

                                                                                                       ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTAATAATCCCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                     <  1:655242/141‑1 (MQ=255)

                                                                                                          ccGTTCACATCATAGCACCAGTCATCAGGGGCTGTGCCATTGCTTTTTTTTCTTATTGGCGAAGTAGACGCAATACAAAATAAGGGACATATAGCGGCACCAATAAAAAACCCCTCACTACCGCAAATAGTGAGGGGATTg                                  <  2:198947/140‑1 (MQ=255)

                                                                                                             ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                               >  1:320366/1‑141 (MQ=255)

                                                                                                                 catcatAGCACCAGTCATCAGGGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTCAGTAGACGCAATACAAAATAGGTGACATATAGCCGCGCCAATAACAATCCCCTCACTACCGCAAAAAGGGGGGGGATTGGTGtata                           >  2:539819/1‑141 (MQ=255)

                                                                                                                    catAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGGTTGGTGTATAAGt                        >  1:441623/1‑141 (MQ=255)

                                                                                                                           ccAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACt                 <  2:65744/141‑1 (MQ=255)

                                                                                                                            cAGTCATCCGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACtt                <  2:394030/141‑1 (MQ=255)

                                                                                                                                 atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATTCCCCTAGCACCGCAAATCGTGAGGGGCTTGGTGTATAAGTAAATAATTAtttt           >  1:315973/1‑141 (MQ=255)

                                                                                                                                 atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtttt           >  2:291509/1‑141 (MQ=255)

                                                                                                                                        gCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTc    <  1:220031/141‑1 (MQ=255)

                                                                                                                                          tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGAAATATAGCCGAACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATGGGTGTATAAGTAAATACTTATGTTCGTTGTCtt  >  1:263186/1‑141 (MQ=255)

GACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTT  >  CP000730/1929779‑1930052

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`