Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,386,303 A→T 17.2% I1105I (ATT→ATA yhdP ← DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,386,3030AT17.2% 679.5 / 126.8 249I1105I (ATT→ATAyhdPDUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base A (116/90);  new base T (27/16);  total (143/106)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

TGGAGTCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGTTTTC  >  CP009273/3386155‑3386451
                                                                                                                                                    |                                                                                                                                                    
tGGAGTCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTatgaa                                                                                                                                                     >  1:2428411/1‑150 (MQ=255)
    gTCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAatgat                                                                                                                                                 >  1:1091602/1‑150 (MQ=255)
    gTCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAatgat                                                                                                                                                 <  2:2443600/150‑1 (MQ=255)
       aaaaTAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGcc                                                                                                                                              >  2:811265/1‑150 (MQ=255)
        aaaTAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCa                                                                                                                                             >  2:3431110/1‑150 (MQ=255)
          aTAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCAtt                                                                                                                                           <  1:4324931/150‑1 (MQ=255)
          aTAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCAtt                                                                                                                                           >  1:969502/1‑150 (MQ=255)
            aGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCa                                                                                                                                         <  1:3943652/150‑1 (MQ=255)
            aGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCa                                                                                                                                         >  1:3555535/1‑150 (MQ=255)
             gAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCaa                                                                                                                                        >  1:4080433/1‑150 (MQ=255)
             gAACCCTTCGACAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCaa                                                                                                                                        <  1:3484498/150‑1 (MQ=255)
              aaCCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAAc                                                                                                                                       >  1:2281848/1‑150 (MQ=255)
                    tCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCt                                                                                                                                 <  2:2905848/150‑1 (MQ=255)
                     cGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCtt                                                                                                                                <  2:2218804/150‑1 (MQ=255)
                      gCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTc                                                                                                                               >  2:183392/1‑150 (MQ=255)
                         aaaaGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATc                                                                                                                            >  1:1024530/1‑150 (MQ=255)
                           aaGTGTCTCTGAAATCAAAACGCAGCGTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGAGGCCATTCAACGTCGCTTCATCTg                                                                                                                          <  2:136891/150‑1 (MQ=255)
                            aGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTgg                                                                                                                         >  1:438124/1‑150 (MQ=255)
                             gtgtCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGc                                                                                                                        <  2:413416/150‑1 (MQ=255)
                              tgtCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCt                                                                                                                       >  1:711936/1‑150 (MQ=255)
                               gtCTCTGAAATCAATACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTg                                                                                                                      <  1:3266965/150‑1 (MQ=255)
                               gtCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTg                                                                                                                      <  1:1938944/150‑1 (MQ=255)
                                   ctGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAg                                                                                                                  >  2:2275032/1‑150 (MQ=255)
                                              aCGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTgcg                                                                                                       <  2:3695134/150‑1 (MQ=255)
                                                 cAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCa                                                                                                    >  1:1087651/1‑150 (MQ=255)
                                                  aGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAg                                                                                                   >  1:416333/1‑150 (MQ=255)
                                                  aGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAg                                                                                                   >  2:4263034/1‑150 (MQ=255)
                                                  aGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAg                                                                                                   <  2:4329141/150‑1 (MQ=255)
                                                   gCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAgt                                                                                                  <  2:2418593/150‑1 (MQ=255)
                                                    cTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGGTTTGCGCCAgtg                                                                                                 >  2:1088190/1‑150 (MQ=255)
                                                    cTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAgtg                                                                                                 >  2:4091763/1‑150 (MQ=255)
                                                      tACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTa                                                                                               >  1:1405893/1‑150 (MQ=255)
                                                       aCGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTaa                                                                                              <  2:1216028/150‑1 (MQ=255)
                                                         gaaTCATGGCATCTACGATCAATATGCGCTGCAATTGCCCGGCATGTCATGTATTGTTTTCGGTTATTTCGCCTTTACCCAGTTGAGTATTTATGATGCCATTCAACGTCGCTTATTCTGGCTGGCAGGGTGATTTGCGCCAGTGTatat                                                                                            <  2:700010/148‑1 (MQ=255)
                                                          cATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATc                                                                                           >  1:1131779/1‑150 (MQ=255)
                                                          cATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATc                                                                                           >  1:1131837/1‑150 (MQ=255)
                                                           aTCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCg                                                                                          >  1:4178453/1‑150 (MQ=255)
                                                           aTCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCg                                                                                          >  2:4118694/1‑150 (MQ=255)
                                                              aGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTaa                                                                                       >  1:2075782/1‑150 (MQ=255)
                                                              aGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTaa                                                                                       >  1:306850/1‑150 (MQ=255)
                                                               gggCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAAt                                                                                      <  2:884421/150‑1 (MQ=255)
                                                                ggCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATc                                                                                     <  1:476446/150‑1 (MQ=255)
                                                                   aTCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCcac                                                                                  >  2:4007958/1‑150 (MQ=255)
                                                                   aTCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCcac                                                                                  >  1:1974816/1‑150 (MQ=255)
                                                                   aTCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCcac                                                                                  >  2:1829488/1‑150 (MQ=255)
                                                                    tCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCcaca                                                                                 >  2:345700/1‑150 (MQ=255)
                                                                     cTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACAt                                                                                <  1:1594732/150‑1 (MQ=255)
                                                                     cTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACAt                                                                                <  1:1595109/150‑1 (MQ=255)
                                                                       aCGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTa                                                                              >  2:1468479/1‑150 (MQ=255)
                                                                          cTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAAt                                                                           <  1:2516608/150‑1 (MQ=255)
                                                                           tCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATg                                                                          <  1:527477/150‑1 (MQ=255)
                                                                             aaTAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATgac                                                                        >  1:1374513/1‑150 (MQ=255)
                                                                                aaGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTACTCCACATTAAATgacgac                                                                     >  2:1741408/1‑150 (MQ=255)
                                                                                  gcgcAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTg                                                                   <  1:134779/150‑1 (MQ=255)
                                                                                   cgcAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTgg                                                                  >  1:4035555/1‑150 (MQ=255)
                                                                                     cagcaATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCg                                                                >  2:630338/1‑150 (MQ=255)
                                                                                     cagcaATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCg                                                                >  2:2540036/1‑150 (MQ=255)
                                                                                     cagcaATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCg                                                                >  2:3131292/1‑150 (MQ=255)
                                                                                      agcaATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGt                                                               >  1:1404683/1‑150 (MQ=255)
                                                                                        caATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTAt                                                             >  2:809764/1‑150 (MQ=255)
                                                                                         aaTTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATg                                                            >  1:1497237/1‑150 (MQ=255)
                                                                                         aaTTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATg                                                            >  2:4201435/1‑150 (MQ=255)
                                                                                              cccGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCTCTTCATCATGCTGCCATGGTGCTTTGCGCCAGTGTAAATCGTAAACCACATTAAATGACGACTGGCGTATGGGcga                                                       >  2:1969564/1‑149 (MQ=255)
                                                                                               ccGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGcgtc                                                      <  1:3374943/150‑1 (MQ=255)
                                                                                                cGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGcgtcg                                                     <  2:2605924/150‑1 (MQ=255)
                                                                                                 ggCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGcgtcgt                                                    >  2:2332182/1‑150 (MQ=255)
                                                                                                  gCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTg                                                   <  1:2917992/150‑1 (MQ=255)
                                                                                                    aTGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGac                                                 <  1:539483/150‑1 (MQ=255)
                                                                                                     tGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGaca                                                <  2:3296527/150‑1 (MQ=255)
                                                                                                        ccGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCa                                             <  2:508686/150‑1 (MQ=255)
                                                                                                             aTTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCaaaaaa                                        >  2:1505828/1‑150 (MQ=255)
                                                                                                             aTTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCaaaaaa                                        <  2:921704/150‑1 (MQ=255)
                                                                                                               tGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAAtt                                      <  2:1302541/150‑1 (MQ=255)
                                                                                                               tGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAAtt                                      >  1:739133/1‑150 (MQ=255)
                                                                                                               tGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAAtt                                      >  2:3602731/1‑150 (MQ=255)
                                                                                                                gATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTc                                     >  2:45899/1‑150 (MQ=255)
                                                                                                                 atttatgTATTTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAAGGTCGCTTAATCTGGCTGCCAGGGTGATTTGCGCCAGTGTATATGGTTATACACATTAAATGACGAATGGCGTATGGGCGTCGTGACACCAAAAAATTCt                                    <  1:1741408/144‑1 (MQ=255)
                                                                                                                 aTTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCt                                    <  1:3187785/150‑1 (MQ=255)
                                                                                                                 aTTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCt                                    >  1:1037334/1‑150 (MQ=255)
                                                                                                                  tttCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTg                                   <  2:1554996/150‑1 (MQ=255)
                                                                                                                   ttCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTgc                                  <  1:4206545/150‑1 (MQ=255)
                                                                                                                    tCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTgcg                                 >  1:2129119/1‑150 (MQ=255)
                                                                                                                    tCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTgcg                                 <  2:4032440/150‑1 (MQ=255)
                                                                                                                            ttCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATc                         >  1:810004/1‑150 (MQ=255)
                                                                                                                            ttCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATc                         <  1:2806802/150‑1 (MQ=255)
                                                                                                                             tCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTACATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCa                        >  2:373293/1‑150 (MQ=255)
                                                                                                                               gCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAAt                      <  2:3036648/150‑1 (MQ=255)
                                                                                                                               gCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAAt                      <  1:2238073/150‑1 (MQ=255)
                                                                                                                                ccTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAAtt                     >  2:4347303/1‑150 (MQ=255)
                                                                                                                                ccTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAAtt                     >  1:2670006/1‑150 (MQ=255)
                                                                                                                                   ttACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTc                  <  1:1875466/150‑1 (MQ=255)
                                                                                                                                    tACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCt                 <  1:953034/150‑1 (MQ=255)
                                                                                                                                        cAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCCATTTTCTGccc             >  2:216309/1‑150 (MQ=255)
                                                                                                                                            tGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCgcgc         >  2:957550/1‑150 (MQ=255)
                                                                                                                                              aGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAg       >  1:37239/1‑150 (MQ=255)
                                                                                                                                              aGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAg       >  1:2211836/1‑150 (MQ=255)
                                                                                                                                               gTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGt      <  2:2749642/150‑1 (MQ=255)
                                                                                                                                                  tgaatgaTGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGtttt   <  1:3336467/150‑1 (MQ=255)
                                                                                                                                                  tgaatgaTGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCGTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGgttt   >  1:389832/1‑150 (MQ=255)
                                                                                                                                                   gaatgaTGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGTTTTc  <  1:1396451/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGAGTCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGTTTTC  >  CP009273/3386155‑3386451

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: