Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,073,654 T→C 5.9% intergenic (‑148/+476) mdtA ← / ← PP_3586 multidrug transporter membrane fusion protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,073,6540TC 5.9% 263.0 / 4.4 85intergenic (‑148/+476)mdtA/PP_3586multidrug transporter membrane fusion protein/transposase
Reads supporting (aligned to +/- strand):  ref base T (40/40);  new base C (3/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

TTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCG  >  NC_002947/4073510‑4073798
                                                                                                                                                |                                                                                                                                                
ttGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGAGGGTGCCAATAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTcc                                                                                                                                             >  2:1419116/1‑150 (MQ=255)
          cTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTAc                                                                                                                                   >  1:1091639/1‑150 (MQ=255)
          cTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTAc                                                                                                                                   >  2:1513853/1‑150 (MQ=255)
          cTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTAc                                                                                                                                   <  2:465991/150‑1 (MQ=255)
             gTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAgg                                                                                                                                <  2:73278/150‑1 (MQ=255)
             gTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAgg                                                                                                                                <  1:1326214/150‑1 (MQ=255)
                cGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgtg                                                                                                                             >  2:1386823/1‑150 (MQ=255)
                 ggAATGACGACTTGAGGAGGATGAACGATAAGCATTTTCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCGTAACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgtgt                                                                                                                            >  1:544365/1‑150 (MQ=255)
                  gAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGACTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGtt                                                                                                                           >  1:1020701/1‑150 (MQ=255)
                      gacgacTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCtt                                                                                                                       <  2:1020701/150‑1 (MQ=255)
                         gacTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAAGGAAAACAGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGa                                                                                                                    >  2:1118562/1‑150 (MQ=255)
                         gacTTGAGGAGGATGAACGATAAGCACTGCAATTGTTTCAGGGAATGGCCTGTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGa                                                                                                                    <  1:1184479/150‑1 (MQ=255)
                            ttGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCg                                                                                                                 >  2:506841/1‑150 (MQ=255)
                                ggaggaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTgg                                                                                                             >  2:1217688/1‑150 (MQ=255)
                                ggaggaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATAGTGTCTACCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTgg                                                                                                             >  2:1277599/1‑150 (MQ=255)
                                    gaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgt                                                                                                                              <  1:1278748/129‑1 (MQ=255)
                                    gaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgt                                                                                                                              >  2:1278748/1‑129 (MQ=255)
                                    gaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGtt                                                                                                                           >  2:1135168/1‑132 (MQ=255)
                                    gaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGtt                                                                                                                           <  1:1135168/132‑1 (MQ=255)
                                    gaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCggg                                                                                                         <  1:1547730/150‑1 (MQ=255)
                                      tGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCgg                                                                                                                                           >  2:627709/1‑114 (MQ=255)
                                      tGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCgg                                                                                                                                           <  1:627709/114‑1 (MQ=255)
                                      tGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGcc                                                                                                       <  1:1513853/150‑1 (MQ=255)
                                      tGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCCGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGcc                                                                                                       >  2:544125/1‑150 (MQ=255)
                                       gAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAgg                                                                                                                                >  1:315833/1‑124 (MQ=255)
                                       gAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAgg                                                                                                                                <  2:315833/124‑1 (MQ=255)
                                           gATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAg                                                                                                  <  1:544125/150‑1 (MQ=255)
                                             tAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACACGTTCCGGGCCGttt                                                                                                                                    >  2:62940/1‑113 (MQ=255)
                                             tAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACACGTTCCGGGCCGttt                                                                                                                                    <  1:62940/114‑2 (MQ=255)
                                              aaGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGtt                                                                                               <  1:1217688/150‑1 (MQ=255)
                                                      ccATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAgg                                                                                                                                >  2:625262/1‑109 (MQ=255)
                                                      ccATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAgg                                                                                                                                <  1:625262/109‑1 (MQ=255)
                                                      ccATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTg                                                                                                                      <  2:579522/119‑1 (MQ=255)
                                                      ccATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTg                                                                                                                      >  1:579522/1‑119 (MQ=255)
                                                      ccATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTgac                                                                                       >  2:636606/1‑148 (MQ=255)
                                                       cATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGATTTGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCt                                                                                      <  1:1621436/150‑1 (MQ=255)
                                                                 gCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTg                                                                                                              <  1:373606/116‑1 (MQ=255)
                                                                 gCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTg                                                                                                              >  2:373606/1‑116 (MQ=255)
                                                                   aaaGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCtt                                                                          >  2:1457167/1‑150 (MQ=255)
                                                                   aaaGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCtt                                                                          >  2:124434/1‑150 (MQ=255)
                                                                    aaGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCttt                                                                         >  1:57241/1‑150 (MQ=255)
                                                                       ggCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTcgc                                                                      <  2:1209600/150‑1 (MQ=255)
                                                                          cTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTCGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCggg                                                                   <  1:1419116/150‑1 (MQ=255)
                                                                           tttACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGt                                                                  >  2:963788/1‑150 (MQ=255)
                                                                           tttACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGt                                                                  <  2:1538591/150‑1 (MQ=255)
                                                                               ccAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGt                                                                                   <  1:237932/129‑1 (MQ=255)
                                                                               ccAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGt                                                                                   >  2:237932/1‑129 (MQ=255)
                                                                                cAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGt                                                                                   <  1:1395449/128‑1 (MQ=255)
                                                                                cAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGt                                                                                   >  2:1395449/1‑128 (MQ=255)
                                                                                     gTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCg                                                                                                                 <  1:1506019/93‑1 (MQ=255)
                                                                                     gTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCg                                                                                                                 >  2:1506019/1‑93 (MQ=255)
                                                                                     gTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCtt                                                        >  2:902287/1‑150 (MQ=255)
                                                                                       ttACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTTCACCCGcgaaa                                                      <  2:98701/150‑5 (MQ=255)
                                                                                        tACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTa                                                     <  2:57241/150‑1 (MQ=255)
                                                                                        tACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACACGTTCCGGGcc                                                                                                                                        >  1:494885/1‑67 (MQ=38)
                                                                                        tACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACACGTTCCGGGcc                                                                                                                                        <  2:494885/67‑1 (MQ=38)
                                                                                         aCCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTAt                                                    >  1:1116559/1‑150 (MQ=255)
                                                                                          ccTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATc                                                   <  2:331338/150‑1 (MQ=255)
                                                                                           cTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGc                                                                                        <  1:37613/112‑1 (MQ=255)
                                                                                           cTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGc                                                                                        >  2:37613/1‑112 (MQ=255)
                                                                                             ttGGCTGAGGTTGCCAAGAGGCGGAATTTTAAAGGAAGTCGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTGGCTGGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAGCCCGCTTCTATcaaa                                                <  2:602168/150‑1 (MQ=255)
                                                                                                ccTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaacaa                                             >  2:1567633/1‑150 (MQ=255)
                                                                                                 cTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGTGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaacaaa                                            >  2:432124/1‑150 (MQ=255)
                                                                                                  tGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACaaaa                                           <  1:636606/150‑1 (MQ=255)
                                                                                                   gACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAg                                          <  1:963788/150‑1 (MQ=255)
                                                                                                    aCGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaa                                                <  2:770955/143‑1 (MQ=255)
                                                                                                    aCGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaa                                                >  1:770955/1‑143 (MQ=255)
                                                                                                    aCGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGa                                         >  2:1220054/1‑150 (MQ=255)
                                                                                                     cGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGAt                                        <  1:1220054/150‑1 (MQ=255)
                                                                                                      gTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGtt                                                                                          <  2:1436136/99‑1 (MQ=255)
                                                                                                      gTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGtt                                                                                          >  1:1436136/1‑99 (MQ=255)
                                                                                                       ttGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTTCCTGGGCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATaa                                      <  2:597181/150‑1 (MQ=255)
                                                                                                        tGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAAc                                     >  1:677270/1‑150 (MQ=255)
                                                                                                         gCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAAcc                                    <  1:1386823/150‑1 (MQ=255)
                                                                                                           cAAGAGGCGGAATTTTAAATGAACACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATCACCCa                                  >  2:175167/1‑150 (MQ=255)
                                                                                                           cAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGt                                                                                   >  2:40268/1‑101 (MQ=255)
                                                                                                           cAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGt                                                                                   <  1:40268/101‑1 (MQ=255)
                                                                                                            aaGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCAc                                 >  2:160711/1‑150 (MQ=255)
                                                                                                                   ggAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGttctt                          <  1:506841/150‑1 (MQ=255)
                                                                                                                     aaTTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGttctttc                        <  1:902287/150‑1 (MQ=255)
                                                                                                                          tAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCg                   <  2:1091639/150‑1 (MQ=255)
                                                                                                                                  aaCGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAAc           >  1:58584/1‑150 (MQ=255)
                                                                                                                                       cccGGAACACGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGgcg      >  1:768586/1‑150 (MQ=255)
                                                                                                                                        ccGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTTCACCCGcgaaa                                                      >  2:1523173/1‑97 (MQ=255)
                                                                                                                                        ccGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTTCACCCGcgaaa                                                      <  1:1523173/101‑5 (MQ=255)
                                                                                                                                               gtgtTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTACCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGcgcgcg  >  1:881157/1‑146 (MQ=255)
                                                                                                                                               gtgtTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTACCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGcgcgcg  <  2:881157/146‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCG  >  NC_002947/4073510‑4073798

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