Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,529,263 A→T 12.8% 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,2630AT12.8% 884.4 / 94.8 285L19* (TTA→TAA) nanCN‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base A (114/134);  new base T (11/26);  total (125/160)
Fisher's exact test for biased strand distribution p-value = 7.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

TATTATCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATA  >  CP009273/4529114‑4529412
                                                                                                                                                     |                                                                                                                                                     
    atCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATgg                                                                                                                                                   >  2:5121332/1‑150 (MQ=255)
     tCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATggg                                                                                                                                                  <  1:1063642/150‑1 (MQ=255)
      cATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATgggg                                                                                                                                                 >  2:4055431/1‑150 (MQ=255)
       aTGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGACATTAATGGGGa                                                                                                                                                >  1:4480174/1‑150 (MQ=255)
       aTGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGa                                                                                                                                                >  2:1131162/1‑150 (MQ=255)
         gAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACg                                                                                                                                              >  1:3382712/1‑150 (MQ=255)
          aaTGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGa                                                                                                                                             >  2:1554839/1‑150 (MQ=255)
           aTGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGaa                                                                                                                                            >  2:2658577/1‑150 (MQ=255)
            tGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGaaa                                                                                                                                           <  1:4974015/150‑1 (MQ=255)
            tGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGaaa                                                                                                                                           <  1:151197/150‑1 (MQ=255)
             ggTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGaaaa                                                                                                                                          <  2:3418590/150‑1 (MQ=255)
                aTTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGca                                                                                                                                       >  2:3038438/1‑150 (MQ=255)
                  tCCAGGTATTACTTTCCATGCTTGCCCAACATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCTTAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGcaca                                                                                                                                     <  1:4542539/150‑1 (MQ=255)
                   ccAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAg                                                                                                                                    >  1:3929554/1‑150 (MQ=255)
                   ccAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAg                                                                                                                                    >  1:200624/1‑150 (MQ=255)
                    cAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGt                                                                                                                                   <  2:3954751/150‑1 (MQ=255)
                     aGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGta                                                                                                                                  <  1:581330/150‑1 (MQ=255)
                     aGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGta                                                                                                                                  <  2:632105/150‑1 (MQ=255)
                     aGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGta                                                                                                                                  >  1:2701262/1‑150 (MQ=255)
                       gTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaat                                                                                                                                >  1:1306847/1‑150 (MQ=255)
                        tATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaata                                                                                                                               <  1:2559382/150‑1 (MQ=255)
                           tACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAataata                                                                                                                            >  1:4403642/1‑150 (MQ=255)
                           tACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAataata                                                                                                                            >  2:1134628/1‑150 (MQ=255)
                              tttCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGc                                                                                                                         <  2:492350/150‑1 (MQ=255)
                                tCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCa                                                                                                                       >  2:4563875/1‑150 (MQ=255)
                                tCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCa                                                                                                                       <  1:2261678/150‑1 (MQ=255)
                                   aTGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGaa                                                                                                                    >  2:1347632/1‑150 (MQ=255)
                                   aTGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGaa                                                                                                                    >  2:4946685/1‑150 (MQ=255)
                                    tGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGaaa                                                                                                                   >  1:1401146/1‑150 (MQ=255)
                                        tGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAt                                                                                                               >  2:3579042/1‑150 (MQ=255)
                                         gCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAtt                                                                                                              <  2:887975/150‑1 (MQ=255)
                                         gCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAtt                                                                                                              <  1:4160022/150‑1 (MQ=255)
                                             accaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAg                                                                                                          >  2:1179667/1‑150 (MQ=255)
                                               catccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGcc                                                                                                        >  1:3763072/1‑150 (MQ=255)
                                                atccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCt                                                                                                       <  1:1172757/150‑1 (MQ=255)
                                                 tccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtt                                                                                                      <  2:2114713/150‑1 (MQ=255)
                                                 tccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtt                                                                                                      <  1:1617124/150‑1 (MQ=255)
                                                  ccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCttt                                                                                                     >  2:72969/1‑150 (MQ=255)
                                                   catTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCttnt                                                                                                    <  1:4950154/150‑1 (MQ=255)
                                                    atTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCttttt                                                                                                   <  1:3133827/150‑1 (MQ=255)
                                                     ttttGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTc                                                                                                  <  1:4467379/150‑1 (MQ=255)
                                                       ttGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCAt                                                                                                <  1:1702972/150‑1 (MQ=255)
                                                          ccATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTTCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATaaa                                                                                             <  1:3278620/150‑1 (MQ=255)
                                                          ccATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATaaa                                                                                             <  2:3189683/150‑1 (MQ=255)
                                                           cATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAt                                                                                            <  2:2827988/150‑1 (MQ=255)
                                                           cATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAt                                                                                            >  2:771950/1‑150 (MQ=255)
                                                             tCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAttt                                                                                          <  1:5081971/150‑1 (MQ=255)
                                                             tCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAttt                                                                                          <  2:4275694/150‑1 (MQ=255)
                                                              cccTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTc                                                                                         <  2:355304/150‑1 (MQ=255)
                                                                 tCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACt                                                                                      <  1:1315867/150‑1 (MQ=255)
                                                                  cACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTc                                                                                     <  1:3879506/150‑1 (MQ=255)
                                                                  cACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTc                                                                                     <  2:1961330/150‑1 (MQ=255)
                                                                  cACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTc                                                                                     <  2:2492482/150‑1 (MQ=255)
                                                                    cTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAt                                                                                   <  2:180510/150‑1 (MQ=255)
                                                                     tGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAtt                                                                                  >  1:133134/1‑150 (MQ=255)
                                                                     tGACTTTGAGCCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAtt                                                                                  <  2:1897366/150‑1 (MQ=255)
                                                                      gACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAttt                                                                                 <  1:2388824/150‑1 (MQ=255)
                                                                            gagtcgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGa                                                                           >  2:1216144/1‑150 (MQ=255)
                                                                                  agtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGa                                                                     <  1:2847336/150‑1 (MQ=255)
                                                                                     ctcATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAg                                                                  <  2:3617089/150‑1 (MQ=255)
                                                                                      tcATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGc                                                                 <  1:4578344/150‑1 (MQ=255)
                                                                                          aGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATa                                                             <  2:4987841/150‑1 (MQ=255)
                                                                                           ggCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATa                                                                         <  1:2256699/137‑1 (MQ=255)
                                                                                            gCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACa                                                           <  2:3141803/150‑1 (MQ=255)
                                                                                               tGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAg                                                        <  1:4222214/150‑1 (MQ=255)
                                                                                                 gCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGcc                                                      <  2:2277171/150‑1 (MQ=255)
                                                                                                  cTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGccc                                                     >  1:4891753/1‑150 (MQ=255)
                                                                                                     ccACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCa                                                  >  2:2357232/1‑150 (MQ=255)
                                                                                                       aCTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATa                                                <  2:2251139/150‑1 (MQ=255)
                                                                                                       aCTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATa                                                >  1:2471149/1‑150 (MQ=255)
                                                                                                       aCTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATa                                                >  2:2435705/1‑150 (MQ=255)
                                                                                                        ctccGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGTCGAAAATCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATATATTTCACTCTTTTGTTGGATACAGAAAGCAATACAAAGCCCGCATaa                                               <  2:2193106/147‑1 (MQ=255)
                                                                                                         tACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATaaa                                              <  1:3698733/150‑1 (MQ=255)
                                                                                                         tACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATGAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATaaa                                              >  2:4450011/1‑150 (MQ=255)
                                                                                                          aCGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAAc                                             <  1:85298/150‑1 (MQ=255)
                                                                                                          aCGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAAc                                             <  1:3490220/150‑1 (MQ=255)
                                                                                                           cGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGGAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACa                                            <  2:16901/150‑1 (MQ=255)
                                                                                                              tatCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAAtt                                         <  2:4627256/150‑1 (MQ=255)
                                                                                                                      acgtacgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTa                                 >  1:1859950/1‑150 (MQ=255)
                                                                                                                        gtacgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtg                               >  1:665439/1‑150 (MQ=255)
                                                                                                                          acgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgtt                             <  2:4357231/150‑1 (MQ=255)
                                                                                                                            gtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgttgt                           <  1:981123/150‑1 (MQ=255)
                                                                                                                            gtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgttgt                           <  1:3552814/150‑1 (MQ=255)
                                                                                                                             tCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTtgtg                          <  1:195720/150‑1 (MQ=255)
                                                                                                                             tCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTtgtg                          >  2:4910830/1‑150 (MQ=255)
                                                                                                                                  tgtCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAata                     <  2:200624/150‑1 (MQ=255)
                                                                                                                                   gtCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAatat                    <  1:1131162/150‑1 (MQ=255)
                                                                                                                                    tCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAtt                   >  1:2522335/1‑150 (MQ=255)
                                                                                                                                       cAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAttttt                <  2:1401146/150‑1 (MQ=255)
                                                                                                                                        aGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAtttttt               >  1:5138948/1‑150 (MQ=255)
                                                                                                                                          ccTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTg             <  1:3980348/150‑1 (MQ=255)
                                                                                                                                           cTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGc            >  1:4717126/1‑150 (MQ=255)
                                                                                                                                             gagaAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCa          >  2:5083612/1‑150 (MQ=255)
                                                                                                                                              agaAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAg         <  1:3595970/150‑1 (MQ=255)
                                                                                                                                              agaAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAATTATTTTAGCCTTTTTCATAAATTTAACTCATTTTTAGTATACAGAAAGAAATACAAATCCAGCATAAACAATTATCATTTATGTTGTGTAATATTTTTTTGACAg         >  1:1565552/1‑150 (MQ=255)
                                                                                                                                                aaaTTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGc       <  1:4512841/150‑1 (MQ=255)
                                                                                                                                                   ttAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTta    <  1:465540/150‑1 (MQ=255)
                                                                                                                                                     aaTGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTtata  >  2:4733279/1‑150 (MQ=255)
                                                                                                                                                     aaTGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTtata  >  1:260342/1‑150 (MQ=255)
                                                                                                                                                     aaTGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTtata  <  1:72969/150‑1 (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: