BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,930	T→A	46.2%	intergenic (‑1013/‑155)	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,930	0	T	A	46.2%	37.5 / 213.7	156	intergenic (‑1013/‑155)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (46/38); new base A (39/33); total (85/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 157 total aligned reads displayed.

ATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATT  >  CP000730/1929793‑1930068

attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTGATCAGTGGCTGTGCTGTTGCGTATTATACttt                                                                                                                                               >  2:513247/8‑140 (MQ=255)

   tGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGACGCCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTACGTATTTCTCTTAtt                                                                                                                                             >  2:175705/1‑141 (MQ=255)

    caagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                           <  1:863092/138‑1 (MQ=255)

    acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATTGAAGAGGCATGGACG‑GGCAT‑‑TAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                           <  1:823026/138‑1 (MQ=255)

    acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑‑TCGACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTGTACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                           <  2:994198/138‑1 (MQ=255)

      cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                                          >  1:1050144/1‑141 (MQ=255)

      agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGCGGAGCTTCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                                         <  1:1040905/140‑1 (MQ=255)

       gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt                                                                                                                                         >  1:542007/1‑141 (MQ=255)

        atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa                                                                                                                                        <  2:931290/141‑1 (MQ=255)

         tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                                      <  1:175163/138‑4 (MQ=255)

         tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑G‑T‑T‑GC‑TAATCATTTTCGTTCACATTATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                      >  2:583068/4‑141 (MQ=255)

         gggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                                      <  1:804605/138‑4 (MQ=255)

          ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                     <  1:655840/139‑1 (MQ=255)

          ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                     <  2:186077/139‑1 (MQ=255)

          ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                     <  2:156225/139‑1 (MQ=255)

          ggtataAGAATCGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                     <  2:33520/139‑1 (MQ=255)

           gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGACTTAGTa                                                                                                                                    >  1:59171/2‑141 (MQ=255)

           gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                                    >  2:18345/2‑141 (MQ=255)

           gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                                    <  2:539936/140‑1 (MQ=255)

            tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                   >  1:337062/1‑140 (MQ=255)

             ataataataGGGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                  <  2:838248/141‑3 (MQ=255)

             ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                  <  2:911961/141‑3 (MQ=255)

              aaataatAGAGCCGCCTATCGCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                                 <  2:873598/140‑4 (MQ=255)

                ataagagagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                               <  2:178958/136‑1 (MQ=255)

                   atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑aa                                                                                                                          >  2:475339/2‑140 (MQ=255)

                    tagaggcctcTATCTCTCAGGTGTCAATT‑‑‑GACGAAGTGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTGCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTCTTGGCTTAGTAAACGCa‑‑aca                                                                                                                         <  1:746053/132‑3 (MQ=255)

                    tagAGTCGCCTATCTCTCAGGCGTCAA‑‑‑TGGACGAAGAGAGGAGGTGCA‑TTA‑G‑T‑TGGC‑TAATCATTTTCGTGCACGTCATAGCACCAGTCACCAGTGGCCGTGCTGTGGCGTATTATACTTGTTGGCATCGTAGACCCa‑‑aca                                                                                                                         >  2:548803/1‑139 (MQ=255)

                     agagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                        >  2:810943/1‑138 (MQ=255)

                     aaagTCGCCTATCTCTCAGGCAGCAA‑‑‑TGGACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                        <  2:32700/139‑4 (MQ=255)

                      gagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                       <  1:203723/141‑5 (MQ=255)

                      gagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                       >  1:218811/1‑137 (MQ=255)

                       agTCGCCTATCTCTCAGGCGGCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCGTAGCACCAGGCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCA‑‑Acacat                                                                                                                      >  1:610664/1‑140 (MQ=255)

                        gTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGGGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCTCCAGTCATCAGTGGCTGTGCTGTTGCGTATGATACTTATTGGCTTAGTAAAAGCa‑‑acaaata                                                                                                                     >  2:970926/1‑135 (MQ=255)

                         tCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataa                                                                                                                    <  1:901024/141‑8 (MQ=255)

                          cGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTCAGTAGACGCA‑‑ATACaa‑‑aa                                                                                                                  >  2:974809/1‑141 (MQ=255)

                          cGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACaa‑‑aa                                                                                                                  >  1:138422/1‑141 (MQ=255)

                           gCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AAt                                                                                                                 >  1:598935/1‑141 (MQ=255)

                           gCCTATCTCTCAGGCGTCAAAT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAGTCATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATCATACTTATTGGCTTAGTACACGCa‑‑acaaataatt                                                                                                                  >  1:645268/1‑132 (MQ=255)

                            ccTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACaa‑‑aaga                                                                                                                >  1:830701/1‑139 (MQ=255)

                            ccTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATa                                                                                                                >  1:767668/1‑141 (MQ=255)

                            ccTATCTCTCAGGCGTCAATT‑‑‑GACGAAGATAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                 <  2:42778/141‑11 (MQ=255)

                             cTATCTCTCAGGCGTCAATT‑‑‑GGCGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                >  1:551886/1‑130 (MQ=255)

                             cTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGCGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATCGCACCAGTCCTCAGTGGCTGTGCTGTTTCGTATTAGACTTATTGGCTTAGTAAACGCacgaagtaattat                                                                                                                  >  2:298574/1‑129 (MQ=255)

                             cTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                >  1:975657/1‑130 (MQ=255)

                             cTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                >  1:699709/1‑130 (MQ=255)

                              tATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑TACTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGGATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaaa                                                                                                               >  1:438285/1‑129 (MQ=255)

                               aTCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGt                                                                                                             >  1:1023345/1‑141 (MQ=255)

                               aTCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGACATTGCGTTTTTTTCTTATTGGTTAAGTAGACGCA‑‑CTATAA‑‑AATAggg                                                                                                             >  2:958987/1‑140 (MQ=255)

                               aTCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCGCATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                              >  2:605352/1‑128 (MQ=255)

                               aTCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                              >  1:189835/1‑128 (MQ=255)

                               aTCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                              <  2:256161/141‑14 (MQ=255)

                                tctctcAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatag                                                                                                             <  2:1062047/141‑15 (MQ=255)

                                tctctcAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatag                                                                                                             >  2:852191/1‑127 (MQ=255)

                                tctctcAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatag                                                                                                             <  2:949714/141‑15 (MQ=255)

                                 ctttcAGGCGTCAATT‑TAGACGCAGAGAGTAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGa                                                                                                           <  1:286633/138‑1 (MQ=255)

                                  tctcAGGCGTCAATT‑TACACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGCGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAc                                                                                                          <  2:709455/141‑1 (MQ=255)

                                     cAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACa‑ta                                                                                                      <  1:399780/141‑1 (MQ=255)

                                      aGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATCCAA‑‑AAAAGGTGCCa‑tat                                                                                                     >  1:76182/1‑141 (MQ=255)

                                           tCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑g                                                                                              >  1:904700/1‑141 (MQ=255)

                                             aaTT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑CTACAA‑‑AATAGGTGCCA‑TAT‑AGCC‑Gcc                                                                                            >  2:708616/1‑140 (MQ=255)

                                             aaTT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCa                                                                                            >  2:925655/1‑141 (MQ=255)

                                               tt‑tAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑CTAGAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCAcc                                                                                          >  1:487394/1‑141 (MQ=255)

                                               tt‑tAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCCTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCAcc                                                                                          <  2:836007/141‑1 (MQ=255)

                                                    gatgCAGAGTGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATa                                                                                      <  2:672060/138‑1 (MQ=255)

                                                     aCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATGGATTTTGTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATaa                                                                                     >  1:771971/1‑141 (MQ=255)

                                                         agagagGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGCCA‑TAT‑AGCA‑GCCCCAATAAAAAt                                                                                 >  1:132380/1‑141 (MQ=255)

                                                            gaggagGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATccc                                                                              <  1:226373/141‑1 (MQ=255)

                                                                 ggTGTA‑TAA‑G‑G‑T‑GA‑TGCTTACTGTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                                         <  2:308173/141‑1 (MQ=255)

                                                                    gtA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTAc                                                                      >  2:735290/1‑141 (MQ=255)

                                                                         aa‑G‑G‑T‑GG‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                                  >  2:350452/1‑141 (MQ=255)

                                                                         aa‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                                  <  2:1030728/141‑1 (MQ=255)

                                                                          a‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTGTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaa                                                                 >  2:908138/1‑141 (MQ=255)

                                                                            g‑g‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaaa                                                                <  2:937701/141‑1 (MQ=255)

                                                                                 ggA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATa                                                              <  2:346207/140‑1 (MQ=255)

                                                                                   gatgttgtttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAg                                                             <  1:496544/133‑1 (MQ=255)

                                                                                   a‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGt                                                            <  2:407981/141‑1 (MQ=255)

                                                                                        ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTATTTTCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATACAAATCCCCTCACTACCGCAAATAGTGAgg                                                        >  2:11435/1‑141 (MQ=255)

                                                                                        ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                        <  1:511082/141‑1 (MQ=255)

                                                                                        ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                        <  2:188102/141‑1 (MQ=255)

                                                                                        ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTGTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                        >  1:1007127/1‑141 (MQ=255)

                                                                                         tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAGTGCGTTTTTTTGTTATTGGCTAAGTAGCCGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCCCTACCGCAAATAGTGAggg                                                       >  1:407725/1‑141 (MQ=255)

                                                                                          aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                      >  2:764792/1‑141 (MQ=255)

                                                                                           ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                     <  2:44583/141‑1 (MQ=255)

                                                                                             ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                                   <  1:238397/141‑1 (MQ=255)

                                                                                              tCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGGGGGGAtgg                                                  >  1:396407/1‑139 (MQ=255)

                                                                                                    acaTCATAGCACCAGTCATTAGTGCCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtat                                            <  1:858724/141‑1 (MQ=255)

                                                                                                         atAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcat                                            >  1:415073/1‑141 (MQ=255)

                                                                                                         atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGTGGTATAAGTa                                       >  2:243478/1‑141 (MQ=255)

                                                                                                           aGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACC‑TAT‑AGCC‑GCACCCATAAAATGCCCCGCCATTCCGCAAATAGTGAGGGGATTGGTGTAtaaataaa                                     >  2:744662/1‑141 (MQ=255)

                                                                                                            gCACCAGTCATCAGCGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGAACTTAGTGGGGGGATGGGTGTATAAGTAAAt                                    >  1:20602/1‑141 (MQ=255)

                                                                                                                    catcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                            >  2:726838/1‑141 (MQ=255)

                                                                                                                    catcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                            <  1:659974/141‑1 (MQ=255)

                                                                                                                    aatcgttgccTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                            <  1:392388/132‑1 (MQ=255)

                                                                                                                      tcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTc                          >  1:3721/1‑141 (MQ=255)

                                                                                                                        aGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGCCG‑TGT‑ATCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGt                        <  1:1042748/141‑1 (MQ=255)

                                                                                                                                  ccATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAtt              >  1:957166/1‑141 (MQ=255)

                                                                                                                                     ttGCCTTGTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑GAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATtata           <  2:147131/141‑1 (MQ=255)

                                                                                                                                       gCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACt         >  1:462097/1‑141 (MQ=255)

                                                                                                                                       gCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGAC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACt         >  2:297797/1‑141 (MQ=255)

                                                                                                                                              tttCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAGTTTCGTTGTCTTAATTATACTGCTAAtt  >  1:599268/1‑141 (MQ=255)

ATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATT  >  CP000730/1929793‑1930068

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: `‑`