Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,994,076 A→T 93.1% 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.1% 809.3 / 38.0 232L114* (TTA→TAA) yigEDUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base A (9/7);  new base T (128/88);  total (137/95)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 235 total aligned reads displayed.

CCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCA  >  CP009273/3993927‑3994222
                                                                                                                                                     |                                                                                                                                                  
   tCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACg                                                                                                                                                 <  1:1500172/150‑1 (MQ=255)
      aCATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCa                                                                                                                                              <  2:2419687/150‑1 (MQ=255)
       cATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCAc                                                                                                                                             >  2:1902748/1‑150 (MQ=255)
        agtGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCAcc                                                                                                                                            <  2:4202000/148‑1 (MQ=255)
        aTTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCAcc                                                                                                                                            <  2:3913855/150‑1 (MQ=255)
          tGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCtt                                                                                                                                          <  2:1743006/150‑1 (MQ=255)
          tGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCtt                                                                                                                                          >  2:2130539/1‑150 (MQ=255)
           ggCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTc                                                                                                                                         >  2:866408/1‑150 (MQ=255)
             cccTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctg                                                                                                                                       <  2:1660056/150‑1 (MQ=255)
              ccTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGACCGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgc                                                                                                                                      <  1:3385251/150‑1 (MQ=255)
              ccTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgc                                                                                                                                      <  1:543428/150‑1 (MQ=255)
                tGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgctg                                                                                                                                    >  1:188854/1‑150 (MQ=255)
                 gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGa                                                                                                                                   >  2:1935087/1‑150 (MQ=255)
                 gACTGCACCGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGa                                                                                                                                   <  2:1149117/150‑1 (MQ=255)
                  aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGAc                                                                                                                                  >  2:1895591/1‑150 (MQ=255)
                      cACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGtt                                                                                                                              >  2:1598779/1‑150 (MQ=255)
                        ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGtttt                                                                                                                            <  1:1929562/150‑1 (MQ=255)
                         cGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTc                                                                                                                           <  1:1020541/150‑1 (MQ=255)
                             caCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATg                                                                                                                       <  2:4081827/149‑1 (MQ=255)
                               cTGAATCTCGTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTa                                                                                                                     <  1:3091165/150‑1 (MQ=255)
                                 gAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACa                                                                                                                   <  1:1902718/150‑1 (MQ=255)
                                 gAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACa                                                                                                                   >  2:3735907/1‑150 (MQ=255)
                                 gAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACa                                                                                                                   >  2:3949123/1‑150 (MQ=255)
                                    tctcTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc                                                                                                                <  2:2295032/150‑1 (MQ=255)
                                     ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc                                                                                                               <  1:620073/150‑1 (MQ=255)
                                        tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg                                                                                                            >  1:137314/1‑150 (MQ=255)
                                         ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc                                                                                                           >  2:727135/1‑150 (MQ=255)
                                          tACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcg                                                                                                          >  2:1663602/1‑150 (MQ=255)
                                          tACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcg                                                                                                          >  1:2642333/1‑150 (MQ=255)
                                            cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc                                                                                                        >  1:383921/1‑150 (MQ=255)
                                             tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggcg                                                                                                       <  1:1316727/150‑1 (MQ=255)
                                              ggTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGgcgc                                                                                                      >  2:2506558/1‑150 (MQ=255)
                                              ggTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGgcgc                                                                                                      >  2:3212862/1‑150 (MQ=255)
                                              ggTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGgcgc                                                                                                      >  2:265808/1‑150 (MQ=255)
                                                 tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATa                                                                                                   >  1:4269623/1‑150 (MQ=255)
                                                 tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATa                                                                                                   <  1:2016749/150‑1 (MQ=255)
                                                  ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAg                                                                                                  >  1:3238883/1‑150 (MQ=255)
                                                  ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAg                                                                                                  >  1:1669851/1‑150 (MQ=255)
                                                  ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAg                                                                                                  <  1:2090787/150‑1 (MQ=255)
                                                    gAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCt                                                                                                <  2:1220466/150‑1 (MQ=255)
                                                     aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                                               <  2:1527867/150‑1 (MQ=255)
                                                     aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                                               >  1:1201505/1‑150 (MQ=255)
                                                     aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                                               >  1:1993763/1‑150 (MQ=255)
                                                     aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                                               >  2:192168/1‑150 (MQ=255)
                                                     aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                                               >  1:1200941/1‑150 (MQ=255)
                                                        gCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTtca                                                                                            <  1:1085796/150‑1 (MQ=255)
                                                         cATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTtcat                                                                                           <  1:1174563/150‑1 (MQ=255)
                                                          aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTtcatc                                                                                          >  2:1025380/1‑150 (MQ=255)
                                                           tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTcatcat                                                                                        >  1:3120875/1‑150 (MQ=255)
                                                             cAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATa                                                                                       <  2:2991217/150‑1 (MQ=255)
                                                             cAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATa                                                                                       <  1:887861/150‑1 (MQ=255)
                                                                    acgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgccg                                                                                <  1:1120322/150‑1 (MQ=255)
                                                                    acgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgccg                                                                                >  1:3167693/1‑150 (MQ=255)
                                                                     cgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgccgc                                                                               <  1:2242061/150‑1 (MQ=255)
                                                                       aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGccgccg                                                                             >  2:2684041/1‑150 (MQ=255)
                                                                       aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGccgccg                                                                             >  2:378947/1‑150 (MQ=255)
                                                                         gCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGtt                                                                           <  1:2097071/150‑1 (MQ=255)
                                                                          ccGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTc                                                                          <  1:2176187/150‑1 (MQ=255)
                                                                           cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCa                                                                         <  2:1971583/150‑1 (MQ=255)
                                                                             aCTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATc                                                                       >  2:3055418/1‑150 (MQ=255)
                                                                                   tctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATc                                                                 <  1:41161/150‑1 (MQ=255)
                                                                                    ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                                <  1:3766028/150‑1 (MQ=255)
                                                                                    ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                                >  1:418758/1‑150 (MQ=255)
                                                                                       ccGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCa                                                             >  1:1582184/1‑150 (MQ=255)
                                                                                          cgACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcct                                                          >  2:190052/1‑150 (MQ=255)
                                                                                           gACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctg                                                         <  1:1035323/150‑1 (MQ=255)
                                                                                            aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctga                                                        >  1:2974411/1‑150 (MQ=255)
                                                                                            aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctga                                                        >  1:800519/1‑150 (MQ=255)
                                                                                              aTAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgacc                                                      >  1:1796318/1‑150 (MQ=255)
                                                                                               tAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccc                                                     >  2:2508963/1‑150 (MQ=255)
                                                                                                aaaaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccct                                                    <  2:712530/150‑1 (MQ=255)
                                                                                                  aaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccctga                                                  >  1:39457/1‑150 (MQ=255)
                                                                                                   acacGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccctgac                                                 >  1:2374032/1‑150 (MQ=255)
                                                                                                    cacGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACt                                                >  1:2411349/1‑150 (MQ=255)
                                                                                                           ccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAata                                         <  1:1450825/150‑1 (MQ=255)
                                                                                                             aGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATc                                       >  1:3787138/1‑150 (MQ=255)
                                                                                                             aGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAGTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATc                                       >  2:1408425/1‑150 (MQ=255)
                                                                                                              ggACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATcc                                      >  1:1573532/1‑150 (MQ=255)
                                                                                                               gACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCg                                     >  2:814415/1‑150 (MQ=255)
                                                                                                                 cGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGcc                                   <  2:1666763/150‑1 (MQ=255)
                                                                                                                 cGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGcc                                   <  2:1079104/150‑1 (MQ=255)
                                                                                                                  ggATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCca                                  >  2:2805174/1‑150 (MQ=255)
                                                                                                                  ggATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCca                                  >  2:634123/1‑150 (MQ=255)
                                                                                                                  ggATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCca                                  >  1:2817315/1‑150 (MQ=255)
                                                                                                                      aaaGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCcagcag                              >  1:490159/1‑150 (MQ=255)
                                                                                                                         gAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGTGc                           >  2:3435773/1‑150 (MQ=255)
                                                                                                                          aaaTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCa                          >  2:1268510/1‑150 (MQ=255)
                                                                                                                          aaaTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTGCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCa                          >  2:1504007/1‑150 (MQ=255)
                                                                                                                           aaTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCAt                         >  1:850685/1‑150 (MQ=255)
                                                                                                                           aaTTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCAt                         >  1:3518714/1‑150 (MQ=255)
                                                                                                                                 cTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACg                   >  1:1054505/1‑150 (MQ=255)
                                                                                                                                        tGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCa            >  2:3257788/1‑150 (MQ=255)
                                                                                                                                        tGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCa            >  1:3296498/1‑150 (MQ=255)
                                                                                                                                        tGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCa            >  2:258952/1‑150 (MQ=255)
                                                                                                                                         gAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCAc           <  1:3930726/150‑1 (MQ=255)
                                                                                                                                          aaGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACg          >  1:995324/1‑150 (MQ=255)
                                                                                                                                              gagaTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTc      >  2:3107746/1‑150 (MQ=255)
                                                                                                                                               agaTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCg     <  2:44881/150‑1 (MQ=255)
                                                                                                                                                 aTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGcc   >  1:2275978/1‑150 (MQ=255)
                                                                                                                                                  ttttACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCa  <  2:3120273/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCA  >  CP009273/3993927‑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: