Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,141,654 G→C 5.2% P905R (CCG→CGG)  rne ← ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,141,6540GC 5.2% 255.7 / 2.4 75P905R (CCG→CGG) rneribonuclease E
Reads supporting (aligned to +/- strand):  ref base G (41/30);  new base C (2/2);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACC  >  NC_000913/1141505‑1141803
                                                                                                                                                     |                                                                                                                                                     
ccGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCg                                                                                                                                                       >  2:163465/1‑150 (MQ=255)
ccGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCg                                                                                                                                                       <  1:758991/150‑1 (MQ=255)
      ctTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCgg                                                                                                                                                 >  1:1053754/1‑150 (MQ=255)
           tGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAAcc                                                                                                                                            >  2:402779/1‑150 (MQ=255)
           tGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAAcc                                                                                                                                            <  2:1239199/150‑1 (MQ=255)
            gCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCt                                                                                                                                           <  1:187342/150‑1 (MQ=255)
             cGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTg                                                                                                                                          <  1:1216677/150‑1 (MQ=255)
             cGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGAACACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCCGAACCGCAAccgg                                                                                                                                          <  1:72569/150‑3 (MQ=255)
                   aaCCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggg                                                                                                                                    >  2:837537/1‑150 (MQ=255)
                   aaCCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggg                                                                                                                                    <  1:1125952/150‑1 (MQ=255)
                    accaccGGTTCTGCTTGTTCTGCAACTTCCTGGGATACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCACCCACTTCCGGTTGCGGTTCGGCAACCTGCACgggg                                                                                                                                   >  2:1428645/1‑150 (MQ=255)
                     ccaccGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggggg                                                                                                                                  >  2:1393552/1‑150 (MQ=255)
                     ccaccGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCAGGTTGCCGAACCGCAACCGGaagtggtt                                                                                                                                  <  1:1880530/150‑7 (MQ=255)
                       accGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCt                                                                                                                                <  1:1060637/150‑1 (MQ=255)
                        ccGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCtt                                                                                                                               >  1:148156/1‑150 (MQ=255)
                             tCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACg                                                                                                                          <  1:1352414/150‑1 (MQ=255)
                                         gCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCAct                                                                                                              <  1:393977/150‑1 (MQ=255)
                                         gCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCAct                                                                                                              >  2:131762/1‑150 (MQ=255)
                                         gCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCAct                                                                                                              >  2:13769/1‑150 (MQ=255)
                                          cAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCActt                                                                                                             >  1:962768/1‑150 (MQ=255)
                                             cTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCActtctt                                                                                                          >  2:141559/1‑150 (MQ=255)
                                                   gggCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTcaacaa                                                                                                    >  2:37731/1‑150 (MQ=255)
                                                   gggCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTcaacaa                                                                                                    >  1:286968/1‑150 (MQ=255)
                                                        aCGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGgc                                                                                               >  1:606324/1‑150 (MQ=255)
                                                         cGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGgcg                                                                                              <  2:1053754/150‑1 (MQ=255)
                                                         cGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGgcg                                                                                              >  1:1668537/1‑150 (MQ=255)
                                                          ggCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGgcgc                                                                                             >  2:1393634/1‑150 (MQ=255)
                                                           gCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGgcgcg                                                                                            >  1:112815/1‑150 (MQ=255)
                                                           gCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTATACGACACCGGCAAGTTGTTCAACAAATGgcggg                                                                                            <  2:1792956/150‑3 (MQ=255)
                                                           gCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGGCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGgcgcg                                                                                            <  1:241238/150‑1 (MQ=255)
                                                               cATCAGACTCGCTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTaa                                                                                        >  2:1671037/1‑150 (MQ=255)
                                                                 tCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTaaca                                                                                      >  2:161547/1‑150 (MQ=255)
                                                                  cAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTaacaa                                                                                     >  2:1221935/1‑150 (MQ=255)
                                                                      cTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCgg                                                                                 <  1:402779/150‑1 (MQ=255)
                                                                        cGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGtt                                                                               >  2:38678/1‑150 (MQ=255)
                                                                        cGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGtt                                                                               >  1:971789/1‑150 (MQ=255)
                                                                          gTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTct                                                                             >  1:1311376/1‑149 (MQ=255)
                                                                          gTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCg                                                                             <  1:837537/150‑1 (MQ=255)
                                                                           tAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGa                                                                            <  1:38678/150‑1 (MQ=255)
                                                                              tCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAg                                                                         <  1:131762/150‑1 (MQ=255)
                                                                                 ccTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAgcgg                                                                      <  2:971789/150‑1 (MQ=255)
                                                                                    gCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAgcggcgg                                                                   >  1:450076/1‑150 (MQ=255)
                                                                                     cGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGGCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAgcggcggc                                                                  <  2:424519/150‑1 (MQ=255)
                                                                                        cTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGAc                                                               >  2:676220/1‑150 (MQ=255)
                                                                                             cAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGa                                                          <  1:1070243/150‑1 (MQ=255)
                                                                                             cAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAAAAACCTGTTCGATAGCGGCGGCGACAGGGa                                                          <  2:962768/150‑1 (MQ=255)
                                                                                              aGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGAc                                                         <  2:297391/150‑1 (MQ=255)
                                                                                              aGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGAc                                                         >  1:220601/1‑150 (MQ=255)
                                                                                              aGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGAc                                                         >  1:615951/1‑150 (MQ=255)
                                                                                                 tACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTc                                                      <  2:148156/150‑1 (MQ=255)
                                                                                                     gcgGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCCCTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTc                                                  >  1:1407890/1‑150 (MQ=255)
                                                                                                     gcgGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTc                                                  <  1:141559/150‑1 (MQ=255)
                                                                                                      cgGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCa                                                 <  2:1668537/150‑1 (MQ=255)
                                                                                                       ggCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCAc                                                <  2:112815/150‑1 (MQ=255)
                                                                                                       ggCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCAc                                                >  2:1484892/1‑150 (MQ=255)
                                                                                                           gCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATc                                            >  1:753596/1‑150 (MQ=255)
                                                                                                           gCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATc                                            <  2:606324/150‑1 (MQ=255)
                                                                                                              aTCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACGCCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGc                                         >  2:1990611/1‑150 (MQ=255)
                                                                                                              aTCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGc                                         <  1:676220/150‑1 (MQ=255)
                                                                                                                 aCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCCCCATCGGCTGc                                      >  2:565025/1‑150 (MQ=255)
                                                                                                                  cTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGca                                     >  1:1090317/1‑150 (MQ=255)
                                                                                                                  cTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGca                                     >  1:1097539/1‑150 (MQ=255)
                                                                                                                    tCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGcaca                                   <  1:465952/150‑1 (MQ=255)
                                                                                                                             gtcgtTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTActtct                          >  2:1790619/1‑150 (MQ=255)
                                                                                                                                gtTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACTACCGGTTCGATAGGGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTcc                       <  1:1671037/150‑1 (MQ=255)
                                                                                                                                    cAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGtt                   >  2:728431/1‑150 (MQ=255)
                                                                                                                                    cAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGtt                   >  2:790222/1‑150 (MQ=255)
                                                                                                                                     aaCCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTc                  >  2:352348/1‑150 (MQ=255)
                                                                                                                                        cACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGtt                   >  2:1828392/1‑146 (MQ=255)
                                                                                                                                        cACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGtt                   <  1:1828392/146‑1 (MQ=255)
                                                                                                                                         cgttCAGGTTGCCGAACCGCAACCGGAACGGGGGCTTCTACGACCCCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGCCCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCgcgc              >  1:503279/3‑150 (MQ=255)
                                                                                                                                         aCTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGt                    <  1:578232/144‑1 (MQ=255)
                                                                                                                                         aCTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGt                    >  2:578232/1‑144 (MQ=255)
                                                                                                                                          gttCAGGTTGCCGAACCGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCgct             >  2:1874775/2‑150 (MQ=255)
                                                                                                                                              cGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGt                    <  1:1225601/139‑1 (MQ=255)
                                                                                                                                              cGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGt                    >  2:1225601/1‑139 (MQ=255)
                                                                                                                                                   gCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCaa    >  2:919976/1‑150 (MQ=255)
                                                                                                                                                     ggTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAAcc  >  2:1729451/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACC  >  NC_000913/1141505‑1141803

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