Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 572,129 T→G 12.4% pseudogene (153/1040 nt) nmpC ← DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273572,1290TG12.4% 747.8 / 107.8 293pseudogene (153/1040 nt)nmpCDLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base T (117/139);  new base G (20/16);  total (137/156)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 294 total aligned reads displayed.

CAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTT  >  CP009273/571981‑572278
                                                                                                                                                    |                                                                                                                                                     
 aaGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGc                                                                                                                                                     >  1:4350669/1‑150 (MQ=255)
 aaGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGc                                                                                                                                                     >  1:4449529/1‑150 (MQ=255)
 aaGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGc                                                                                                                                                     >  2:3777217/1‑150 (MQ=255)
 aaGACGGGTTTTATCTTTGGAGGAACCTCGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGc                                                                                                                                                     >  2:3683228/1‑150 (MQ=255)
  aGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGca                                                                                                                                                    >  1:2237450/1‑150 (MQ=255)
  aGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGca                                                                                                                                                    >  1:1737191/1‑150 (MQ=255)
  aGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGca                                                                                                                                                    >  1:2517335/1‑150 (MQ=255)
     cGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGcatca                                                                                                                                                 >  2:3487456/1‑150 (MQ=255)
       ggTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCAtt                                                                                                                                               >  1:4178381/1‑150 (MQ=255)
         ttttATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTag                                                                                                                                             <  1:1000544/150‑1 (MQ=255)
         ttttATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTag                                                                                                                                             <  1:4502115/150‑1 (MQ=255)
           ttATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTagag                                                                                                                                           >  2:945479/1‑150 (MQ=255)
             aTCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAgagga                                                                                                                                         >  1:126258/1‑150 (MQ=255)
              tCTTTGGATGAACCTTGAGATTCAGCTCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGTTCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAgaggag                                                                                                                                        <  1:1970595/150‑1 (MQ=255)
              tCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAgaggag                                                                                                                                        <  2:3797971/150‑1 (MQ=255)
                 ttGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAag                                                                                                                                     <  2:2731658/150‑1 (MQ=255)
                   ggaggaACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGGCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAagta                                                                                                                                   >  1:862918/1‑150 (MQ=255)
                    gaggaACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGgagagagta                                                                                                                                  >  1:3829163/1‑145 (MQ=255)
                      ggaACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTg                                                                                                                                >  2:1880594/1‑150 (MQ=255)
                       gaACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGc                                                                                                                               >  1:1443635/1‑150 (MQ=255)
                         aCCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCtt                                                                                                                             >  2:2931478/1‑150 (MQ=255)
                            ttGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGc                                                                                                                          >  1:1824293/1‑150 (MQ=255)
                             tGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCa                                                                                                                         <  2:1261826/150‑1 (MQ=255)
                              gagaTTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCAt                                                                                                                        >  1:1144531/1‑150 (MQ=255)
                              gagaTTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCAt                                                                                                                        <  1:3469572/150‑1 (MQ=255)
                               agaTTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCAtt                                                                                                                       >  2:4278451/1‑150 (MQ=255)
                                 aTTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTAACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTaa                                                                                                                     <  1:1066238/150‑1 (MQ=255)
                                   tCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTTGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATAACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACt                                                                                                                   >  1:737054/1‑150 (MQ=255)
                                     aGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACttt                                                                                                                 <  2:533312/150‑1 (MQ=255)
                                     aGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACttt                                                                                                                 >  2:2082851/1‑150 (MQ=255)
                                     aGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACttt                                                                                                                 <  2:1606814/150‑1 (MQ=255)
                                        gcgGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGTTTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTccc                                                                                                              <  1:1113805/150‑1 (MQ=255)
                                          ggTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGt                                                                                                            >  1:3051431/1‑150 (MQ=255)
                                          ggTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGt                                                                                                            >  1:1624390/1‑150 (MQ=255)
                                             tGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACa                                                                                                         <  2:3907471/150‑1 (MQ=255)
                                              gCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAg                                                                                                        >  1:1368502/1‑150 (MQ=255)
                                               ccTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGa                                                                                                       >  2:1169111/1‑150 (MQ=255)
                                               ccTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGa                                                                                                       >  1:374082/1‑150 (MQ=255)
                                                cTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGAt                                                                                                      <  1:2886614/150‑1 (MQ=255)
                                                 tttGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATc                                                                                                     <  1:481771/150‑1 (MQ=255)
                                                 tttGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATc                                                                                                     <  1:4015851/150‑1 (MQ=255)
                                                 tttGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATc                                                                                                     <  1:3202334/150‑1 (MQ=255)
                                                   tGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCa                                                                                                   <  2:1512087/150‑1 (MQ=255)
                                                       ttCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCtt                                                                                               >  1:370812/1‑150 (MQ=255)
                                                       ttCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCtt                                                                                               <  1:356980/150‑1 (MQ=255)
                                                        tCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTg                                                                                              <  2:3351941/150‑1 (MQ=255)
                                                        tCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTg                                                                                              <  1:651086/150‑1 (MQ=255)
                                                         cATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGt                                                                                             <  1:1644568/150‑1 (MQ=255)
                                                          atatTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGtt                                                                                            <  1:1547628/150‑1 (MQ=255)
                                                           tatTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGATGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTa                                                                                           <  2:2822171/150‑1 (MQ=255)
                                                           tatTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTa                                                                                           <  1:2576512/150‑1 (MQ=255)
                                                             ttCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACt                                                                                         <  2:3834028/150‑1 (MQ=255)
                                                                ccACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTc                                                                                      >  2:4384764/1‑150 (MQ=255)
                                                                  aCTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCtt                                                                                    <  1:3045989/150‑1 (MQ=255)
                                                                   cTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCttt                                                                                   >  1:4436717/1‑150 (MQ=255)
                                                                     gACCGAATCCTTTCTGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCTTGACGTGCATAAGTATTATTACCTTCATATGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATACAGCTTGTTACGGTCTttat                                                                                 <  2:3983102/150‑1 (MQ=255)
                                                                        cGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTtattat                                                                              <  2:418137/150‑1 (MQ=255)
                                                                        cGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTtattat                                                                              >  2:262779/1‑150 (MQ=255)
                                                                             ccAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAAtt                                                                         >  1:1530517/1‑150 (MQ=255)
                                                                              cagtcagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAAttt                                                                        <  1:780867/150‑1 (MQ=255)
                                                                               agtcagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTc                                                                       >  2:174632/1‑150 (MQ=255)
                                                                                gtcagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCa                                                                      >  1:3261202/1‑150 (MQ=255)
                                                                                  cagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGc                                                                    >  1:1205892/1‑150 (MQ=255)
                                                                                  cagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGc                                                                    <  2:1002338/150‑1 (MQ=255)
                                                                                  cagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGc                                                                    <  1:387264/150‑1 (MQ=255)
                                                                                   agtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCt                                                                   <  1:1941242/150‑1 (MQ=255)
                                                                                       gATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCt                                                               <  1:3589594/150‑1 (MQ=255)
                                                                                       gATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCt                                                               <  1:2753369/150‑1 (MQ=255)
                                                                                        aTCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTGGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTg                                                              >  1:2037099/1‑150 (MQ=255)
                                                                                        aTCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTg                                                              <  2:759866/150‑1 (MQ=255)
                                                                                          cGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAg                                                            >  1:3367075/1‑150 (MQ=255)
                                                                                           gTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGc                                                           >  2:1950005/1‑150 (MQ=255)
                                                                                            ttGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCa                                                          <  2:3966268/150‑1 (MQ=255)
                                                                                                  tGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATc                                                    <  2:307829/150‑1 (MQ=255)
                                                                                                       ttCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCAt                                               >  1:341973/1‑150 (MQ=255)
                                                                                                         cACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTa                                             >  1:2270819/1‑150 (MQ=255)
                                                                                                           ccTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAATTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGt                                           <  2:2044062/150‑1 (MQ=255)
                                                                                                           ccTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGt                                           <  2:2649979/150‑1 (MQ=255)
                                                                                                            cTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTa                                          >  1:1694825/1‑150 (MQ=255)
                                                                                                             tttGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTAc                                         <  1:312534/150‑1 (MQ=255)
                                                                                                              ttGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACt                                        <  1:4576564/150‑1 (MQ=255)
                                                                                                                   gCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGc                                   >  2:366580/1‑150 (MQ=255)
                                                                                                                    ccAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCa                                  <  2:2818380/150‑1 (MQ=255)
                                                                                                                     cAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAg                                 >  2:444092/1‑150 (MQ=255)
                                                                                                                       aGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCt                               >  1:2638809/1‑150 (MQ=255)
                                                                                                                       aGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCt                               <  2:4084686/150‑1 (MQ=255)
                                                                                                                       aGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCt                               >  1:2837754/1‑150 (MQ=255)
                                                                                                                          cGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTAca                            >  2:687649/1‑150 (MQ=255)
                                                                                                                          cGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTAca                            >  2:3977311/1‑150 (MQ=255)
                                                                                                                           gggCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTAcag                           >  2:42051/1‑150 (MQ=255)
                                                                                                                            ggCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTAcagc                          <  1:4327944/150‑1 (MQ=255)
                                                                                                                               aTAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGa                       <  1:377753/150‑1 (MQ=255)
                                                                                                                                tAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGaa                      >  2:1642108/1‑150 (MQ=255)
                                                                                                                                   gtagtaTCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAAtt                   >  1:1910064/1‑150 (MQ=255)
                                                                                                                                    tagtaTCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTg                  <  2:1368502/150‑1 (MQ=255)
                                                                                                                                         tCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACt             <  2:4350669/150‑1 (MQ=255)
                                                                                                                                            ccATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGtt          <  2:2233698/150‑1 (MQ=255)
                                                                                                                                               tcatcTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAAt       <  1:4269412/150‑1 (MQ=255)
                                                                                                                                               tcatcTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAAt       <  2:686017/150‑1 (MQ=255)
                                                                                                                                                catcTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAAtt      >  2:4514676/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTT  >  CP009273/571981‑572278

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: