Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,068,854 C→A 41.1% E230* (GAG→TAG)  yggF ← fructose 1,6 bisphosphatase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,068,8540CA41.1% 327.7 / inf 282E230* (GAG→TAG) yggFfructose 1,6 bisphosphatase isozyme
Reads supporting (aligned to +/- strand):  ref base C (86/80);  new base A (63/53);  total (149/133)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 285 total aligned reads displayed.

AAAGAGGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAA  >  CP009273/3068705‑3069003
                                                                                                                                                     |                                                                                                                                                     
aaaGAGGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTc                                                                                                                                                       <  2:1412044/150‑1 (MQ=255)
    agGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGcc                                                                                                                                                   >  1:1073414/1‑150 (MQ=255)
    agGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGcc                                                                                                                                                   >  1:2543584/1‑150 (MQ=255)
    agGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGcc                                                                                                                                                   <  1:806528/150‑1 (MQ=255)
    agGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGcc                                                                                                                                                   <  2:4353635/150‑1 (MQ=255)
    agGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGcc                                                                                                                                                   <  2:3909152/150‑1 (MQ=255)
         atCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCat                                                                                                                                              <  1:1254791/150‑1 (MQ=255)
         atCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCat                                                                                                                                              <  1:2938553/150‑1 (MQ=255)
            gTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATc                                                                                                                                           >  2:2178027/1‑150 (MQ=255)
            gTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATc                                                                                                                                           >  1:1988602/1‑150 (MQ=255)
            gTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATc                                                                                                                                           >  1:3578510/1‑150 (MQ=255)
               cccccTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATcccc                                                                                                                                        >  2:1051411/1‑150 (MQ=255)
                 accTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGc                                                                                                                                      >  2:4513811/2‑150 (MQ=255)
                     cACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTaa                                                                                                                                  >  2:557642/1‑150 (MQ=255)
                       ccAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATg                                                                                                                                <  1:888764/150‑1 (MQ=255)
                        cAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGc                                                                                                                               <  1:1023395/150‑1 (MQ=255)
                         aGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCt                                                                                                                              <  2:168114/150‑1 (MQ=255)
                          gTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCtt                                                                                                                             <  1:1487871/150‑1 (MQ=255)
                           ttaaTCGAGCGAGTACACACGGTTAACGTCAACGCCCAGTGGCTTGCTACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCGGCCTACTGCttt                                                                                                                            <  2:3132107/147‑1 (MQ=255)
                           ttCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCttt                                                                                                                            <  2:3483017/150‑1 (MQ=255)
                           ttCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCttt                                                                                                                            <  2:1114763/150‑1 (MQ=255)
                            tCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTa                                                                                                                           >  1:1958768/1‑150 (MQ=255)
                            tCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTa                                                                                                                           >  1:2621637/1‑150 (MQ=255)
                              aTCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTTCATATCCCCGCCTATTGCTTTAAc                                                                                                                         <  2:578399/150‑1 (MQ=255)
                              aTCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAAc                                                                                                                         >  1:3131572/1‑150 (MQ=255)
                              aTCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAAc                                                                                                                         <  1:4342870/150‑1 (MQ=255)
                              aTCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAAc                                                                                                                         >  2:2075246/1‑150 (MQ=255)
                                cgagcgagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACgg                                                                                                                       <  2:3045668/150‑1 (MQ=255)
                                cgagcgagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACgg                                                                                                                       <  2:1460464/150‑1 (MQ=255)
                                 gagcgagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGc                                                                                                                      >  2:3618779/1‑150 (MQ=255)
                                 gagcgagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGc                                                                                                                      >  1:179410/1‑150 (MQ=255)
                                      agTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAgg                                                                                                                 <  1:4096078/150‑1 (MQ=255)
                                       gTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGc                                                                                                                <  2:4406452/150‑1 (MQ=255)
                                        tACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCa                                                                                                               >  2:2912090/1‑150 (MQ=255)
                                         acacacGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAg                                                                                                              <  1:4491659/150‑1 (MQ=255)
                                          cacacGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGa                                                                                                             >  2:956754/1‑150 (MQ=255)
                                           acacGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaa                                                                                                            >  1:990697/1‑150 (MQ=255)
                                            cacGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaaa                                                                                                           <  2:3239713/150‑1 (MQ=255)
                                             acGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAAt                                                                                                          <  2:2131421/150‑1 (MQ=255)
                                             acGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAAt                                                                                                          >  2:4163479/1‑150 (MQ=255)
                                               ggTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCa                                                                                                        >  1:4516980/1‑150 (MQ=255)
                                                gTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCAc                                                                                                       <  2:1229549/150‑1 (MQ=255)
                                                gTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCAc                                                                                                       <  2:3547660/150‑1 (MQ=255)
                                                gTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCAc                                                                                                       <  2:2038655/150‑1 (MQ=255)
                                                gTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCCCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCAc                                                                                                       >  1:2391493/1‑150 (MQ=255)
                                                     cGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCtt                                                                                                  <  2:2568894/150‑1 (MQ=255)
                                                      gTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTc                                                                                                 <  2:3185251/150‑1 (MQ=255)
                                                           cGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGgcg                                                                                            <  1:2834241/150‑1 (MQ=255)
                                                            gCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGgcgc                                                                                           >  1:2479647/1‑150 (MQ=255)
                                                               cATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAcc                                                                                        >  1:3416499/1‑150 (MQ=255)
                                                                aTTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccg                                                                                       >  1:2328958/1‑150 (MQ=255)
                                                                    ccTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGa                                                                                   >  2:4216035/1‑150 (MQ=255)
                                                                    ccTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGa                                                                                   <  1:2806730/150‑1 (MQ=255)
                                                                    ccTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGa                                                                                   >  1:142711/1‑150 (MQ=255)
                                                                     cTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGAAATCCCCATAAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCAAAACCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGAt                                                                                  >  2:4129154/1‑150 (MQ=255)
                                                                      ttGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATg                                                                                 <  2:4012698/150‑1 (MQ=255)
                                                                       tGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATgg                                                                                <  2:168462/150‑1 (MQ=255)
                                                                        gCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgt                                                                               <  2:3679351/150‑1 (MQ=255)
                                                                            cGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCTCCTAATGATTTAACGGCACAGGCAGAAAACACGCCTTCTGGAGCACCGCCGAAGGTGTac                                                                           >  1:3301611/1‑150 (MQ=255)
                                                                              cTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacat                                                                         >  2:601733/1‑150 (MQ=255)
                                                                                  cgcTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcaca                                                                     <  1:3208699/150‑1 (MQ=255)
                                                                                  cgcTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCGGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcaca                                                                     <  1:3208781/150‑1 (MQ=255)
                                                                                    cTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcacatc                                                                   >  1:1693463/1‑150 (MQ=255)
                                                                                     tCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCg                                                                  >  2:3468781/1‑150 (MQ=255)
                                                                                      ccTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGt                                                                 >  2:3146611/1‑150 (MQ=255)
                                                                                       cTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTCCATCACATCGTa                                                                >  1:3884492/1‑150 (MQ=255)
                                                                                       cTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTa                                                                <  1:4571095/150‑1 (MQ=255)
                                                                                           tctGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGa                                                            >  1:1278407/1‑150 (MQ=255)
                                                                                            ctGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGAt                                                           >  2:4380227/1‑150 (MQ=255)
                                                                                                aaTCTGCCGATTTTCCNTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCNGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATc                                                       >  2:2360655/1‑150 (MQ=255)
                                                                                                 aTCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATcc                                                      >  1:1364760/1‑150 (MQ=255)
                                                                                                   cTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCGCATCGTAAGGATTATCCTg                                                    <  1:3699457/150‑1 (MQ=255)
                                                                                                     gCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGcc                                                  <  1:1314744/150‑1 (MQ=255)
                                                                                                     gCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGcc                                                  <  2:1584953/150‑1 (MQ=255)
                                                                                                      ccGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCca                                                 <  2:1756071/150‑1 (MQ=255)
                                                                                                        gATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCcagc                                               >  2:2720397/1‑150 (MQ=255)
                                                                                                        gATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCcagc                                               >  2:85453/1‑150 (MQ=255)
                                                                                                         aTTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCcagca                                              >  1:2911800/1‑150 (MQ=255)
                                                                                                             tCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTc                                          <  1:1583911/150‑1 (MQ=255)
                                                                                                              ccGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCa                                         >  2:1373632/1‑150 (MQ=255)
                                                                                                                 tgtAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCa                                      >  2:820089/1‑150 (MQ=255)
                                                                                                                 tgtAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCa                                      <  1:483473/150‑1 (MQ=255)
                                                                                                                  gtAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCATTCTGGCGCACCGCCGATGGTGTACATCACATCGGAAGGATTATCCTGCCAGCAGGTCAGCAc                                     >  1:4471561/1‑150 (MQ=255)
                                                                                                                    aaTCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGAGTATCCTGCCAGCAGGTCAGCACGc                                   >  2:225173/1‑150 (MQ=255)
                                                                                                                      tCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTg                                 <  1:673338/150‑1 (MQ=255)
                                                                                                                       ccccTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTgg                                >  1:3841039/1‑150 (MQ=255)
                                                                                                                          cTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTggcgg                             >  2:3385378/1‑150 (MQ=255)
                                                                                                                             tAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGa                          >  1:1603632/1‑150 (MQ=255)
                                                                                                                              aGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGAc                         >  1:4013945/1‑150 (MQ=255)
                                                                                                                                   ggCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGc                    >  2:720472/1‑150 (MQ=255)
                                                                                                                                       aaaGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATc                >  1:781350/1‑150 (MQ=255)
                                                                                                                                         aGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCg              >  1:3614529/1‑150 (MQ=255)
                                                                                                                                         aGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACACCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCg              >  1:1471549/1‑150 (MQ=255)
                                                                                                                                            cAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCa           >  2:3770119/1‑150 (MQ=255)
                                                                                                                                             aaTCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAg          >  1:4162049/1‑150 (MQ=255)
                                                                                                                                                 aGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCa      <  2:3869758/150‑1 (MQ=255)
                                                                                                                                                   cTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCaaa    >  2:1888201/1‑150 (MQ=255)
                                                                                                                                                   cTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCaaa    <  2:1958768/150‑1 (MQ=255)
                                                                                                                                                    tCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCaaaa   >  1:2987653/1‑150 (MQ=255)
                                                                                                                                                    tCGGCATGCATATCACCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCaaaa   >  2:1960624/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAGAGGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAA  >  CP009273/3068705‑3069003

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: