Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,931,021 C→A 88.6% A123E (GCG→GAG)  rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,931,0210CA88.6% 475.4 / 59.5 193A123E (GCG→GAG) rbsKribokinase
Reads supporting (aligned to +/- strand):  ref base C (14/8);  new base A (98/73);  total (112/81)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 198 total aligned reads displayed.

CCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCG  >  CP009273/3930875‑3931160
                                                                                                                                                  |                                                                                                                                           
ccAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGgcg                                                                                                                                            >  1:395661/1‑148 (MQ=255)
 cagcagCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGgagc                                                                                                                                           <  2:121782/148‑1 (MQ=255)
 cagcagCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGgagc                                                                                                                                           <  2:2348394/148‑1 (MQ=255)
    cagcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtgg                                                                                                                                        <  2:1246006/148‑1 (MQ=255)
    cagcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtgg                                                                                                                                        >  2:2974923/1‑148 (MQ=255)
    cagcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtgg                                                                                                                                        <  2:1937635/148‑1 (MQ=255)
        tCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGa                                                                                                                                    <  1:720774/148‑1 (MQ=255)
        tCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGa                                                                                                                                    <  1:1515857/148‑1 (MQ=255)
         cGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGaa                                                                                                                                   >  2:559224/1‑148 (MQ=255)
             aCTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAgcgc                                                                                                                               >  2:3201852/1‑148 (MQ=255)
              cTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAATCGCa                                                                                                                              >  1:1734602/1‑148 (MQ=255)
              cTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCa                                                                                                                              >  1:3101557/1‑148 (MQ=255)
              cTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCa                                                                                                                              >  1:1332205/1‑148 (MQ=255)
               tGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCaa                                                                                                                             <  1:2418759/148‑1 (MQ=255)
                gATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAAc                                                                                                                            >  1:2459161/1‑148 (MQ=255)
                 aTAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACg                                                                                                                           <  2:2078858/148‑1 (MQ=255)
                 aTAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACg                                                                                                                           <  2:1836864/148‑1 (MQ=255)
                   aaCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTg                                                                                                                         <  1:861776/148‑1 (MQ=255)
                    aCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGa                                                                                                                        >  2:2946586/1‑148 (MQ=255)
                     cATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAg                                                                                                                       >  2:344270/1‑148 (MQ=255)
                      aTTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGc                                                                                                                      >  2:651133/1‑148 (MQ=255)
                        tGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                    >  2:2173330/1‑148 (MQ=255)
                        tGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                    >  1:1104190/1‑148 (MQ=255)
                        tGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                    <  1:1221283/148‑1 (MQ=255)
                        tGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                    >  1:2755594/1‑148 (MQ=255)
                        tGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                    >  2:2801552/1‑148 (MQ=255)
                        tGATAGTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTGTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                    <  1:3303982/148‑1 (MQ=255)
                             ttACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGc                                                                                                               >  1:2198283/1‑148 (MQ=255)
                             ttACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGc                                                                                                               >  1:887488/1‑148 (MQ=255)
                             ttACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGc                                                                                                               >  1:2130759/1‑148 (MQ=255)
                              tACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGcc                                                                                                              >  1:2926779/1‑148 (MQ=255)
                              tACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGcc                                                                                                              >  2:2660989/1‑148 (MQ=255)
                                 tCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAAc                                                                                                           >  2:1736104/1‑148 (MQ=255)
                                  ccGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAAcg                                                                                                          >  1:317622/1‑148 (MQ=255)
                                        aGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAg                                                                                                    <  2:2110096/148‑1 (MQ=255)
                                         gCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGc                                                                                                   <  1:1759819/148‑1 (MQ=255)
                                            tGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCatt                                                                                                >  2:539674/1‑148 (MQ=255)
                                            tGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAAGGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCatt                                                                                                <  2:786629/148‑1 (MQ=255)
                                             gATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCatta                                                                                               >  2:1822652/1‑148 (MQ=255)
                                                cAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCAttatta                                                                                            >  2:3257707/1‑148 (MQ=255)
                                                 aaaGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTaa                                                                                           <  1:378796/148‑1 (MQ=255)
                                                 aaaGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTaa                                                                                           <  1:1133459/148‑1 (MQ=255)
                                                 aaaGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTaa                                                                                           >  2:3071316/1‑148 (MQ=255)
                                                  aaGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAAt                                                                                          <  2:1304405/148‑1 (MQ=255)
                                                  aaGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAAt                                                                                          <  1:1267223/148‑1 (MQ=255)
                                                     gCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCa                                                                                       <  1:1897573/148‑1 (MQ=255)
                                                      cGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAg                                                                                      <  1:1409308/148‑1 (MQ=255)
                                                         aTCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTg                                                                                   >  2:1956512/1‑148 (MQ=255)
                                                          tCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTgg                                                                                  >  2:3125677/1‑148 (MQ=255)
                                                           cAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCCGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGa                                                                                 >  1:3207938/1‑148 (MQ=255)
                                                            aaCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGaa                                                                                <  2:3105187/148‑1 (MQ=255)
                                                            aaCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGaa                                                                                <  1:2880480/148‑1 (MQ=255)
                                                             aCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAAt                                                                               <  2:1052253/148‑1 (MQ=255)
                                                              cAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                                                              >  2:1983384/1‑148 (MQ=255)
                                                                 gtgtGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcacc                                                                           >  2:68679/1‑148 (MQ=255)
                                                                 gtgtGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcacc                                                                           >  2:1286751/1‑148 (MQ=255)
                                                                    tgGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCATCTGGAATCACCACt                                                                        >  1:2954485/1‑148 (MQ=255)
                                                                        gcTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGaa                                                                    <  1:1878209/148‑1 (MQ=255)
                                                                         cTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGaaa                                                                   >  1:670899/1‑148 (MQ=255)
                                                                          tGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAg                                                                  >  1:2250208/1‑148 (MQ=255)
                                                                           gATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgt                                                                 >  2:20573/1‑148 (MQ=255)
                                                                             tttttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgtgt                                                               <  2:2308426/148‑1 (MQ=255)
                                                                             tttttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgtgt                                                               >  2:361616/1‑148 (MQ=255)
                                                                              ttttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGNTGGAAGCGCAACGAGAGCGTATTGCCAACGCGTCAGCATTAATAATGCAGCTGGAATCACCACTCGAAAGtgtg                                                              >  1:1680212/1‑148 (MQ=255)
                                                                              ttttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGtgtg                                                              >  2:465855/1‑148 (MQ=255)
                                                                                     aaTGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAg                                                       <  1:2083110/148‑1 (MQ=255)
                                                                                       tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcg                                                     >  2:2198769/1‑148 (MQ=255)
                                                                                       tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcg                                                     >  1:984997/1‑148 (MQ=255)
                                                                                        ggCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcgg                                                    >  2:3304924/1‑148 (MQ=255)
                                                                                               gTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaatt                                             <  2:422246/148‑3 (MQ=255)
                                                                                               gTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAAt                                             <  1:2660989/148‑1 (MQ=255)
                                                                                                tGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATc                                            >  2:1007012/1‑148 (MQ=255)
                                                                                                tGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATc                                            >  2:704576/1‑148 (MQ=255)
                                                                                                tGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATc                                            >  2:2578120/1‑148 (MQ=255)
                                                                                                  agaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGc                                          <  2:112480/148‑1 (MQ=255)
                                                                                                      tGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCAt                                      >  1:915210/1‑148 (MQ=255)
                                                                                                          aTCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCaaa                                  <  1:2801588/148‑1 (MQ=255)
                                                                                                           tCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCaaaa                                 >  1:1363515/1‑148 (MQ=255)
                                                                                                           tCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCaaaa                                 <  1:559224/148‑1 (MQ=255)
                                                                                                           tCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCaaaa                                 >  1:2137128/1‑148 (MQ=255)
                                                                                                              gTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATaa                              >  1:514272/1‑148 (MQ=255)
                                                                                                               tATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAg                             >  2:2261291/1‑148 (MQ=255)
                                                                                                               tATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAg                             >  1:2531454/1‑148 (MQ=255)
                                                                                                                aTTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGa                            <  1:948116/148‑1 (MQ=255)
                                                                                                                 ttCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGAc                           >  2:152403/1‑148 (MQ=255)
                                                                                                                  tCATGCCGGCGCTAATGCTGCCCTGTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACt                          <  2:2909479/148‑1 (MQ=255)
                                                                                                                   cATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTa                         <  2:284915/148‑1 (MQ=255)
                                                                                                                    aTGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTAt                        <  2:887488/148‑1 (MQ=255)
                                                                                                                      gCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCg                      >  1:756950/1‑148 (MQ=255)
                                                                                                                      gCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCg                      >  1:2393950/1‑148 (MQ=255)
                                                                                                                           cgcTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTgcg                 >  1:1034217/1‑148 (MQ=255)
                                                                                                                            gcTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTgcgc                <  1:2812040/148‑1 (MQ=255)
                                                                                                                              tAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCtt              >  1:2766075/1‑148 (MQ=255)
                                                                                                                                aTGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTaa            >  1:1007747/1‑148 (MQ=255)
                                                                                                                                 tgctgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAc           >  2:1377384/1‑148 (MQ=255)
                                                                                                                                 tgctgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAc           >  2:839619/1‑148 (MQ=255)
                                                                                                                                  gctgcCCTGTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAcc          <  2:1687483/148‑1 (MQ=255)
                                                                                                                                      cccTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggc      >  1:2062768/1‑148 (MQ=255)
                                                                                                                                      cccTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggc      >  1:1596624/1‑148 (MQ=255)
                                                                                                                                          ttCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggctccg  >  1:1438666/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
CCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCG  >  CP009273/3930875‑3931160

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