Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,207,552 G→T 53.0% M172I (ATG→ATT aceA → isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,207,5520GT53.0% 29.9 / inf 232M172I (ATG→ATTaceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (46/63);  new base T (68/55);  total (114/118)
Fisher's exact test for biased strand distribution p-value = 4.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

GATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCTGAC  >  CP009273/4207405‑4207699
                                                                                                                                                   |                                                                                                                                                   
gATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGaa                                                                                                                                                   >  2:2426370/1‑150 (MQ=255)
  tCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAg                                                                                                                                                 <  1:3477277/150‑1 (MQ=255)
  tCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAg                                                                                                                                                 >  1:1636349/1‑150 (MQ=255)
  tCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAg                                                                                                                                                 >  1:4072189/1‑150 (MQ=255)
   caacaaCACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGc                                                                                                                                                <  1:3575530/150‑1 (MQ=255)
    aacaacACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGct                                                                                                                                               <  2:2247723/150‑2 (MQ=255)
    aacaacACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGct                                                                                                                                               >  1:1739403/1‑149 (MQ=255)
     acaacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCga                                                                                                                                              >  1:1230387/1‑150 (MQ=255)
      caacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCtat                                                                                                                                             <  1:262016/150‑1 (MQ=255)
      caacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCtat                                                                                                                                             <  1:3203295/150‑1 (MQ=255)
        acacCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgatga                                                                                                                                           >  2:697431/1‑150 (MQ=255)
          acCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGAtt                                                                                                                                         >  1:3129056/1‑150 (MQ=255)
          acCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGAtt                                                                                                                                         >  1:3129033/1‑150 (MQ=255)
            cTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGa                                                                                                                                       <  2:1172892/150‑1 (MQ=255)
              tCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAg                                                                                                                                     >  1:1696185/1‑150 (MQ=255)
              tCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAg                                                                                                                                     <  2:652226/150‑1 (MQ=255)
                    gtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTg                                                                                                                               <  2:1369302/150‑1 (MQ=255)
                    gtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTg                                                                                                                               <  2:3462621/150‑1 (MQ=255)
                      gCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCa                                                                                                                             <  2:106541/150‑1 (MQ=255)
                      gCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCa                                                                                                                             <  1:2781858/150‑1 (MQ=255)
                      gCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCa                                                                                                                             >  1:257600/1‑150 (MQ=255)
                       ccGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAg                                                                                                                            <  1:1179268/150‑1 (MQ=255)
                       ccGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGAAg                                                                                                                            >  2:1550182/1‑150 (MQ=255)
                       ccGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGGTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAg                                                                                                                            <  1:288280/150‑1 (MQ=255)
                         gATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCg                                                                                                                          <  2:1022716/150‑1 (MQ=255)
                         gATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCg                                                                                                                          >  2:2085245/1‑150 (MQ=255)
                          aTCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCgg                                                                                                                         <  2:937259/150‑1 (MQ=255)
                          aTCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCgg                                                                                                                         >  1:3265008/1‑150 (MQ=255)
                              gATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGt                                                                                                                     <  2:3949195/150‑1 (MQ=255)
                                    aTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACtt                                                                                                               <  1:1537079/150‑1 (MQ=255)
                                    aTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACtt                                                                                                               >  1:68463/1‑150 (MQ=255)
                                       gTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGa                                                                                                            <  2:2231707/150‑1 (MQ=255)
                                         ccGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAg                                                                                                          <  1:1844245/150‑1 (MQ=255)
                                          cgcgGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGa                                                                                                         <  2:3032598/150‑1 (MQ=255)
                                           gcgGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCTGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGTAGCCGGTGCAGCGGCAGTTCACTTCGAAGAt                                                                                                        <  1:561093/150‑1 (MQ=255)
                                           gcgGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGAt                                                                                                        <  2:1835250/150‑1 (MQ=255)
                                            cgGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGACTGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATc                                                                                                       <  1:1781935/150‑1 (MQ=255)
                                             gggCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCa                                                                                                      >  2:228000/1‑150 (MQ=255)
                                             gggCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCa                                                                                                      >  2:1006295/1‑150 (MQ=255)
                                             gggCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCa                                                                                                      <  2:1605205/150‑1 (MQ=255)
                                                cATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCt                                                                                                   <  1:1496745/150‑1 (MQ=255)
                                                  ttGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTgg                                                                                                 <  1:946497/150‑1 (MQ=255)
                                                   tGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGc                                                                                                >  1:1203006/1‑150 (MQ=255)
                                                   tGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGc                                                                                                <  2:164862/150‑1 (MQ=255)
                                                     aGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGt                                                                                              >  1:4056293/1‑150 (MQ=255)
                                                        cGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAg                                                                                           >  1:249913/1‑150 (MQ=255)
                                                         gggCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGt                                                                                          >  2:2928451/1‑150 (MQ=255)
                                                          ggCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTg                                                                                         <  2:667509/150‑1 (MQ=255)
                                                            cGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaa                                                                                       >  1:3353778/1‑150 (MQ=255)
                                                            cGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaa                                                                                       >  1:2795799/1‑150 (MQ=255)
                                                             gATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaag                                                                                      >  1:3350590/1‑150 (MQ=255)
                                                             gATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaag                                                                                      >  2:2013991/1‑150 (MQ=255)
                                                             caacagcgcTATATCGAGAACTTCCTGACGATCTTTGCCGATGCGGAAGCCGGTATTGGCGATGTCCTGAATGCCATTGAACTGAAGAAAGAAATGATAGAAGCCCGTTCAGCGGCAGTTCACTTCGAAGATCACCTGGCTTCAGTgaag                                                                                      >  2:107009/6‑150 (MQ=255)
                                                              aTCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaaga                                                                                     >  2:575123/1‑150 (MQ=255)
                                                              aTCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaaga                                                                                     <  2:2600796/150‑1 (MQ=255)
                                                                   cgcTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGc                                                                                >  2:1483927/1‑150 (MQ=255)
                                                                   cgcTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGc                                                                                >  2:785442/1‑150 (MQ=255)
                                                                     cTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCgg                                                                              >  1:774312/1‑150 (MQ=255)
                                                                        tGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTca                                                                           <  1:1673723/150‑1 (MQ=255)
                                                                        tGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTca                                                                           >  1:1038539/1‑150 (MQ=255)
                                                                             aTTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATgg                                                                      <  2:1056837/150‑1 (MQ=255)
                                                                               tACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGggc                                                                    >  1:145089/1‑150 (MQ=255)
                                                                                aCTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGggcg                                                                   >  2:2744512/1‑150 (MQ=255)
                                                                                aCTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGggcg                                                                   <  1:2243304/150‑1 (MQ=255)
                                                                                   tCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCa                                                                >  2:3887908/1‑150 (MQ=255)
                                                                                   tCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCa                                                                <  2:598351/150‑1 (MQ=255)
                                                                                    ccTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCaa                                                               >  1:2843963/1‑150 (MQ=255)
                                                                                     cTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCaaa                                                              <  2:4168547/150‑1 (MQ=255)
                                                                                      tGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAg                                                             <  2:2465875/150‑1 (MQ=255)
                                                                                      tGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAg                                                             <  2:2069529/150‑1 (MQ=255)
                                                                                          gATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                         >  2:2555452/1‑150 (MQ=255)
                                                                                          gATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                         >  1:3894335/1‑150 (MQ=255)
                                                                                           aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTa                                                        <  2:3550540/150‑1 (MQ=255)
                                                                                           aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTa                                                        <  2:3775445/150‑1 (MQ=255)
                                                                                              gTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTg                                                     >  2:3777235/1‑150 (MQ=255)
                                                                                                 gCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCa                                                  >  1:1252768/1‑150 (MQ=255)
                                                                                                   cGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAAc                                                <  2:110175/150‑1 (MQ=255)
                                                                                                   cGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAAc                                                >  2:2901842/1‑150 (MQ=255)
                                                                                                      tGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCa                                             >  2:2633938/1‑150 (MQ=255)
                                                                                                        cGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAgg                                           >  2:729528/1‑150 (MQ=255)
                                                                                                         ggAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGa                                          >  2:1679822/1‑150 (MQ=255)
                                                                                                            aGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGc                                       >  2:930827/1‑150 (MQ=255)
                                                                                                             gCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCt                                      <  1:827606/150‑1 (MQ=255)
                                                                                                              ccGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTa                                     <  2:4137185/150‑1 (MQ=255)
                                                                                                               cGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTAt                                    >  2:252420/1‑150 (MQ=255)
                                                                                                                ggTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTAtt                                   <  2:214060/150‑1 (MQ=255)
                                                                                                                ggTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTAtt                                   <  2:1634095/150‑1 (MQ=255)
                                                                                                                ggTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTAtt                                   <  2:1634052/150‑1 (MQ=255)
                                                                                                                 gTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTc                                  >  2:2789720/1‑150 (MQ=255)
                                                                                                                  ttttGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCa                                 <  2:3110372/150‑1 (MQ=255)
                                                                                                                    ttGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGa                               >  2:3407164/1‑150 (MQ=255)
                                                                                                                     tGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTTGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGaa                              <  2:2114824/150‑1 (MQ=255)
                                                                                                                        cGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACt                           >  1:2325720/1‑150 (MQ=255)
                                                                                                                          gtgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTgg                         >  1:2359304/1‑150 (MQ=255)
                                                                                                                          gtgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTgg                         <  1:2651736/150‑1 (MQ=255)
                                                                                                                                 gAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCgcggc                  >  1:2730152/1‑150 (MQ=255)
                                                                                                                                           gAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTTAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCa        <  1:648593/150‑1 (MQ=255)
                                                                                                                                            aaCTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGAGAAGAAATGCGGTCACATGTGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAg       >  1:3865793/1‑150 (MQ=255)
                                                                                                                                               tgattaAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtg    <  2:3993399/150‑1 (MQ=255)
                                                                                                                                                 atgaAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtgac  <  2:3228039/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCTGAC  >  CP009273/4207405‑4207699

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: