Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,107,658 C→A 28.4% E51* (GAA→TAA)  wcaN ← putative regulatory subunit for GalU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,107,6580CA28.4% 596.9 / 203.1 271E51* (GAA→TAA) wcaNputative regulatory subunit for GalU
Reads supporting (aligned to +/- strand):  ref base C (72/122);  new base A (40/37);  total (112/159)
Fisher's exact test for biased strand distribution p-value = 2.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 273 total aligned reads displayed.

CATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTC  >  CP009273/2107509‑2107806
                                                                                                                                                     |                                                                                                                                                    
cATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTc                                                                                                                                                      >  2:3570314/1‑150 (MQ=255)
  tAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCtt                                                                                                                                                    <  2:1491091/150‑1 (MQ=255)
   aaTGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTCCTTGGACGCGTGAGTTACCAGGAGGATTTCttt                                                                                                                                                   <  1:701726/150‑1 (MQ=255)
    aTGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTg                                                                                                                                                  >  1:291026/1‑150 (MQ=255)
       gTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATc                                                                                                                                               <  1:2464552/150‑1 (MQ=255)
        tCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATcc                                                                                                                                              >  1:2378568/1‑150 (MQ=255)
          aCGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCt                                                                                                                                            <  1:2718232/150‑1 (MQ=255)
            gCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGc                                                                                                                                          <  1:2562067/150‑1 (MQ=255)
                  cggACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCcac                                                                                                                                    >  1:2218299/1‑150 (MQ=255)
                    gACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACaa                                                                                                                                  >  2:1304488/1‑150 (MQ=255)
                       aGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAAtct                                                                                                                               <  2:613448/150‑1 (MQ=255)
                       aGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAAtct                                                                                                                               <  2:613296/150‑1 (MQ=255)
                        gATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAAtctc                                                                                                                              <  2:3716970/150‑1 (MQ=255)
                           ggACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTc                                                                                                                           <  1:1627915/150‑1 (MQ=255)
                           ggACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTc                                                                                                                           >  1:3362544/1‑150 (MQ=255)
                            gACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTca                                                                                                                          <  1:3003250/150‑1 (MQ=255)
                             aCTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTcaa                                                                                                                         >  1:4324730/1‑150 (MQ=255)
                              cTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTcaac                                                                                                                        <  2:3368093/150‑1 (MQ=255)
                              cTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTcaac                                                                                                                        <  1:4293029/150‑1 (MQ=255)
                              cTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTcaac                                                                                                                        <  1:3558164/150‑1 (MQ=255)
                               tGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTcaaca                                                                                                                       <  2:3980089/150‑1 (MQ=255)
                                   cTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGt                                                                                                                   <  2:4020986/150‑1 (MQ=255)
                                      ccGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACt                                                                                                                <  1:570571/150‑1 (MQ=255)
                                        gCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGa                                                                                                              <  2:952381/150‑1 (MQ=255)
                                        gCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGa                                                                                                              <  2:3233203/150‑1 (MQ=255)
                                          cagcagTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAAt                                                                                                            <  1:4441178/150‑1 (MQ=255)
                                           agcagTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                                                                           <  1:1910684/150‑1 (MQ=255)
                                           agcagTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                                                                           <  1:3068635/150‑1 (MQ=255)
                                           agcagTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                                                                           <  2:3415776/150‑1 (MQ=255)
                                              agTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCAtt                                                                                                        >  1:3904097/1‑150 (MQ=255)
                                               gTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTg                                                                                                       >  2:3912707/1‑150 (MQ=255)
                                                   aCGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCtt                                                                                                   <  2:3874629/150‑1 (MQ=255)
                                                    cGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTg                                                                                                  >  1:230262/1‑150 (MQ=255)
                                                     gCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGt                                                                                                 >  1:1316399/1‑150 (MQ=255)
                                                     gCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGt                                                                                                 >  2:3908993/1‑150 (MQ=255)
                                                       ttCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTcg                                                                                               <  1:2128190/150‑1 (MQ=255)
                                                        tCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTcga                                                                                              >  2:4501302/1‑150 (MQ=255)
                                                               ctgctCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTgg                                                                                       <  1:985890/150‑1 (MQ=255)
                                                                  ctCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAg                                                                                    <  1:1093961/150‑1 (MQ=255)
                                                                  ctCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAg                                                                                    >  1:1007155/1‑150 (MQ=255)
                                                                    cAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCa                                                                                  >  2:1148042/1‑150 (MQ=255)
                                                                    cAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCa                                                                                  >  2:4555677/1‑150 (MQ=255)
                                                                     aaGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAt                                                                                 >  2:1506704/1‑150 (MQ=255)
                                                                     aaGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAt                                                                                 >  2:1506339/1‑150 (MQ=255)
                                                                     aaGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGATACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATTTAATGAATAATTGGCTTGTCGACGATTGGTAGCAt                                                                                 >  1:1195172/1‑150 (MQ=255)
                                                                      aGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtc                                                                                >  2:3014535/1‑150 (MQ=255)
                                                                         aGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATctct                                                                             >  1:1661522/1‑150 (MQ=255)
                                                                           tGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCttt                                                                           <  2:4098365/150‑1 (MQ=255)
                                                                            gATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTg                                                                          >  2:1420/1‑150 (MQ=255)
                                                                              ttCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTggg                                                                        <  1:1845188/150‑1 (MQ=255)
                                                                                cTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTa                                                                      <  2:3178934/150‑1 (MQ=255)
                                                                                 tAACTCATAAGAGGTTTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTAt                                                                     <  2:3770505/150‑1 (MQ=255)
                                                                                 tAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTag                                                                     <  2:2507851/150‑2 (MQ=255)
                                                                                  aaCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATc                                                                    >  1:4074984/1‑150 (MQ=255)
                                                                                   aCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCg                                                                   <  1:2323085/150‑1 (MQ=255)
                                                                                     tCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGcc                                                                 <  1:1567380/150‑1 (MQ=255)
                                                                                         aaGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAg                                                             >  1:3386641/1‑150 (MQ=255)
                                                                                          agagGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGt                                                            >  1:781213/1‑150 (MQ=255)
                                                                                             ggTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTggc                                                         <  1:1456346/150‑1 (MQ=255)
                                                                                                   gAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCaa                                                   >  2:2177554/1‑150 (MQ=255)
                                                                                                   gAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCaa                                                   <  1:890690/150‑1 (MQ=255)
                                                                                                   gAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCaa                                                   >  1:584554/1‑150 (MQ=255)
                                                                                                    aaGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAAc                                                  >  1:2824938/1‑150 (MQ=255)
                                                                                                     aGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACa                                                 >  2:76241/1‑150 (MQ=255)
                                                                                                       tGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACata                                               >  2:2606126/1‑150 (MQ=255)
                                                                                                           tttCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCa                                           >  2:4529647/1‑150 (MQ=255)
                                                                                                            ttCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCAt                                          <  2:503172/150‑1 (MQ=255)
                                                                                                            ttCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCAt                                          >  2:4403107/1‑150 (MQ=255)
                                                                                                            ttCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCAt                                          <  1:3834578/150‑1 (MQ=255)
                                                                                                             tCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATc                                         <  2:3396381/150‑1 (MQ=255)
                                                                                                             tCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATc                                         <  2:4202843/150‑1 (MQ=255)
                                                                                                              cGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATcc                                        <  2:2259009/150‑1 (MQ=255)
                                                                                                               gACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATccc                                       <  2:1033623/150‑1 (MQ=255)
                                                                                                                aCCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCg                                      >  2:2481020/1‑150 (MQ=255)
                                                                                                                aCCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCg                                      <  2:1980024/150‑1 (MQ=255)
                                                                                                                 ccGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCga                                     >  2:528629/1‑150 (MQ=255)
                                                                                                                     gTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTTAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGAcc                                 >  2:2804548/1‑150 (MQ=255)
                                                                                                                     gTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGAcc                                 >  2:2950468/1‑150 (MQ=255)
                                                                                                                        cTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGc                              <  2:291026/150‑1 (MQ=255)
                                                                                                                        cTTGGACGCGTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGc                              >  1:2088862/1‑150 (MQ=255)
                                                                                                                         ttGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCt                             <  1:2337222/150‑1 (MQ=255)
                                                                                                                              cgcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAgg                        <  2:1214731/150‑1 (MQ=255)
                                                                                                                               gcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGa                       >  1:3755266/1‑150 (MQ=255)
                                                                                                                                 gTGAGTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAAt                     >  1:1099944/1‑150 (MQ=255)
                                                                                                                                    aGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAAc                  <  1:3984026/150‑1 (MQ=255)
                                                                                                                                     gTTACCAGGAGGATTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACt                 <  1:4161587/150‑1 (MQ=255)
                                                                                                                                         ccAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCtt             <  1:3169944/150‑1 (MQ=255)
                                                                                                                                         ccAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCtt             <  2:2280343/150‑1 (MQ=255)
                                                                                                                                          cAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCttt            <  2:2378741/150‑1 (MQ=255)
                                                                                                                                           aggaggATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCtttt           >  2:1524405/1‑150 (MQ=255)
                                                                                                                                            ggaggaTTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTa          >  1:4037444/1‑150 (MQ=255)
                                                                                                                                             gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTaa         >  1:2775692/1‑150 (MQ=255)
                                                                                                                                               ggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAAt       <  2:2513496/150‑1 (MQ=255)
                                                                                                                                               ggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAAt       <  2:19229/150‑1 (MQ=255)
                                                                                                                                                gaTTTATTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAAtt      <  1:2830853/150‑1 (MQ=255)
                                                                                                                                                 aTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTc     <  2:2126335/150‑1 (MQ=255)
                                                                                                                                                  tttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCg    >  2:2718765/1‑150 (MQ=255)
                                                                                                                                                  tttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCg    >  1:3460594/1‑150 (MQ=255)
                                                                                                                                                    tCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTc  <  1:4628902/150‑1 (MQ=255)
                                                                                                                                                    tCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTc  >  2:1769299/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTC  >  CP009273/2107509‑2107806

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: