Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,171,856 T→A 5.2% E104V (GAG→GTG)  PP_3670 ← DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,171,8560TA 5.2% 322.2 / 5.6 109E104V (GAG→GTG) PP_3670DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (47/54);  new base A (3/3);  total (51/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 111 total aligned reads displayed.

ACATCGCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATCAGCAACGTG  >  NC_002947/4171707‑4172005
                                                                                                                                                     |                                                                                                                                                     
     gCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGcc                                                                                                                                                  >  1:1815176/1‑150 (MQ=255)
       cccGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAg                                                                                                                                                >  2:350300/1‑150 (MQ=255)
        ccGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGt                                                                                                                                               <  1:611096/150‑1 (MQ=255)
        ccGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGt                                                                                                                                               <  2:1662566/150‑1 (MQ=255)
        ccGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGAGGCCACGGccgtc                                                                                                                                               <  2:110644/150‑4 (MQ=255)
            cAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCa                                                                                                                                           <  1:765851/150‑1 (MQ=255)
            cAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCa                                                                                                                                           <  2:2129876/150‑1 (MQ=255)
             aGGGTGCACCCCCACCAGCACGCCGATACATCCCATCAACACCGCCACTCACTGACCCACTGTCACTCGCTGTGTGAGCAGCGGCGCCGACAGCGCGGTGCCCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAg                                                                                                                                          <  2:216068/150‑1 (MQ=255)
             aGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAg                                                                                                                                          <  1:1921459/150‑1 (MQ=255)
                 tGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTAt                                                                                                                                      >  2:1936643/1‑150 (MQ=255)
                     ccaccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCa                                                                                                                                  <  1:350300/150‑1 (MQ=255)
                      caccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCaa                                                                                                                                 <  1:1846583/150‑1 (MQ=255)
                      caccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCaa                                                                                                                                 >  1:2171102/1‑150 (MQ=255)
                      caccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCaa                                                                                                                                 >  1:2165986/1‑150 (MQ=255)
                      caccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCaa                                                                                                                                 <  1:1155001/150‑1 (MQ=255)
                          accaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCg                                                                                                                             <  1:2146443/150‑1 (MQ=255)
                           ccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTa                                                                                                                                       >  1:1285151/1‑139 (MQ=255)
                                aCGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACGGTCACTCGCTCTTTGAGCAGCGGCGCCGAGAGCGCGGTGACAAGCAAAGGGGCAAGGAACTTGACGCCCGAGTCCTCGGCCAGTGGCAGGTATTGAAAGCCGGTCGTg                                                                                                                       >  2:1331354/1‑150 (MQ=255)
                                 cGCCGATAAAGCCCATCACCACCGACACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGa                                                                                                                      >  1:290941/1‑150 (MQ=255)
                                  gCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGaa                                                                                                                     >  2:1992110/1‑150 (MQ=255)
                                  gCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGAGGCCACGGccgtc                                                                                                                                               <  1:823079/124‑4 (MQ=255)
                                  gCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGAGGCCACGGccgtc                                                                                                                                               >  2:823079/1‑121 (MQ=255)
                                   ccGATAAAGCCCATCACCACCGCCAACCACTGCACGACTTGCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTCTTGCAAGCCGGTCGTGAAc                                                                                                                    >  1:1145733/1‑150 (MQ=255)
                                    cGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAAca                                                                                                                   <  1:221297/150‑1 (MQ=255)
                                    cGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAAca                                                                                                                   <  2:1687714/150‑1 (MQ=255)
                                    cGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAAc                                                                                                                    <  2:2192195/149‑1 (MQ=255)
                                    cGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAAc                                                                                                                    >  1:2192195/1‑149 (MQ=255)
                                      aTAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGAGGTGAAcagc                                                                                                                 >  1:100930/1‑150 (MQ=255)
                                          aGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGc                                                                                                             <  2:2165986/150‑1 (MQ=255)
                                           gCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGTGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCt                                                                                                            >  1:361894/1‑150 (MQ=255)
                                           gCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCt                                                                                                            >  1:1405362/1‑150 (MQ=255)
                                           gCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCt                                                                                                            >  2:613102/1‑150 (MQ=255)
                                           gCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCACCCGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCt                                                                                                            >  1:1240328/1‑150 (MQ=255)
                                            cccATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTg                                                                                                           >  1:337795/1‑150 (MQ=255)
                                            cccATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACAGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTg                                                                                                           <  1:1477122/150‑1 (MQ=255)
                                             ccATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTgg                                                                                                          >  1:2041960/1‑150 (MQ=255)
                                              cATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCATTGCCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGt                                                                                                         >  1:122835/1‑150 (MQ=255)
                                               aTCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGc                                                                                                                               <  1:930880/127‑1 (MQ=255)
                                               aTCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGc                                                                                                                               >  2:930880/1‑127 (MQ=255)
                                               aTCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTg                                                                                                        <  1:1539136/150‑1 (MQ=255)
                                               aTCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTg                                                                                                        >  2:291201/1‑150 (MQ=255)
                                                   ccaccGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTca                                                                                                    >  2:756179/1‑150 (MQ=255)
                                                   ccaccGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTca                                                                                                    >  2:1423647/1‑150 (MQ=255)
                                                       cGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTcagcag                                                                                                >  1:1106655/1‑150 (MQ=255)
                                                         ccacccacTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCGGCAggg                                                                                              >  2:1707655/1‑149 (MQ=255)
                                                                      cACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGGTCGGCAGGCGCGCGGCCCGCaa                                                                                 >  2:1880022/1‑150 (MQ=255)
                                                                      cACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCaa                                                                                 >  2:1789640/1‑150 (MQ=255)
                                                                      cACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCaa                                                                                 >  2:1416934/1‑150 (MQ=255)
                                                                         tGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGt                                                                              <  1:850351/150‑1 (MQ=255)
                                                                             aCTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGc                                                                          <  1:801843/150‑1 (MQ=255)
                                                                             aCTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGc                                                                          >  2:1982669/1‑150 (MQ=255)
                                                                             aCTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGc                                                                          >  1:1457694/1‑150 (MQ=255)
                                                                              cTCNCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCa                                                                         >  2:1498012/1‑150 (MQ=255)
                                                                                 gCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATa                                                                      <  2:1815176/150‑1 (MQ=255)
                                                                                  ctctTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATAc                                                                     >  2:597554/1‑150 (MQ=255)
                                                                                    ctTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATAccg                                                                   <  2:290941/150‑1 (MQ=255)
                                                                                     tttGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGTAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGccc                                                                                     <  1:1860729/131‑1 (MQ=255)
                                                                                     tttGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATAccgg                                                                  >  2:1685590/1‑150 (MQ=255)
                                                                                          gcagcGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACcggccggc                                                             <  2:478514/150‑1 (MQ=255)
                                                                                              cggcgCCGACCGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCGCGCAAGGTCTGCAATACCGGCCGGCGGGt                                                         <  2:718099/150‑1 (MQ=255)
                                                                                               ggcgCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCa                                                                                           <  2:1929721/115‑1 (MQ=255)
                                                                                               ggcgCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCa                                                                                           >  1:1929721/1‑115 (MQ=255)
                                                                                                    cGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCaa                                                   <  1:1498012/150‑1 (MQ=255)
                                                                                                        aGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACg                                               <  2:1683537/150‑1 (MQ=255)
                                                                                                         gcgcGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGt                                              <  2:1457694/150‑1 (MQ=255)
                                                                                                          cgcgGTGACCAGCACCGGTGCAAGGAAGTTGACGGCCGTGCCCTCTTCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGtt                                             <  1:959372/150‑1 (MQ=255)
                                                                                                            cgGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGa                                           <  2:1614800/150‑1 (MQ=255)
                                                                                                            cgGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGa                                           <  1:1789640/150‑1 (MQ=255)
                                                                                                             ggTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCCATGTCTTCCGGCGGGTGCGCAACACGTTGAg                                          >  1:603413/1‑150 (MQ=255)
                                                                                                              gTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGc                                         <  1:609748/150‑1 (MQ=255)
                                                                                                                 aCCAGCACCGGGGCACGGAAGTTGCCGGACGTGGCGTCGGCCAGTGGCACGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCa                                      <  2:122835/150‑1 (MQ=255)
                                                                                                                 aCCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCa                                      <  2:2203472/150‑1 (MQ=255)
                                                                                                                   cAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGt                                                         <  1:1314311/129‑1 (MQ=255)
                                                                                                                   cAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGt                                                         >  2:1314311/1‑129 (MQ=255)
                                                                                                                   cAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGa                                    >  2:773983/1‑150 (MQ=255)
                                                                                                                   cAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGa                                    >  2:2379863/1‑150 (MQ=255)
                                                                                                                   cAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGGCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGa                                    <  1:1880022/150‑1 (MQ=255)
                                                                                                                           gggCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGTGCaa                                                   >  1:351453/1‑127 (MQ=255)
                                                                                                                           gggCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGTGCaa                                                   <  2:351453/127‑1 (MQ=255)
                                                                                                                           gggCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCa                            <  2:2041960/150‑1 (MQ=255)
                                                                                                                           gggCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCa                            <  2:337795/150‑1 (MQ=255)
                                                                                                                               aaGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGaa                        <  1:243537/150‑1 (MQ=255)
                                                                                                                                aGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTgg                               <  1:698524/142‑1 (MQ=255)
                                                                                                                                aGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTgg                               >  2:698524/1‑142 (MQ=255)
                                                                                                                                  gAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCATCACGTTGAGCCCAGACTTGGGCAGGAAGAt                     <  2:1405362/150‑1 (MQ=255)
                                                                                                                                      ttGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCg                 >  1:1672951/1‑150 (MQ=255)
                                                                                                                                       cactggCCGAGGCCACGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGATCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCTCAGACTTGGGCAg                           <  2:1417266/135‑1 (MQ=255)
                                                                                                                                       cactggCCGAGGCCACGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGATCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCTCAGACTTGGGCAg                           >  1:1417266/5‑139 (MQ=255)
                                                                                                                                        actggCCGAGGCCACGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGc               >  1:1430130/4‑150 (MQ=255)
                                                                                                                                         aCGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGACCAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGcc              <  1:2379863/150‑1 (MQ=255)
                                                                                                                                          cGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCa             <  1:1504700/150‑1 (MQ=255)
                                                                                                                                           ggCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCAt            >  1:359405/1‑150 (MQ=255)
                                                                                                                                            gCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACCGAAATACGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATc           <  1:1331354/150‑1 (MQ=255)
                                                                                                                                            ccactGGCCGAGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATc           >  2:223627/5‑150 (MQ=255)
                                                                                                                                                  gCCGCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATCAGCAAc     <  1:928186/150‑1 (MQ=255)
                                                                                                                                                   ccTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTg                                                                                                           >  1:1502648/1‑47 (MQ=255)
                                                                                                                                                   ccTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTg                                                                                                           <  2:1502648/47‑1 (MQ=255)
                                                                                                                                                   ccTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCgg                >  1:1733753/1‑138 (MQ=255)
                                                                                                                                                   ccTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCgg                <  2:1733753/138‑1 (MQ=255)
                                                                                                                                                     tCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATCAGCAACgtg  <  1:1975894/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACATCGCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATCAGCAACGTG  >  NC_002947/4171707‑4172005

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: