Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,085,145 A→C 5.9% A194A (GCT→GCG)  frvB ← putative PTS enzyme, IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,085,1450AC 5.9% 926.9 / 4.5 370A194A (GCT→GCG) frvBputative PTS enzyme, IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base A (177/158);  new base C (3/31);  total (180/190)
Fisher's exact test for biased strand distribution p-value = 3.16e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 375 total aligned reads displayed.

ACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGCCACCAGGATACCACCGCC  >  CP009273/4084996‑4085294
                                                                                                                                                     |                                                                                                                                                     
aCGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATaa                                                                                                                                                       >  1:4205133/1‑150 (MQ=255)
   aaaTAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAgcg                                                                                                                                                    <  2:4093325/150‑1 (MQ=255)
   aaaTAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAgcg                                                                                                                                                    >  2:4425125/1‑150 (MQ=255)
   aaaTAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAgcg                                                                                                                                                    >  1:4086290/1‑150 (MQ=255)
    aaTAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAgcgc                                                                                                                                                   >  1:4680954/1‑150 (MQ=255)
     aTAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGcc                                                                                                                                                  <  2:3627541/150‑1 (MQ=255)
      tAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGccc                                                                                                                                                 >  1:4317772/1‑150 (MQ=255)
       aGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCa                                                                                                                                                >  2:1637461/1‑150 (MQ=255)
            aTCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCatgatg                                                                                                                                           >  2:3060795/1‑150 (MQ=255)
                ccAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGACAATCGACGAGGCAATATAAGCGCCCATGATGGGGa                                                                                                                                       >  1:493534/1‑150 (MQ=255)
                 cAGCCCCAGCACCACCGCGCCGTGTAATCCCGCACCCGACTGGGTGCCGAGATGCGATTTGTGGTTCGCCATATAGCAAACCATAAACGCGGGAGCGAAAGCCGTTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGAt                                                                                                                                      <  1:3947307/150‑1 (MQ=255)
                 cAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGAt                                                                                                                                      >  2:4114424/1‑150 (MQ=255)
                    ccccAGCACCCCCGCGCCGCGTAATCCCGCACCCGACTGGGTGCCGCGCCGCGTTTTGGGGTTCGCAAGATAGAAAACCAGAAACGCCGGAGCTAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCAt                                                                                                                                   <  1:266049/150‑1 (MQ=255)
                    ccccAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCAt                                                                                                                                   <  2:1795007/150‑1 (MQ=255)
                      ccAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATaa                                                                                                                                 >  2:1687605/1‑150 (MQ=255)
                      ccAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATaa                                                                                                                                 >  2:256417/1‑150 (MQ=255)
                       cAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATaaa                                                                                                                                >  1:3573479/1‑150 (MQ=255)
                       cAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATaaa                                                                                                                                >  1:726621/1‑150 (MQ=255)
                        aGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAg                                                                                                                               >  1:3744126/1‑150 (MQ=255)
                         gCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAgg                                                                                                                              >  1:1678150/1‑150 (MQ=255)
                          caccacCGCGCCGAGTAATCCCGCACCCGACTGGGTGCCGAGCAGCGATTTGTCGTTCGCCAGATAGAAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGt                                                                                                                             <  2:1054376/150‑1 (MQ=255)
                          caccacCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGt                                                                                                                             >  2:708372/1‑150 (MQ=255)
                           accaccGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTg                                                                                                                            >  1:1146972/1‑150 (MQ=255)
                             caccGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCACTATAAGCGCCCATGATGGGGATCATAAAGGTGaa                                                                                                                          >  1:3594726/1‑150 (MQ=255)
                              accGCGCCGTGGAATCCCGCCCCGGACTGGGTGCCGCGCAGCGCTTTGTCGTTCGCCTTATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAg                                                                                                                         <  1:4722383/150‑1 (MQ=255)
                                cgcgCCGAGGAATCCCGCACCCGACGGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGcc                                                                                                                       <  2:290148/150‑1 (MQ=255)
                                 gcgcCGTGTAACCCCGCACCCGCGGGGGTGCCGAGCTGGGTTTTTTCTTTCGCCTTTTCCAAAACCACAAACGCCGGAGCTAAAGCCGTTTTATCGGAAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCa                                                                                                                      <  1:1619697/150‑1 (MQ=255)
                                   gcCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAg                                                                                                                    <  2:4509460/150‑1 (MQ=255)
                                        ggATTCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAAc                                                                                                               <  2:1958670/150‑1 (MQ=255)
                                        ggAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAAc                                                                                                               >  2:4226325/1‑150 (MQ=255)
                                         gAATCCCGCCCCCGAGTGGGTGCCGAGCAGCGGTTTTTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCTAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAAcc                                                                                                              <  2:804988/150‑1 (MQ=255)
                                         gAACCCCGCACCGGTGGGGGTGCGTAGCCGCTATTTTTGGTTCGCAATATACAAAACCATAACCGCGGGAGCTAACGCGGTTTTATGGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAAcc                                                                                                              <  1:1273216/150‑1 (MQ=255)
                                          aaTCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCg                                                                                                             >  1:4527399/1‑150 (MQ=255)
                                          aaTCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCg                                                                                                             <  2:1426531/150‑1 (MQ=255)
                                           aTCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGa                                                                                                            >  1:2824907/1‑150 (MQ=255)
                                           aTCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGa                                                                                                            >  1:2271016/1‑150 (MQ=255)
                                            gcccGCCCCCGACGGGGTGCCGAGCCGCGCTTTGTCGTTCGCCATATAGCAAACCAGAACCGCCGGAGCGAAAGCCGGTTTAGCGGAAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGAt                                                                                                           <  1:4552044/149‑1 (MQ=255)
                                              ccGCACCCGACTGGGTGCCGCGCCGCGTTTTGTCGTTCGCAATATAGCAAACCACAACCGCCGGAGCTAAAGCCGTTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATgg                                                                                                         <  1:583644/150‑1 (MQ=255)
                                               cGCCCCCGACTGTGTGCCGAGCAGCGCTTTGTCGTGCGCCAGATAGCAAACCTGACACGCCGGAGCGACAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGa                                                                                                        <  2:612288/150‑1 (MQ=255)
                                               cGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGa                                                                                                        >  1:3677446/1‑150 (MQ=255)
                                                 cACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCATAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGAtt                                                                                                      <  2:1234181/150‑1 (MQ=255)
                                                   cccGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCCAGGTAACCGCTCGATTCa                                                                                                    >  2:4705635/1‑150 (MQ=255)
                                                     cGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCACTATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAAc                                                                                                  >  1:2453302/1‑150 (MQ=255)
                                                     cGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAAc                                                                                                  >  1:834823/1‑150 (MQ=255)
                                                      gACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGAGTCAACt                                                                                                 >  2:2724254/1‑150 (MQ=255)
                                                       aCTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTa                                                                                                >  1:1833377/1‑150 (MQ=255)
                                                         tGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTAcc                                                                                              >  1:1997582/1‑150 (MQ=255)
                                                           ggTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCa                                                                                            >  1:951539/1‑150 (MQ=255)
                                                           ggTGCCGAGCAGCGCTTTGTCGTTCGCAAGATAGCAAACCAGACACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCa                                                                                            <  1:1110926/150‑1 (MQ=255)
                                                           ggTGCCGAGCAGCGCTTGGTCGTTCGCCAGAGAGCAAACCAGCACCGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCACTATAAGCGCCCAGGAGGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCa                                                                                            <  2:4332402/150‑1 (MQ=255)
                                                            gTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGACACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAg                                                                                           <  1:2703630/150‑1 (MQ=255)
                                                             tGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAgg                                                                                          <  2:3556407/150‑1 (MQ=255)
                                                               ccGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAACCGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTa                                                                                        <  2:2259319/150‑1 (MQ=255)
                                                                cGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTaa                                                                                       <  2:3821001/150‑1 (MQ=255)
                                                                cGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTaa                                                                                       <  2:3337552/150‑1 (MQ=255)
                                                                cGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTaa                                                                                       >  2:2114367/1‑150 (MQ=255)
                                                                 gAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTaaa                                                                                      >  2:3456507/1‑150 (MQ=255)
                                                                  agcagcGTTTTGTCGTTCGCCATATACAAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAac                                                                                     <  1:4583660/150‑1 (MQ=255)
                                                                   gcagcGCTTTGTCGTTCGCCAGATAGAAACCCAGACACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATACGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAacg                                                                                    <  1:4619960/150‑1 (MQ=255)
                                                                                 ttCGCCAGATAGCAAACCAGAACCGCGGGAGCGAAAGCCGTTTTATGGGCAATCGACGAGGCAATAAAAGCGCCCATGAGGGGGATCATAAAGGTGAAGCCAAGGTAACCGTTGGATTCAACTACCCCGGTAAACGACGGCGCGCCTTTc                                                                      <  1:796390/150‑1 (MQ=255)
                                                                                   cGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGTTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGa                                                                    <  2:4127458/150‑1 (MQ=255)
                                                                                      cAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACat                                                                 >  2:370599/1‑150 (MQ=255)
                                                                                      atgATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGAAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACat                                                                 <  2:4590488/148‑1 (MQ=255)
                                                                                       ttttAGAAAACCAGAACCGCGGTTGGTAACGCGGTTTTATCGGAAATCGACGAGGAAATATAAGCGCCATTGTTGGGGATCATAAAGGTGAAGCCATGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACata                                                                <  1:1017918/147‑1 (MQ=255)
                                                                                       aGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACata                                                                >  1:3522457/1‑150 (MQ=255)
                                                                                         aTAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATc                                                              >  2:1693622/1‑150 (MQ=255)
                                                                                          tAGCAAACCAGACACGCCGGCGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATACGCGCCCATGATGGGGATCATAAAGGTGACGCCAAGGTAACCGATGGATTCAACTACCCCGGTAAACGACGGCGCGCCTTTCGACATATCg                                                             <  2:2604300/150‑1 (MQ=255)
                                                                                           aGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCgg                                                            >  2:3674729/1‑150 (MQ=255)
                                                                                              aaaCCACAACCGCGGGAGAGAATGCGGGTTTAGCGGCAATCGCCGAGTCAATATACGCGCCCATGATGGGGATCATAACGGTGAAGCCAAGGTAACCGATGTATTCAACTACCCCGGTAAACGACGGCGCGCCTTTCGACATATCGGtat                                                         <  1:1347119/150‑1 (MQ=255)
                                                                                               aaCCATAACCGCGGGTGCGAAAGCCGTTTTTTCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGtata                                                        <  1:4706986/150‑1 (MQ=255)
                                                                                                aCCAGAAACGCCGTAGCGAAAGCGGGTTTATCGGAAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATAc                                                       <  2:1142622/150‑1 (MQ=255)
                                                                                                 cccgAACCGCCGGAGCGAAAGCCGTTTTATCGGCAATCGACGAGGAAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACg                                                      <  2:3241040/147‑1 (MQ=255)
                                                                                                     aaaCGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCaa                                                  >  1:3353477/1‑150 (MQ=255)
                                                                                                      aaCGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCaaa                                                 >  2:1357397/1‑150 (MQ=255)
                                                                                                       aCGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAAc                                                >  1:1324475/1‑150 (MQ=255)
                                                                                                        cGCCGGCGCGAAAGCCGGTTTATCGGCAATCGCCGCGGAAATATCAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCCGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAAcc                                               <  1:1529952/150‑1 (MQ=255)
                                                                                                              aGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTg                                         >  1:3409372/1‑150 (MQ=255)
                                                                                                                cGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCACACCGAACTGCa                                       >  2:2367593/1‑150 (MQ=255)
                                                                                                                cGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCa                                       >  1:662294/1‑150 (MQ=255)
                                                                                                                cGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCa                                       >  1:2936426/1‑150 (MQ=255)
                                                                                                                  aaaGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCAcc                                     >  1:14820/1‑150 (MQ=255)
                                                                                                                     gCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGc                                  >  1:1794775/1‑150 (MQ=255)
                                                                                                                       ggttttATGGTCAATCGACGAGGAAATAAACGCGCCCATGATGGGGATCATAAAGGTGAAGCCACGGTAACCTATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCAt                                <  1:2124671/147‑1 (MQ=255)
                                                                                                                       cGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCAt                                >  2:568994/1‑150 (MQ=255)
                                                                                                                       cGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACAGAACTGCACCAGCAt                                >  1:4007611/1‑150 (MQ=255)
                                                                                                                        ggTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATg                               >  2:1751739/1‑150 (MQ=255)
                                                                                                                         gTTTATCGGCAATCGACGAGGCAATATACGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGGAAACGACGGCGCGCCTTTCGACATAGCGGTATACGGCAAACCGAACTGCACCAGCATGt                              >  2:3531267/1‑150 (MQ=255)
                                                                                                                           ttATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTg                            >  1:617213/1‑150 (MQ=255)
                                                                                                                            tATGGGCAATCGACGTGGAAATATAAGCGCCCATTATGGGGATAATAAAGGTGAAGCCTAGGTAACCGTTGTATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTgg                           <  2:3991116/150‑1 (MQ=255)
                                                                                                                             aTCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGc                          >  2:2044949/1‑150 (MQ=255)
                                                                                                                              tCGGCAATCGACGAGGCAATATACGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCa                         <  2:1537574/150‑1 (MQ=255)
                                                                                                                              tCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCa                         >  2:4635971/1‑150 (MQ=255)
                                                                                                                              ggggCAATCGACGCGGCAATATACGCGCCCATGATGGGGATCATAAAGTTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCa                         <  2:13248/148‑1 (MQ=255)
                                                                                                                                   aaTCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAgcgc                    >  2:1847649/1‑150 (MQ=255)
                                                                                                                                     tCGACGTGGAAATATACGCCCCCATGATGGGTATCATAAAGGTGAAGCCATGGTAACCTATGGATTCAACTACCCCGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGcca                  <  1:4051599/150‑1 (MQ=255)
                                                                                                                                            ggCAATATAAGCGCCCATGATGGGGACCATAAAGGTTAAGCCAAGGTAACCTATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGCCACCAGGAt           <  2:1295454/150‑1 (MQ=255)
                                                                                                                                              cAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGCCACCAGGATac         >  2:64180/1‑150 (MQ=255)
                                                                                                                                                 tataCGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGCCACCAGGATaccac      <  1:3060795/150‑1 (MQ=255)
                                                                                                                                                  ataAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGCCACCAGGATaccacc     >  2:4407842/1‑150 (MQ=255)
                                                                                                                                                  ataAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGCCACCAGGATaccacc     >  2:4738549/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGCCACCAGGATACCACCGCC  >  CP009273/4084996‑4085294

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑