BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,946	G→A	30.4%	intergenic (‑1029/‑139)	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,946	0	G	A	30.4%	164.6 / 84.0	126	intergenic (‑1029/‑139)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (51/36); new base A (23/15); total (75/51)
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.72e-01

Only 100 of 127 total aligned reads displayed.

ATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929809‑1930086

ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                  <  2:944018/141‑1 (MQ=255)

   atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                                                             <  1:955160/140‑2 (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)

        gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACaa                                                                                                                                         >  1:858012/1‑141 (MQ=255)

        gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                                        <  1:592195/141‑7 (MQ=255)

        gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                                        >  1:40418/1‑135 (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)

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

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

          cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGGAAACGCa‑‑acaaataat                                                                                                                                      >  2:344020/1‑133 (MQ=255)

           gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑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)

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

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

            ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                                    <  2:985247/141‑11 (MQ=255)

             cTATCTCTCAGGCGTCAATT‑‑GGCGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                                   >  1:551886/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)

               aTCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGt                                                                                                                                >  2:846315/1‑141 (MQ=255)

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

               aTCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑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‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                 >  1:418314/1‑128 (MQ=255)

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

               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAAGCa‑‑acaaataattaata                                                                                                                                 >  1:712638/1‑128 (MQ=255)

               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGACTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                 >  2:89379/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‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTTTTGTCTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaagag                                                                                                                                <  2:437983/141‑15 (MQ=255)

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

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

                 ttctgAGGCTTCAAGTTAGACGCAGAGAGGAGGGGTA‑GAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCTCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGa                                                                                                                              <  1:441641/140‑1 (MQ=255)

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

                  tctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAc                                                                                                                             >  2:238397/1‑141 (MQ=255)

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

                     cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑a‑ta                                                                                                                        <  2:178970/141‑1 (MQ=255)

                     cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑a‑ta                                                                                                                        <  1:399780/141‑1 (MQ=255)

                      aGGCGTCAATTTAGACGCGGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑CTACAA‑‑AATAGGTGAA‑a‑tat                                                                                                                       >  1:210056/1‑141 (MQ=255)

                      aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATCCAA‑‑AAAAGGTGCC‑a‑tat                                                                                                                       >  1:76182/1‑141 (MQ=255)

                             aaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCa                                                                                                              >  1:608925/1‑141 (MQ=255)

                               tttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑CTAGAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCAcc                                                                                                            >  1:487394/1‑141 (MQ=255)

                                  aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAAt                                                                                                         <  1:794324/141‑1 (MQ=255)

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

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

                                      gCAGGGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCTCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATaaaa                                                                                                     <  1:345224/141‑1 (MQ=255)

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

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

                                             ggaggGGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                                              <  2:292634/141‑1 (MQ=255)

                                                ggTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                                                           <  1:23388/141‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                      cgtATGTTCTTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCCTTTTTTTCTTATTGGCCAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCCCTCTACCGCAAATAGTGAg                                                                           <  2:724028/139‑1 (MQ=255)

                                                                      cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                                           >  2:468686/1‑141 (MQ=255)

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

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

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

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

                                                                        ttttttCGTTCCCATCATAGCACCAGTCATCAGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                         <  1:175705/139‑1 (MQ=255)

                                                                        tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                         >  1:55913/1‑141 (MQ=255)

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

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

                                                                             tCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AGACAA‑‑GATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAGTg                                                                    >  2:528525/1‑141 (MQ=255)

                                                                              cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTgg                                                                   >  2:183628/1‑141 (MQ=255)

                                                                                   acaTCATAGCCCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtat                                                              <  2:301165/141‑1 (MQ=255)

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

                                                                                   acaTCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTT‑TATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGAGGGGATTgg                                                                   >  2:17004/1‑141 (MQ=255)

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

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

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

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

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

                                                                                                    atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCGACTAGTGAAGGGGGTGGTGGAGAGGTAAGTACtgcttat                                             >  1:325314/1‑135 (MQ=255)

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

                                                                                                         tGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTg                                        <  2:1050144/141‑1 (MQ=255)

                                                                                                          ggCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGATGTAGTAAATACTTATTTTCGTTGt                                       <  2:1023864/141‑1 (MQ=255)

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

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

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

                                                                                                                                     tGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtttt            <  2:608925/141‑1 (MQ=255)

                                                                                                                                         aaaGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTa        <  1:774601/140‑1 (MQ=255)

                                                                                                                                            gTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAACATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTCATTTTTCTTTTTGTaaaa     >  2:567359/1‑141 (MQ=255)

                                                                                                                                            gTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTGATTATACTGCTAATTTTTCTTTTTGTaaaa     >  1:947578/1‑141 (MQ=255)

                                                                                                                                            gTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaaaa     >  2:158463/1‑141 (MQ=255)

                                                                                                                                            gTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGAC‑GCGCCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaaaa     >  1:944107/1‑141 (MQ=255)

                                                                                                                                             tAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACGAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGCTTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTTTTTTTGTAAAat    >  2:695969/1‑141 (MQ=255)

                                                                                                                                               gACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAatat  >  2:689046/1‑141 (MQ=255)

ATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929809‑1930086

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