Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,186,708 G→T 9.5% T245T (ACC→ACA)  purB ← adenylosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,186,7080GT 9.5% 724.2 / 31.9 205T245T (ACC→ACA) purBadenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base G (92/93);  new base T (8/12);  total (100/105)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

ACCAATCTCACCAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCACCTGGTTAAGCTGG  >  CP009273/1186561‑1186855
                                                                                                                                                   |                                                                                                                                                   
aCCAATCTCACCAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGggtg                                                                                                                                                     <  1:2549586/148‑1 (MQ=255)
 ccAATCTCACCAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGggtgg                                                                                                                                                    <  1:810354/148‑1 (MQ=255)
   aaTCTCACCAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGgtggtg                                                                                                                                                  <  1:504464/148‑1 (MQ=255)
    aTCTCACCAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGgtgt                                                                                                                                                 <  1:2201104/148‑1 (MQ=255)
    aTCTCACCAGCAATGGTTTTCGGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGgtgt                                                                                                                                                 >  2:1497315/1‑148 (MQ=255)
      ctcACCAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTAc                                                                                                                                               <  1:1003174/148‑1 (MQ=255)
      ctcACCAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTAc                                                                                                                                               <  1:3200504/148‑1 (MQ=255)
        cACCAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACgg                                                                                                                                             >  1:1779818/1‑148 (MQ=255)
           cAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGtt                                                                                                                                          <  1:1594076/148‑1 (MQ=255)
            aGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTc                                                                                                                                         <  2:2055205/148‑1 (MQ=255)
              cAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTTGTGTACGGGTTCCa                                                                                                                                       <  1:2868852/148‑1 (MQ=255)
              cAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTTGTGTACGGGTTCCa                                                                                                                                       <  1:79869/148‑1 (MQ=255)
                aTGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTTGTGTACGGGTTCCACt                                                                                                                                     >  1:2191181/1‑148 (MQ=255)
                aTGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACt                                                                                                                                     >  2:2556858/1‑148 (MQ=255)
                   gTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGaa                                                                                                                                  <  1:3443226/148‑1 (MQ=255)
                       tCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTTGTGTACGGGTTCCACTGAATAcc                                                                                                                              <  1:583386/148‑1 (MQ=255)
                           tttGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGc                                                                                                                          >  1:2486550/1‑148 (MQ=255)
                           tttGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGc                                                                                                                          >  1:2486255/1‑148 (MQ=255)
                             tGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGa                                                                                                                        <  2:1921390/148‑1 (MQ=255)
                                  tGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGa                                                                                                                   <  2:398591/148‑1 (MQ=255)
                                  tGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGa                                                                                                                   >  2:1594042/1‑148 (MQ=255)
                                     ttAAGGGCGATATAACCCCAGACGTCACTGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGa                                                                                                                <  1:2095520/148‑1 (MQ=255)
                                     ttAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGa                                                                                                                >  1:2011486/1‑148 (MQ=255)
                                     ttAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGa                                                                                                                <  2:2080429/148‑1 (MQ=255)
                                         gggCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAActc                                                                                                            >  2:322879/1‑148 (MQ=255)
                                         gggCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAActc                                                                                                            >  1:12303/1‑148 (MQ=255)
                                         gggCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAActc                                                                                                            >  1:2202844/1‑148 (MQ=255)
                                         gggCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATACTGTTGAAGCGCGCAACCCAATCAAACAGGTCGGCAATCTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATAACAAGCGAGGTGACGAActc                                                                                                            >  1:3010905/1‑148 (MQ=255)
                                            cGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTc                                                                                                         >  2:3170202/1‑148 (MQ=255)
                                             gATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCg                                                                                                        <  2:1131035/148‑1 (MQ=255)
                                                ataACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTg                                                                                                     >  2:690536/1‑148 (MQ=255)
                                                  aaCCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGaa                                                                                                   <  2:3183389/148‑1 (MQ=255)
                                                  aaCCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGaa                                                                                                   <  1:2351103/148‑1 (MQ=255)
                                                    ccccAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACt                                                                                                 >  2:1781475/1‑148 (MQ=255)
                                                     cccAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTTGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTg                                                                                                <  2:2817747/148‑1 (MQ=255)
                                                     cccAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTg                                                                                                >  1:177946/1‑148 (MQ=255)
                                                              cACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTTGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTc                                                                                       >  2:2005214/1‑148 (MQ=255)
                                                              cACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTc                                                                                       >  2:3422090/1‑148 (MQ=255)
                                                                cGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAAATGATGCCAGTCaa                                                                                     >  2:956189/1‑148 (MQ=255)
                                                                   tCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTTGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACtt                                                                                  <  2:742684/148‑1 (MQ=255)
                                                                   tCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACtt                                                                                  >  2:3436315/1‑148 (MQ=255)
                                                                   tCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACtt                                                                                  >  2:2430142/1‑148 (MQ=255)
                                                                    cAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTc                                                                                 <  1:3087612/148‑1 (MQ=255)
                                                                    cAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTc                                                                                 <  1:1202668/148‑1 (MQ=255)
                                                                       aGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCgg                                                                              >  2:797260/1‑148 (MQ=255)
                                                                            aTCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAg                                                                         >  1:1680222/1‑148 (MQ=255)
                                                                            aTCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAg                                                                         >  2:3066740/1‑148 (MQ=255)
                                                                              cAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCACCTTCCGGGTAAgcg                                                                       >  2:3186610/1‑148 (MQ=255)
                                                                                gAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAgcggc                                                                     >  1:1239767/1‑148 (MQ=255)
                                                                                 aaTAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAgcggcg                                                                    <  2:1194266/148‑1 (MQ=255)
                                                                                  aTAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGa                                                                   <  2:1376515/148‑1 (MQ=255)
                                                                                    aGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTTGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGAtg                                                                 >  1:390011/1‑148 (MQ=255)
                                                                                       gtTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTgg                                                              <  1:1075139/148‑1 (MQ=255)
                                                                                       gtTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTgg                                                              >  1:2627206/1‑148 (MQ=255)
                                                                                         tGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGc                                                            >  1:2892347/1‑148 (MQ=255)
                                                                                          gAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCg                                                           <  1:227992/148‑1 (MQ=255)
                                                                                            aGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGtt                                                         <  2:301909/148‑1 (MQ=255)
                                                                                              cgcgcAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTtat                                                       <  2:2831025/148‑1 (MQ=255)
                                                                                              cgcgcAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTtat                                                       <  2:614269/148‑1 (MQ=255)
                                                                                                cgcaacgcaaTCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAg                                                     <  1:1742443/148‑1 (MQ=255)
                                                                                                cgcaacgcaaTCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAg                                                     <  1:3499754/148‑1 (MQ=255)
                                                                                                  caccgcaaTCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGtt                                                   <  2:2452234/145‑1 (MQ=255)
                                                                                                    acgcaaTCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTac                                                 <  1:3094645/148‑1 (MQ=255)
                                                                                                      gcaaTCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTaccg                                               <  1:2037423/148‑1 (MQ=255)
                                                                                                       caaTCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTaccga                                              <  2:3392908/148‑1 (MQ=255)
                                                                                                        aaTCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTaccgac                                             <  2:1731472/148‑1 (MQ=255)
                                                                                                         aTCAAACAGTTCGGCATTGTAGTCGTGGGTTTCGATGTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGTGTTGACGAAATCTTCGCTGAACTGTTGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTaccgacc                                            <  1:3314697/148‑1 (MQ=255)
                                                                                                          tCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTaccgaccg                                           >  2:268165/1‑148 (MQ=255)
                                                                                                            aaaCAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACcgcg                                         >  2:2150105/1‑148 (MQ=255)
                                                                                                             aaCAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCgcgc                                        >  1:3508407/1‑148 (MQ=255)
                                                                                                             aaCAGTTCGGCAAGGTCGGCGTGCGGGTCGATCTGGGTGGTGTACGGGTACCACGGAATCCCCCGCGAGGTGACGAACTCTTCGCGGAACGGATGCCAGTCAACTTCCGGGTCAGCGGCGATGTGGGCGTTATAGTTACCGACCgcgc                                        <  1:373137/148‑1 (MQ=255)
                                                                                                               cAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCg                                      >  1:1416451/1‑148 (MQ=255)
                                                                                                               cAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCg                                      >  1:2671832/1‑148 (MQ=255)
                                                                                                               cAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCg                                      <  2:2612891/148‑1 (MQ=255)
                                                                                                                aGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGt                                     <  2:696718/148‑1 (MQ=255)
                                                                                                                 gTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGtt                                    >  1:1472460/1‑148 (MQ=255)
                                                                                                                  ttCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTg                                   >  1:2202774/1‑148 (MQ=255)
                                                                                                                   tCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGa                                  >  2:2097620/1‑148 (MQ=255)
                                                                                                                    cGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGAt                                 <  1:1912465/148‑1 (MQ=255)
                                                                                                                      gCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGAttt                               <  2:2229433/148‑1 (MQ=255)
                                                                                                                      gCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGAttt                               >  2:2161023/1‑148 (MQ=255)
                                                                                                                       cAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGAtttt                              >  2:1103339/1‑148 (MQ=255)
                                                                                                                         aTGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGc                            >  1:2826291/1‑148 (MQ=255)
                                                                                                                         aTGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGc                            >  2:2317284/1‑148 (MQ=255)
                                                                                                                         aTGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGc                            >  1:2571524/1‑148 (MQ=255)
                                                                                                                            tAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGa                         >  2:1712336/1‑148 (MQ=255)
                                                                                                                              gTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAgg                       <  1:3350195/148‑1 (MQ=255)
                                                                                                                                   gCGGTTCGATCTGGGTTGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGAtctc                  <  1:681981/148‑1 (MQ=255)
                                                                                                                                   gCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGAtctc                  <  1:1878218/148‑1 (MQ=255)
                                                                                                                                   gCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGAtctc                  <  2:1904710/148‑1 (MQ=255)
                                                                                                                                      gTTCGATCTGGGTTGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCAc               <  2:3350879/148‑1 (MQ=255)
                                                                                                                                      gTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCAc               <  2:1259490/148‑1 (MQ=255)
                                                                                                                                      gTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCAc               <  2:2500066/148‑1 (MQ=255)
                                                                                                                                         cGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCACCTg            >  1:1078028/1‑148 (MQ=255)
                                                                                                                                          gATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCACCTgg           <  2:2892347/148‑1 (MQ=255)
                                                                                                                                           aTCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCACCTGGt          >  2:1942467/1‑148 (MQ=255)
                                                                                                                                           aTCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCACCTGGt          >  2:1557783/1‑148 (MQ=255)
                                                                                                                                            tCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCACCTGGtt         >  2:437165/1‑148 (MQ=255)
                                                                                                                                            tCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCACCTGGtt         >  1:2920187/1‑148 (MQ=255)
                                                                                                                                             cTGGGTTGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCACCTGGTTa        <  1:956099/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCAATCTCACCAGCAATGGTTTTCTGTTTGAAGTGGTTAAGGGCGATATAACCCCAGACGTCACGGTCAAAGTCGATCAGAATAGTGTTGAAGCGCGCAACGCAATCAAACAGTTCGGCAATGTAGTCGTGCGGTTCGATCTGGGTGGTGTACGGGTTCCACTGAATACCCAGCGAGGTGACGAACTCTTCGCTGAACTGATGCCAGTCAACTTCCGGGTAAGCGGCGATGTGGGCGTTATAGTTACCGACCGCGCCGTTGATTTTGCCGAGGATCTCCACCTGGTTAAGCTGG  >  CP009273/1186561‑1186855

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: ‑