Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,994,076 A→T 92.9% L114* (TTA→TAA)  yigE ← DUF2233 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,994,0760AT92.9% 909.6 / 43.3 276L114* (TTA→TAA) yigEDUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base A (6/14);  new base T (146/110);  total (152/124)
Fisher's exact test for biased strand distribution p-value = 3.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

ATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTG  >  CP009273/3993929‑3994225
                                                                                                                                                   |                                                                                                                                                     
     cATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCAc                                                                                                                                                >  2:4208549/1‑150 (MQ=255)
       ttGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCt                                                                                                                                              >  1:1451726/1‑150 (MQ=255)
         ggCCCTGACTGCACCGCAAACTGAATCTCTTTACGGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTc                                                                                                                                            >  2:137798/1‑150 (MQ=255)
          gCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTct                                                                                                                                           <  1:546399/150‑1 (MQ=255)
          gCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTct                                                                                                                                           <  1:2736125/150‑1 (MQ=255)
          gCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTct                                                                                                                                           <  2:489693/150‑1 (MQ=255)
            ccTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgc                                                                                                                                         >  2:952294/1‑150 (MQ=255)
             cTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgct                                                                                                                                        <  1:4016532/150‑1 (MQ=255)
              tGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgctg                                                                                                                                       <  1:3172051/150‑1 (MQ=255)
                aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGAc                                                                                                                                     <  2:2359914/150‑1 (MQ=255)
                 cTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGAcc                                                                                                                                    <  1:3622598/150‑1 (MQ=255)
                  tGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATACCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCg                                                                                                                                   >  2:2315066/1‑150 (MQ=255)
                   gCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGt                                                                                                                                  >  1:1429134/1‑150 (MQ=255)
                    cACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGtt                                                                                                                                 >  2:1119504/1‑150 (MQ=255)
                          aaaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGAt                                                                                                                           <  1:815493/150‑1 (MQ=255)
                          aaaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGAt                                                                                                                           >  2:1982191/1‑150 (MQ=255)
                          aaaCTGAATAACTTTAATGGTTTTGAAGGCATCCAGACGAACGATGACGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAAATCACTTCACCTGAAGCGAGATATTACGACACATACTGATGAACGATTTCTAt                                                                                                                           >  2:2590710/1‑150 (MQ=255)
                                aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGATATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACaa                                                                                                                     >  1:3987265/1‑150 (MQ=255)
                                  tctcTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc                                                                                                                   <  2:1893947/150‑1 (MQ=255)
                                  tctcTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc                                                                                                                   <  1:2839491/150‑1 (MQ=255)
                                   ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCCGACCGTTTTCGATGTACAAAcc                                                                                                                  <  1:2271531/150‑1 (MQ=255)
                                     ctgtACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGa                                                                                                                <  1:703944/147‑1 (MQ=255)
                                       ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc                                                                                                              <  1:4509016/150‑1 (MQ=255)
                                       ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc                                                                                                              >  1:4431258/1‑150 (MQ=255)
                                       ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc                                                                                                              >  2:3203588/1‑150 (MQ=255)
                                        tACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcg                                                                                                             >  1:3223817/1‑150 (MQ=255)
                                         aCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcgg                                                                                                            >  1:152874/1‑150 (MQ=255)
                                         aCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcgg                                                                                                            >  1:152983/1‑150 (MQ=255)
                                          cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc                                                                                                           <  1:1142857/150‑1 (MQ=255)
                                              ttttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAt                                                                                                       <  1:2324332/150‑1 (MQ=255)
                                                 tGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCGTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGc                                                                                                    <  2:3166710/150‑1 (MQ=255)
                                                   aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                                                  >  1:579922/1‑150 (MQ=255)
                                                   aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                                                  >  2:3205849/1‑150 (MQ=255)
                                                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTtcatc                                                                                             <  1:4336663/150‑1 (MQ=255)
                                                           cAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATa                                                                                          >  1:2521367/1‑150 (MQ=255)
                                                            aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAg                                                                                         <  2:2925508/150‑1 (MQ=255)
                                                                gaacgaTGCCGACTTTATCTCCCGCGACATTAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCGTCTGCTGACCGTTTTCGATGTACAAACCTAGCGGCGCATAGCTTTCATCATAGATgc                                                                                     >  1:2922557/1‑150 (MQ=255)
                                                                   cgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgccgc                                                                                  >  1:3053861/1‑150 (MQ=255)
                                                                     aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGccgccg                                                                                >  1:1130492/1‑150 (MQ=255)
                                                                     aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGccgccg                                                                                >  1:3409001/1‑150 (MQ=255)
                                                                     aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGccgccg                                                                                >  1:1005846/1‑150 (MQ=255)
                                                                       gCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGtt                                                                              <  1:1800959/150‑1 (MQ=255)
                                                                         cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCa                                                                            <  1:4088228/150‑1 (MQ=255)
                                                                         cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCa                                                                            <  1:1834541/150‑1 (MQ=255)
                                                                          gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCCTTCATCATAGATGCCGCCGTTCAt                                                                           >  2:3665180/1‑150 (MQ=255)
                                                                            cTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCg                                                                         >  2:1963475/1‑150 (MQ=255)
                                                                            cTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCg                                                                         >  1:3779632/1‑150 (MQ=255)
                                                                               tATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCa                                                                      >  2:4291827/1‑150 (MQ=255)
                                                                                  ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                                   <  1:2636676/150‑1 (MQ=255)
                                                                                  ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                                   <  2:1022733/150‑1 (MQ=255)
                                                                                        cgACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcct                                                             >  1:350745/1‑150 (MQ=255)
                                                                                         gACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctg                                                            <  1:3518909/150‑1 (MQ=255)
                                                                                         gACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctg                                                            >  1:2040945/1‑150 (MQ=255)
                                                                                          aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctga                                                           >  1:588054/1‑150 (MQ=255)
                                                                                             tAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccc                                                        >  1:2434526/1‑150 (MQ=255)
                                                                                             tAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccc                                                        >  2:615060/1‑150 (MQ=255)
                                                                                              aaaaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccct                                                       <  2:3383221/150‑1 (MQ=255)
                                                                                                aaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccctga                                                     <  2:1724270/150‑1 (MQ=255)
                                                                                                     gccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTAtt                                                <  1:3988063/150‑1 (MQ=255)
                                                                                                     gccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGAGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTAtt                                                >  2:3371953/1‑150 (MQ=255)
                                                                                                      ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTa                                               <  2:2191189/150‑1 (MQ=255)
                                                                                                        gccAGGACGGATAAAGAAATTCCCTTCACCTGAAGTGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAat                                             <  1:817803/150‑1 (MQ=255)
                                                                                                        gccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAat                                             <  2:2602884/150‑1 (MQ=255)
                                                                                                         ccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAata                                            >  1:1336177/1‑150 (MQ=255)
                                                                                                          cAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAatat                                           >  1:4398630/1‑150 (MQ=255)
                                                                                                          cAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAatat                                           <  2:514046/150‑1 (MQ=255)
                                                                                                          cAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAatat                                           >  1:2154245/1‑150 (MQ=255)
                                                                                                          cAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAatat                                           >  2:3334858/1‑150 (MQ=255)
                                                                                                           aGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATc                                          <  2:3374564/150‑1 (MQ=255)
                                                                                                               aggATAAAGAAATTCCCTTCACCTGAAGAGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGcc                                      <  2:3119536/149‑1 (MQ=255)
                                                                                                                ggATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCca                                     >  2:3692820/1‑150 (MQ=255)
                                                                                                                   tAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACGTGACCCTGACTATTAATATCCGCcagca                                  >  2:2872902/1‑150 (MQ=255)
                                                                                                                   tAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCcagca                                  <  1:650057/150‑1 (MQ=255)
                                                                                                                    aaaGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCcagcag                                 >  2:626421/1‑150 (MQ=255)
                                                                                                                    aaaGAAATTCCCTTCACCTGAAGAGAGATTTTACGCCACCTTCTGCTGACCGTTTTAGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGTCCCTGACTATTAATATCCGCcagcag                                 <  2:230373/150‑1 (MQ=255)
                                                                                                                     aaGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCaga                                <  2:3806850/150‑1 (MQ=255)
                                                                                                                      aGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCagag                               >  2:2451914/1‑150 (MQ=255)
                                                                                                                          aTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATg                           >  1:3942152/1‑150 (MQ=255)
                                                                                                                          aTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATg                           >  1:3863723/1‑150 (MQ=255)
                                                                                                                            tCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTa                         >  1:1263588/1‑150 (MQ=255)
                                                                                                                             cccTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTaa                        >  2:2030486/1‑150 (MQ=255)
                                                                                                                              ccTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAAc                       <  2:3198415/150‑1 (MQ=255)
                                                                                                                              ccTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAAc                       >  2:4439172/1‑150 (MQ=255)
                                                                                                                                ttCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGt                     >  1:844783/1‑150 (MQ=255)
                                                                                                                                  cACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTc                   >  1:2709344/1‑150 (MQ=255)
                                                                                                                                    ccTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTccc                 <  1:3247443/150‑1 (MQ=255)
                                                                                                                                      tGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCa               <  2:1448418/150‑1 (MQ=255)
                                                                                                                                      tGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCa               >  1:855432/1‑150 (MQ=255)
                                                                                                                                      tGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCa               <  1:3331175/150‑1 (MQ=255)
                                                                                                                                        aaGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACg             >  2:1789513/1‑150 (MQ=255)
                                                                                                                                        aaGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACg             >  2:4040814/1‑150 (MQ=255)
                                                                                                                                          gCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCt           >  1:59764/1‑150 (MQ=255)
                                                                                                                                           cGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCtt          <  1:200888/150‑1 (MQ=255)
                                                                                                                                            gagaTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTc         <  2:2094211/150‑1 (MQ=255)
                                                                                                                                            gagaTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTc         >  1:2389723/1‑150 (MQ=255)
                                                                                                                                            gagaTTTTACGCAACCTTCTGCTGACAGTTTTCGATGTACAAACAGAGAGTAGCATAGATTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTACCAAAGCttt         <  2:4029940/150‑2 (MQ=255)
                                                                                                                                             agaTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCg        <  2:3779632/150‑1 (MQ=255)
                                                                                                                                             agaTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCg        <  2:1988478/150‑1 (MQ=255)
                                                                                                                                              gaTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGc       <  1:3944674/150‑1 (MQ=255)
                                                                                                                                                ttttACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCa     >  1:1475389/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTG  >  CP009273/3993929‑3994225

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