Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,529,263 A→T 9.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,2630AT 9.2% 722.8 / 49.1 218L19* (TTA→TAA) nanCN‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base A (70/128);  new base T (7/13);  total (77/141)
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.85e-01

Only 100 of 220 total aligned reads displayed.

TATTATCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTAT  >  CP009273/4529114‑4529411
                                                                                                                                                     |                                                                                                                                                    
tattatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTa                                                                                                                                                      >  1:399123/1‑150 (MQ=255)
tattatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTa                                                                                                                                                      >  2:777178/1‑150 (MQ=255)
  ttatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAAt                                                                                                                                                    <  2:1691320/150‑1 (MQ=255)
   tatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATg                                                                                                                                                   <  1:2515821/150‑1 (MQ=255)
    atCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATgg                                                                                                                                                  >  1:2812665/1‑150 (MQ=255)
      cATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATgggg                                                                                                                                                <  2:4728913/150‑1 (MQ=255)
       aTGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGa                                                                                                                                               >  1:957949/1‑150 (MQ=255)
        tGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGAc                                                                                                                                              >  1:3419526/1‑150 (MQ=255)
         gAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACg                                                                                                                                             >  2:513816/1‑150 (MQ=255)
         gAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACg                                                                                                                                             >  2:2035424/1‑150 (MQ=255)
         gAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACg                                                                                                                                             <  2:2198939/150‑1 (MQ=255)
          aaTGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGa                                                                                                                                            >  2:4191424/1‑150 (MQ=255)
           aTGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGaa                                                                                                                                           >  1:569170/1‑150 (MQ=255)
           aTGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGaa                                                                                                                                           <  1:3676738/150‑1 (MQ=255)
           aTGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACAGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCAACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATTGGGACGaa                                                                                                                                           >  1:2958200/1‑150 (MQ=255)
            tGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGaac                                                                                                                                          >  2:1192717/1‑149 (MQ=255)
             ggTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGACATTAATGGGGACGaaaa                                                                                                                                         >  2:1726201/1‑150 (MQ=255)
                   ccAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAg                                                                                                                                   >  2:3507949/1‑150 (MQ=255)
                   ccAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAg                                                                                                                                   >  2:1660758/1‑150 (MQ=255)
                     aGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGta                                                                                                                                 >  1:3632069/1‑150 (MQ=255)
                      ggTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaa                                                                                                                                >  2:1255019/1‑150 (MQ=255)
                      ggTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaa                                                                                                                                <  2:3999474/150‑1 (MQ=255)
                       gTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaat                                                                                                                               <  2:4120614/150‑1 (MQ=255)
                        tATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaata                                                                                                                              <  1:4211881/150‑1 (MQ=255)
                        tATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaata                                                                                                                              <  1:3297808/150‑1 (MQ=255)
                        tATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaata                                                                                                                              <  2:3599028/150‑1 (MQ=255)
                         aTTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaataa                                                                                                                             >  1:1174836/1‑150 (MQ=255)
                         aTTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaataa                                                                                                                             <  2:2126069/150‑1 (MQ=255)
                            aCTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATAc                                                                                                                          <  1:2227567/150‑1 (MQ=255)
                                tCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCa                                                                                                                      >  2:2445041/1‑150 (MQ=255)
                                     gCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAg                                                                                                                 <  2:317861/150‑1 (MQ=255)
                                       ttGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTa                                                                                                               <  1:3150254/150‑1 (MQ=255)
                                       ttGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTa                                                                                                               >  1:3393678/1‑150 (MQ=255)
                                       ttGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTa                                                                                                               >  1:1719277/1‑150 (MQ=255)
                                        tGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAt                                                                                                              >  2:4292387/1‑150 (MQ=255)
                                        tGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAt                                                                                                              <  1:1599416/150‑1 (MQ=255)
                                          cccACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAttt                                                                                                            <  2:1819245/150‑1 (MQ=255)
                                            caccaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTa                                                                                                          <  1:1087739/150‑1 (MQ=255)
                                            caccaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTa                                                                                                          <  2:1920672/150‑1 (MQ=255)
                                            caccaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTa                                                                                                          <  1:1907841/150‑1 (MQ=255)
                                                   catTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtttt                                                                                                   <  2:3845289/150‑1 (MQ=255)
                                                   catTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtttt                                                                                                   <  2:2590194/150‑1 (MQ=255)
                                                   catTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtttt                                                                                                   <  1:148214/150‑1 (MQ=255)
                                                    atTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCttttt                                                                                                  <  2:3615097/150‑1 (MQ=255)
                                                    atTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCttttt                                                                                                  <  2:3948008/150‑1 (MQ=255)
                                                      tttGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCa                                                                                                >  1:4388633/1‑150 (MQ=255)
                                                      tttGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCa                                                                                                >  1:3529417/1‑150 (MQ=255)
                                                         gCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATaa                                                                                             <  2:2105138/150‑1 (MQ=255)
                                                         gCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATaa                                                                                             <  1:3430522/150‑1 (MQ=255)
                                                           cATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAt                                                                                           <  1:1090561/150‑1 (MQ=255)
                                                            aTCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAtt                                                                                          >  1:2747465/1‑150 (MQ=255)
                                                             tCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAttt                                                                                         <  1:2625574/150‑1 (MQ=255)
                                                                 tCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACt                                                                                     >  2:1155914/1‑150 (MQ=255)
                                                                  cACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTc                                                                                    >  2:1538775/1‑150 (MQ=255)
                                                                      gACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAttt                                                                                <  1:1450956/150‑1 (MQ=255)
                                                                       aCTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTg                                                                               >  2:862743/1‑150 (MQ=255)
                                                                       aCTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTg                                                                               <  1:4575466/150‑1 (MQ=255)
                                                                             agtcgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAGTTTCACTCAGTTGTAGGAt                                                                         >  1:2200476/1‑150 (MQ=255)
                                                                               tcgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATAc                                                                       <  2:1035975/150‑1 (MQ=255)
                                                                                 gagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAg                                                                     >  1:300795/1‑150 (MQ=255)
                                                                                 gagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAg                                                                     <  2:3832828/150‑1 (MQ=255)
                                                                                   gtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGaa                                                                   <  1:3266397/150‑1 (MQ=255)
                                                                                   gtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGaa                                                                   >  1:4184201/1‑150 (MQ=255)
                                                                                     ctcATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAg                                                                 <  1:3625094/150‑1 (MQ=255)
                                                                                      tcATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGc                                                                <  2:1629134/150‑1 (MQ=255)
                                                                                      tcATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGc                                                                <  1:4557255/150‑1 (MQ=255)
                                                                                          aGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATa                                                            >  2:1437859/1‑150 (MQ=255)
                                                                                               tGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAg                                                       <  1:2512888/150‑1 (MQ=255)
                                                                                                  cTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGccc                                                    >  1:4457983/1‑150 (MQ=255)
                                                                                                  cTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGccc                                                    <  2:2687888/150‑1 (MQ=255)
                                                                                                  cTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGccc                                                    <  2:2687659/150‑1 (MQ=255)
                                                                                                   ttCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCg                                                   <  2:2196831/150‑1 (MQ=255)
                                                                                                      cACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCAt                                                <  2:3419552/150‑1 (MQ=255)
                                                                                                             atatCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAAt                                         <  2:3073204/150‑1 (MQ=255)
                                                                                                                tCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAg                                      <  1:1255019/150‑1 (MQ=255)
                                                                                                                 ccaccaCGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGc                                     <  1:1371554/150‑1 (MQ=255)
                                                                                                                 ccaccaCGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGc                                     <  1:4074810/150‑1 (MQ=255)
                                                                                                                   accacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCAt                                   >  1:4417107/1‑150 (MQ=255)
                                                                                                                    ccacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCAtt                                  >  2:1753060/1‑150 (MQ=255)
                                                                                                                      acgtacgtCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTa                                <  1:1151925/150‑1 (MQ=255)
                                                                                                                        gtacgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtg                              >  1:3784403/1‑150 (MQ=255)
                                                                                                                          acgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgtt                            >  1:2458095/1‑150 (MQ=255)
                                                                                                                          acgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgtt                            >  2:908811/1‑150 (MQ=255)
                                                                                                                          acgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgtt                            <  2:1759179/150‑1 (MQ=255)
                                                                                                                            gtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgttgt                          >  1:468146/1‑150 (MQ=255)
                                                                                                                              ccAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACACTTAGCAGTTATGTTgtgt                        >  2:3702502/1‑150 (MQ=255)
                                                                                                                               cAGTGTCGCATCATTAGAAATTAATGTGGAAGAAAATCACATTAATAATAAGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTa                       <  2:2820786/150‑1 (MQ=255)
                                                                                                                               cAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTa                       <  1:225206/150‑1 (MQ=255)
                                                                                                                               cAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTa                       <  2:1588074/150‑1 (MQ=255)
                                                                                                                               cAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTa                       <  2:1588688/150‑1 (MQ=255)
                                                                                                                               cAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTa                       <  1:4028956/150‑1 (MQ=255)
                                                                                                                                aGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTaa                      <  1:777178/150‑1 (MQ=255)
                                                                                                                                 gtgtCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAat                     <  1:1329022/150‑1 (MQ=255)
                                                                                                                                   gtCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCGTTTATGTTGTGTAatat                   <  2:3233751/150‑1 (MQ=255)
                                                                                                                                   gtCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAatat                   <  2:280336/150‑1 (MQ=255)
                                                                                                                                         gCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAttttttt             <  1:2442972/150‑1 (MQ=255)
                                                                                                                                               gaAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAgg       >  2:2920063/1‑150 (MQ=255)
                                                                                                                                                 aaTTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCt     <  2:4079424/150‑1 (MQ=255)
                                                                                                                                                  aTTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAATCAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCtt    >  2:3450623/1‑150 (MQ=255)
                                                                                                                                                  aTTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCtt    >  2:3857239/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATTATCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTAT  >  CP009273/4529114‑4529411

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: