Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,994,076 A→T 93.4% 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,0760AT93.4% 894.2 / 40.4 258L114* (TTA→TAA) yigEDUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base A (7/10);  new base T (126/115);  total (133/125)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 260 total aligned reads displayed.

CATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCA  >  CP009273/3993928‑3994222
                                                                                                                                                    |                                                                                                                                                  
  tCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACg                                                                                                                                                 <  2:2298619/150‑1 (MQ=255)
     aCATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCCGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCa                                                                                                                                              <  1:1028073/150‑1 (MQ=255)
      cATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCAc                                                                                                                                             >  1:3645178/1‑150 (MQ=255)
       aTTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCAcc                                                                                                                                            >  2:4349604/1‑150 (MQ=255)
        ttGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCt                                                                                                                                           <  2:910733/150‑1 (MQ=255)
         tGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCtt                                                                                                                                          <  2:756422/150‑1 (MQ=255)
            cccTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctg                                                                                                                                       >  2:2920204/1‑150 (MQ=255)
             ccTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgc                                                                                                                                      >  2:4517564/1‑150 (MQ=255)
                 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGAc                                                                                                                                  <  1:3182854/150‑1 (MQ=255)
                 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGAc                                                                                                                                  >  2:2199427/1‑150 (MQ=255)
                   tGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCg                                                                                                                                >  1:3112929/1‑150 (MQ=255)
                      aCCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGttt                                                                                                                             <  1:4288919/150‑1 (MQ=255)
                      aCCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGttt                                                                                                                             >  1:2621384/1‑150 (MQ=255)
                       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGtttt                                                                                                                            <  1:2220071/150‑1 (MQ=255)
                       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGtttt                                                                                                                            <  2:597469/150‑1 (MQ=255)
                        cGCAAACTGAATCTCTGTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTc                                                                                                                           <  2:3966577/150‑1 (MQ=255)
                         gCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCg                                                                                                                          <  2:2519112/150‑1 (MQ=255)
                          cAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGa                                                                                                                         >  1:701641/1‑150 (MQ=255)
                             aCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGt                                                                                                                      >  1:809282/1‑150 (MQ=255)
                             aCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGt                                                                                                                      <  2:2787010/150‑1 (MQ=255)
                              cTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTa                                                                                                                     <  1:4048402/150‑1 (MQ=255)
                               ttaaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTAc                                                                                                                    <  1:3365900/148‑1 (MQ=255)
                                gAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACa                                                                                                                   >  1:1801145/1‑150 (MQ=255)
                                  aTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa                                                                                                                 >  1:3923462/1‑150 (MQ=255)
                                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg                                                                                                              >  1:2563775/1‑150 (MQ=255)
                                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg                                                                                                            <  1:1882579/150‑1 (MQ=255)
                                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg                                                                                                            <  1:3395779/150‑1 (MQ=255)
                                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg                                                                                                            <  1:1160023/150‑1 (MQ=255)
                                           cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc                                                                                                        <  1:2783083/150‑1 (MQ=255)
                                              gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa                                                                                                     >  1:827978/1‑150 (MQ=255)
                                                tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATa                                                                                                   <  1:545959/150‑1 (MQ=255)
                                                tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATa                                                                                                   >  1:404827/1‑150 (MQ=255)
                                                 ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAg                                                                                                  <  1:75077/150‑1 (MQ=255)
                                                  tGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGc                                                                                                 >  1:2803109/1‑150 (MQ=255)
                                                    aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                                               >  1:3278407/1‑150 (MQ=255)
                                                    aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                                               <  1:4784668/150‑1 (MQ=255)
                                                        cATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTtcat                                                                                           <  1:4349604/150‑1 (MQ=255)
                                                            cAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATa                                                                                       <  2:2946982/150‑1 (MQ=255)
                                                               acgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGAt                                                                                    <  2:1701603/150‑1 (MQ=255)
                                                                cgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATg                                                                                   >  1:1233143/1‑150 (MQ=255)
                                                                cgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATg                                                                                   >  1:3001478/1‑150 (MQ=255)
                                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgcc                                                                                 <  2:3730910/150‑1 (MQ=255)
                                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGAGgcc                                                                                 <  2:4788189/150‑1 (MQ=255)
                                                                   acgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgccg                                                                                >  2:1337133/1‑150 (MQ=255)
                                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgccgcc                                                                              >  2:4040621/1‑150 (MQ=255)
                                                                      aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGccgccg                                                                             <  2:4074198/150‑1 (MQ=255)
                                                                         ccGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTc                                                                          <  2:2197061/150‑1 (MQ=255)
                                                                         ccGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTc                                                                          <  2:701547/150‑1 (MQ=255)
                                                                         ccGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTc                                                                          >  1:3433668/1‑150 (MQ=255)
                                                                          cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCa                                                                         >  1:2413337/1‑150 (MQ=255)
                                                                          cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCa                                                                         >  2:1377033/1‑150 (MQ=255)
                                                                           gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCAt                                                                        <  2:2447538/150‑1 (MQ=255)
                                                                            aCTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATc                                                                       <  1:2517169/150‑1 (MQ=255)
                                                                             cTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCg                                                                      >  2:674742/1‑150 (MQ=255)
                                                                               ttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGcc                                                                    <  2:2645853/150‑1 (MQ=255)
                                                                               ttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGcc                                                                    <  2:2821035/150‑1 (MQ=255)
                                                                                  tctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATc                                                                 <  1:199651/150‑1 (MQ=255)
                                                                                   ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                                <  1:1198157/150‑1 (MQ=255)
                                                                                   ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                                <  2:1789784/150‑1 (MQ=255)
                                                                                   ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                                >  2:2081825/1‑150 (MQ=255)
                                                                                   ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                                <  2:4193675/150‑1 (MQ=255)
                                                                                      ccGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCa                                                             <  1:1005759/150‑1 (MQ=255)
                                                                                        gcgACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcc                                                           <  2:1006436/150‑1 (MQ=255)
                                                                                          gACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctg                                                         <  1:1352173/150‑1 (MQ=255)
                                                                                           aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctga                                                        <  2:4006041/150‑1 (MQ=255)
                                                                                                    acGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTa                                               <  2:1617366/150‑1 (MQ=255)
                                                                                                    acGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTa                                               <  2:4262482/150‑1 (MQ=255)
                                                                                                     cgccgcCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTAt                                              >  1:2167158/1‑150 (MQ=255)
                                                                                                     cgccgcCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTAt                                              >  2:1027664/1‑150 (MQ=255)
                                                                                                       ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTa                                            >  2:35954/1‑150 (MQ=255)
                                                                                                        cgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTaa                                           >  1:1762045/1‑150 (MQ=255)
                                                                                                         gccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAat                                          <  1:2518746/150‑1 (MQ=255)
                                                                                                           cAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAatat                                        <  2:4551626/150‑1 (MQ=255)
                                                                                                           cAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAatat                                        <  1:2434279/150‑1 (MQ=255)
                                                                                                           cAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAatat                                        <  1:3577860/150‑1 (MQ=255)
                                                                                                             ggACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATcc                                      >  1:626482/1‑150 (MQ=255)
                                                                                                             ggACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATcc                                      >  1:2844480/1‑150 (MQ=255)
                                                                                                                cGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGcc                                   <  2:1916606/150‑1 (MQ=255)
                                                                                                                cGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGcc                                   >  1:4481222/1‑150 (MQ=255)
                                                                                                                    tAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCcagca                               >  1:2249621/1‑150 (MQ=255)
                                                                                                                     aaaGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCcagcag                              >  2:4672247/1‑150 (MQ=255)
                                                                                                                      aaGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCaga                             >  1:579696/1‑150 (MQ=255)
                                                                                                                         aaaTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCa                          >  2:3306465/1‑150 (MQ=255)
                                                                                                                          aaTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCAt                         <  2:3853967/150‑1 (MQ=255)
                                                                                                                          aaTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCAt                         <  2:1445283/150‑1 (MQ=255)
                                                                                                                           aTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATg                        <  2:4124003/150‑1 (MQ=255)
                                                                                                                            ttCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGt                       <  2:3030109/150‑1 (MQ=255)
                                                                                                                            ttCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGt                       >  2:839265/1‑150 (MQ=255)
                                                                                                                            ttCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGt                       <  2:1509882/150‑1 (MQ=255)
                                                                                                                              cccTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTaa                     >  1:1409758/1‑150 (MQ=255)
                                                                                                                               ccTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAAc                    <  2:4598254/150‑1 (MQ=255)
                                                                                                                                 ttCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGt                  >  2:325987/1‑150 (MQ=255)
                                                                                                                                      cTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTcccc             <  2:2172449/150‑1 (MQ=255)
                                                                                                                                        gAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCAc           <  2:1233143/150‑1 (MQ=255)
                                                                                                                                         aaGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACg          <  2:4602984/150‑1 (MQ=255)
                                                                                                                                         aaGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACg          >  2:2812441/1‑150 (MQ=255)
                                                                                                                                             gagaTTTTACGCCACCTTCTGCTGACCGTTGTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTc      <  1:194161/150‑1 (MQ=255)
                                                                                                                                              agaTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCg     >  1:4708283/1‑150 (MQ=255)
                                                                                                                                              agaTTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCg     <  1:905742/150‑1 (MQ=255)
                                                                                                                                               gaTTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGc    >  1:2730638/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCA  >  CP009273/3993928‑3994222

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: