Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,529,263 A→T 17.2% L19* (TTA→TAA)  nanC ← N‑acetylnuraminic acid outer membrane channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,529,2630AT17.2% 594.3 / 109.3 221L19* (TTA→TAA) nanCN‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base A (65/118);  new base T (20/18);  total (85/136)
Fisher's exact test for biased strand distribution p-value = 6.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 227 total aligned reads displayed.

TATTATCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATA  >  CP009273/4529114‑4529412
                                                                                                                                                     |                                                                                                                                                     
tattatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTa                                                                                                                                                       <  2:34601/150‑1 (MQ=255)
tattatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTa                                                                                                                                                       >  2:3416973/1‑150 (MQ=255)
  ttatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTAt                                                                                                                                                     >  2:2493142/1‑150 (MQ=255)
  ttatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAAt                                                                                                                                                     <  2:3276418/150‑1 (MQ=255)
    atCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATgg                                                                                                                                                   <  2:352750/150‑1 (MQ=255)
     tCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATggg                                                                                                                                                  <  2:2577134/150‑1 (MQ=255)
     tCATGAATGGTATTCCAGGTATTACTTTCCATGCGTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATggg                                                                                                                                                  <  2:401359/150‑1 (MQ=255)
     tCATGAATGGTATTCCAGGTATTACTGTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATggg                                                                                                                                                  <  2:879178/150‑1 (MQ=255)
        tGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGAc                                                                                                                                               <  2:1664746/150‑1 (MQ=255)
         gAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACg                                                                                                                                              <  1:3578298/150‑1 (MQ=255)
         gAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACg                                                                                                                                              >  1:576821/1‑150 (MQ=255)
             ggTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGaaaa                                                                                                                                          >  1:4353487/1‑150 (MQ=255)
               tATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGc                                                                                                                                        >  2:3359312/1‑150 (MQ=255)
                aTTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGca                                                                                                                                       >  2:3741398/1‑150 (MQ=255)
                    cAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGt                                                                                                                                   >  2:203395/1‑150 (MQ=255)
                    cAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGt                                                                                                                                   >  1:1037440/1‑150 (MQ=255)
                     aGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGta                                                                                                                                  <  2:2244874/150‑1 (MQ=255)
                        tATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaata                                                                                                                               <  1:2757808/150‑1 (MQ=255)
                         aTTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaataa                                                                                                                              >  1:1447331/1‑150 (MQ=255)
                         aTTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaataa                                                                                                                              >  2:3600406/1‑150 (MQ=255)
                         aTTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaataa                                                                                                                              <  1:348527/150‑1 (MQ=255)
                         aTTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaataa                                                                                                                              <  1:284754/150‑1 (MQ=255)
                          ttACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTaataat                                                                                                                             >  2:3797358/1‑150 (MQ=255)
                          ttACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTaataat                                                                                                                             >  1:1630954/1‑150 (MQ=255)
                            aCTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATAc                                                                                                                           <  1:454723/150‑1 (MQ=255)
                               ttCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGcc                                                                                                                        >  1:3206135/1‑150 (MQ=255)
                               ttCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGcc                                                                                                                        <  2:1547253/150‑1 (MQ=255)
                                tCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCa                                                                                                                       <  1:2579294/150‑1 (MQ=255)
                                  cATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGa                                                                                                                     <  1:3567770/150‑1 (MQ=255)
                                   aTGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGaa                                                                                                                    >  1:1263304/1‑150 (MQ=255)
                                   aTGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGaa                                                                                                                    >  2:4347965/1‑150 (MQ=255)
                                     gCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAg                                                                                                                  <  1:159336/150‑1 (MQ=255)
                                      cTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGt                                                                                                                 >  1:528313/1‑150 (MQ=255)
                                      cTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGt                                                                                                                 >  2:3564014/1‑150 (MQ=255)
                                       ttGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTa                                                                                                                <  2:215476/150‑1 (MQ=255)
                                        tGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAt                                                                                                               >  2:2284624/1‑150 (MQ=255)
                                        tGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAt                                                                                                               >  1:1291144/1‑150 (MQ=255)
                                          cccACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAttt                                                                                                             >  1:1837054/1‑150 (MQ=255)
                                            caccaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTa                                                                                                           >  1:1579671/1‑150 (MQ=255)
                                            caccaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTa                                                                                                           >  1:3768097/1‑150 (MQ=255)
                                            caccaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTa                                                                                                           >  1:1049414/1‑150 (MQ=255)
                                             accaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAg                                                                                                          <  1:4351116/150‑1 (MQ=255)
                                                atccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATGTTAGCCt                                                                                                       >  1:2222172/1‑150 (MQ=255)
                                                 tccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtt                                                                                                      >  1:1183920/1‑150 (MQ=255)
                                                   catTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtttt                                                                                                    <  2:4291622/150‑1 (MQ=255)
                                                     ttttGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAt                                                                                                               <  2:705010/137‑1 (MQ=255)
                                                     ttttGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAt                                                                                                               >  1:705010/1‑137 (MQ=255)
                                                     ttgtGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTc                                                                                                  <  2:963257/147‑1 (MQ=255)
                                                          ccATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATaaa                                                                                             >  1:470300/1‑150 (MQ=255)
                                                           cATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAt                                                                                            >  1:2509102/1‑150 (MQ=255)
                                                               ccTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCa                                                                                        <  1:4044072/150‑1 (MQ=255)
                                                                 tCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACt                                                                                      >  2:1851647/1‑150 (MQ=255)
                                                                 tCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACt                                                                                      <  2:1039763/150‑1 (MQ=255)
                                                                  cACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTc                                                                                     >  2:1910600/1‑150 (MQ=255)
                                                                      gACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAttt                                                                                 >  1:2129747/1‑150 (MQ=255)
                                                                         tttGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGAGATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTa                                                                              <  2:3069698/150‑1 (MQ=255)
                                                                           tGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCCCGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAgg                                                                            <  2:1126556/150‑1 (MQ=255)
                                                                             agtcgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGAt                                                                          <  2:4049795/150‑1 (MQ=255)
                                                                                cgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACa                                                                       >  1:1072486/1‑150 (MQ=255)
                                                                                     ctcATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAg                                                                  >  1:1677825/1‑150 (MQ=255)
                                                                                        aTAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCaa                                                               >  1:2914259/1‑150 (MQ=255)
                                                                                           ggCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATAc                                                            >  1:3050070/1‑150 (MQ=255)
                                                                                              gTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACaaa                                                         <  1:4007103/150‑1 (MQ=255)
                                                                                              gTGGCTTACACTATGATATCCACCAAGTACGTCCAGTGTCGCTGCCTGAGAAATTAATGGTGACGAAAAGCAAAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCTTTTATTTCACTCATTTGTAGGATACAGAAAGAAATACaaa                                                         <  1:4308236/150‑1 (MQ=255)
                                                                                                ggCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGc                                                       <  1:942909/150‑1 (MQ=255)
                                                                                                  cTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGccc                                                     <  1:2010919/150‑1 (MQ=255)
                                                                                                  cTTCAACTACGATATCCACAACGTACGTGCAGTGTCGCAGCCTGATAAAATAATGGGGACGAATATCACAGTTATATTACTCAAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATTCAGTAAGCATTAAAAAGccc                                                     <  2:3151918/150‑1 (MQ=255)
                                                                                                   ttCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCg                                                    >  2:2217048/1‑150 (MQ=255)
                                                                                                   ttCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCg                                                    >  1:2456912/1‑150 (MQ=255)
                                                                                                     ccACTACGATATCCACCACGTACTTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCa                                                  <  2:573680/150‑1 (MQ=255)
                                                                                                     ccACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAGTTGTAGGATACAGAAAGCAATACAAAGCCCGCa                                                  >  1:2225687/1‑150 (MQ=255)
                                                                                                      cACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCAt                                                 <  2:1881361/150‑1 (MQ=255)
                                                                                                      cACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCAt                                                 <  1:300705/150‑1 (MQ=255)
                                                                                                      cACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCAt                                                 <  2:1226419/150‑1 (MQ=255)
                                                                                                      cACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCAt                                                 <  1:299721/150‑1 (MQ=255)
                                                                                                      cACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCAt                                                 >  1:4135691/1‑150 (MQ=255)
                                                                                                                tCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAg                                       <  2:2969367/150‑1 (MQ=255)
                                                                                                                tCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAg                                       <  1:3460610/150‑1 (MQ=255)
                                                                                                                   accacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCAt                                    <  2:470300/150‑1 (MQ=255)
                                                                                                                   accacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCAt                                    <  1:4247489/150‑1 (MQ=255)
                                                                                                                   accacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCAt                                    <  1:2011573/150‑1 (MQ=255)
                                                                                                                    ccacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCAtt                                   >  2:2828227/1‑150 (MQ=255)
                                                                                                                     cacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCAttt                                  <  2:4139318/150‑1 (MQ=255)
                                                                                                                         tacgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgt                              >  1:1012349/1‑150 (MQ=255)
                                                                                                                         tacgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAGTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgt                              >  1:2379032/1‑150 (MQ=255)
                                                                                                                              ccAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTgtgt                         <  1:4347965/150‑1 (MQ=255)
                                                                                                                                   gtCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAatat                    <  1:336857/150‑1 (MQ=255)
                                                                                                                                      gCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAtttt                 <  1:3793121/150‑1 (MQ=255)
                                                                                                                                         gCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAttttttt              <  2:2743578/150‑1 (MQ=255)
                                                                                                                                          ccTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTg             <  2:989260/150‑1 (MQ=255)
                                                                                                                                            tGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGcc           <  2:2509102/150‑1 (MQ=255)
                                                                                                                                              agaAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAg         <  1:3979031/150‑1 (MQ=255)
                                                                                                                                              agaAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAg         <  1:1127143/150‑1 (MQ=255)
                                                                                                                                              agaAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTCCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAg         <  1:2724829/150‑1 (MQ=255)
                                                                                                                                               gaAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAgg        <  1:3026714/150‑1 (MQ=255)
                                                                                                                                                 aaTTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCATTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACTAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCt      <  1:1944857/150‑1 (MQ=255)
                                                                                                                                                   ttAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTta    <  1:1475028/150‑1 (MQ=255)
                                                                                                                                                    tAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTtat   <  2:3640383/150‑1 (MQ=255)
                                                                                                                                                    tAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTtat   <  1:755976/150‑1 (MQ=255)
                                                                                                                                                     aaTGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTtata  >  2:1584701/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTATCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATA  >  CP009273/4529114‑4529412

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: