Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 819,747 G→A *151* (TAG→TAA moaE → molybdopterin synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913819,7470GA98.0% 308.6 / ‑1.5 99coding (453/453 nt)moaEmolybdopterin synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (44/53);  minor base . (1/1);  total (45/54)
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.53e-01

CGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGGAGA  >  NC_000913/819614‑819878
                                                                                                                                     |                                                                                                                                     
cGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAAttttt                                                                                                                                  >  1:1625579/1‑139 (MQ=255)
cGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCAGGTAAttttt                                                                                                                                  >  2:953867/1‑139 (MQ=255)
           gggCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGa                                                                                                                       <  2:1671148/139‑1 (MQ=255)
            ggCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGAt                                                                                                                      >  1:75533/1‑139 (MQ=255)
                     ttatGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCg                                                                                                             <  1:843140/139‑1 (MQ=255)
                     ttatGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCg                                                                                                             <  1:506347/139‑1 (MQ=255)
                        tGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATa                                                                                                                     >  2:1210059/1‑128 (MQ=255)
                        tGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATa                                                                                                                     <  1:1210059/128‑1 (MQ=255)
                           aTTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGt                                                                                                       >  1:67327/1‑139 (MQ=255)
                           aTTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGt                                                                                                       >  2:1796605/1‑139 (MQ=255)
                             tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc                                                                                                     >  1:415257/1‑139 (MQ=255)
                             tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc                                                                                                     >  2:174289/1‑139 (MQ=255)
                             tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc                                                                                                     >  1:405780/1‑139 (MQ=255)
                               tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg                                                                                                   >  2:529174/1‑139 (MQ=255)
                               tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg                                                                                                   >  2:1062326/1‑139 (MQ=255)
                                ctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGc                                                                                                  >  2:1786546/1‑139 (MQ=255)
                                ctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGc                                                                                                  >  2:1774066/1‑139 (MQ=255)
                                    aaaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGcc                                                                                                                          <  2:1681648/111‑1 (MQ=255)
                                    aaaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGcc                                                                                                                          >  1:1681648/1‑111 (MQ=255)
                                    aaaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAtt                                                                                              >  1:1452953/1‑139 (MQ=255)
                                    aaaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAtt                                                                                              >  1:888389/1‑139 (MQ=255)
                                    aaaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAtt                                                                                              >  1:170049/1‑139 (MQ=255)
                                         cgcgCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCAcc                                                                                         <  1:329211/139‑1 (MQ=255)
                                         cgcgCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCAcc                                                                                         >  2:7133/1‑139 (MQ=255)
                                          gcgcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCg                                                                                        <  2:1625579/139‑1 (MQ=255)
                                            gcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                      <  1:1172318/139‑1 (MQ=255)
                                            gcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                      <  2:894002/139‑1 (MQ=255)
                                              aCCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCa                                                                                    <  2:888389/139‑1 (MQ=255)
                                                cGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGACAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGt                                                                                  >  1:295833/1‑139 (MQ=255)
                                                 gTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTACGGCAATTTCACCGCTCAGTa                                                                                 >  2:1055084/1‑139 (MQ=255)
                                                  ttCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGc                                                                                                  >  2:1335784/1‑121 (MQ=255)
                                                  ttCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGc                                                                                                  <  1:1335784/121‑1 (MQ=255)
                                                  ttCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGTAATTTCACCGCTCAGTAc                                                                                >  1:120521/1‑139 (MQ=255)
                                                     tGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATa                                                                             <  1:953867/139‑1 (MQ=255)
                                                            gcgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCTGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                            <  1:1036473/85‑1 (MQ=255)
                                                            gcgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                            >  2:1036473/1‑85 (MQ=255)
                                                            gcgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCAcc                                                                                         >  1:170411/1‑120 (MQ=255)
                                                            gcgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCAcc                                                                                         <  2:170411/120‑1 (MQ=255)
                                                             cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                            <  2:92602/84‑1 (MQ=255)
                                                             cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                            >  1:92602/1‑84 (MQ=255)
                                                             cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATATGTTTACg                                                                     >  1:595195/1‑139 (MQ=255)
                                                             cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg                                                                     >  1:116616/1‑139 (MQ=255)
                                                             cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg                                                                     >  1:1657/1‑139 (MQ=255)
                                                             cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg                                                                     >  1:1524468/1‑139 (MQ=255)
                                                             cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg                                                                     >  1:1504687/1‑139 (MQ=255)
                                                                aGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTg                                                                  <  1:1796605/139‑1 (MQ=255)
                                                                    aCGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTc                                                              <  2:415257/139‑1 (MQ=255)
                                                                     cGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCg                                                             <  2:67327/139‑1 (MQ=255)
                                                                          gAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGca                                                      <  2:405780/139‑1 (MQ=255)
                                                                          gAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGca                                                      <  2:116616/139‑1 (MQ=255)
                                                                             ggCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAg                                                   <  2:550447/139‑1 (MQ=255)
                                                                                 aCCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCtc                                               >  1:864302/1‑139 (MQ=255)
                                                                                  ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCtca                                              <  2:75533/139‑1 (MQ=255)
                                                                                    gCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                            <  2:1578876/61‑1 (MQ=255)
                                                                                    gCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                            >  1:1578876/1‑61 (MQ=255)
                                                                                    gCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGca                                                      <  2:399221/129‑1 (MQ=255)
                                                                                    gCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGca                                                      <  2:285592/129‑1 (MQ=255)
                                                                                    gCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGca                                                      >  1:285592/1‑129 (MQ=255)
                                                                                    gCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGca                                                      >  1:399221/1‑129 (MQ=255)
                                                                                       gggTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCAtt                                         <  2:1452953/139‑1 (MQ=255)
                                                                                            gAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCa                                    <  2:1657/139‑1 (MQ=255)
                                                                                            gAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCa                                    <  2:832555/139‑1 (MQ=255)
                                                                                            gAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCa                                    <  2:864302/139‑1 (MQ=255)
                                                                                            gAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCa                                    <  2:1570343/139‑1 (MQ=255)
                                                                                            gAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCa                                    <  2:1524468/139‑1 (MQ=255)
                                                                                            gAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCa                                    <  1:1774066/139‑1 (MQ=255)
                                                                                             aaGCTCGGGAGAGCGATCAGCAGGCGACAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAg                                   <  2:295833/139‑1 (MQ=255)
                                                                                              aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGTCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGt                                  <  1:7133/139‑1 (MQ=255)
                                                                                                cTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGTGCACAGCCTCATCATTTTGCAGt                                  <  1:529174/139‑1 (MQ=255)
                                                                                                  cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGt                                                                          <  2:1650715/97‑1 (MQ=255)
                                                                                                  cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGt                                                                          >  1:1650715/1‑97 (MQ=255)
                                                                                                  cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTg                                                                  >  1:1375956/1‑105 (MQ=255)
                                                                                                  cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTg                                                                  <  2:1375956/105‑1 (MQ=255)
                                                                                                    ggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                      >  1:682943/1‑83 (MQ=255)
                                                                                                    ggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                      <  2:682943/83‑1 (MQ=255)
                                                                                                     gagaGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACTGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTa                           <  1:1062326/139‑1 (MQ=255)
                                                                                                        agCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGt                                  >  2:1650882/1‑129 (MQ=255)
                                                                                                        agCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGt                                  <  1:1650882/129‑1 (MQ=255)
                                                                                                          ccaTCAGCAGGCGGCAAAACGGTGGTAATTTTTTGTTAGGCGGATAAGGCACCGGTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCAATTTGCAGTATCCTTAAGata                      <  1:1055084/137‑1 (MQ=255)
                                                                                                          cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGTAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGata                      <  2:120521/139‑1 (MQ=255)
                                                                                                          cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGata                      <  1:1786546/139‑1 (MQ=255)
                                                                                                            aTCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGTTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATAtt                    >  2:586638/1‑139 (MQ=255)
                                                                                                               agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGca                                                      >  2:809492/1‑102 (MQ=255)
                                                                                                               agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGca                                                      <  1:809492/102‑1 (MQ=255)
                                                                                                               agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGt                                  >  2:915984/1‑122 (MQ=255)
                                                                                                               agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGt                                  <  1:915984/122‑1 (MQ=255)
                                                                                                                gcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCtt                <  2:170049/139‑1 (MQ=255)
                                                                                                                  agGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTtat              <  1:174289/139‑1 (MQ=255)
                                                                                                                  agGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTtat              <  2:1504687/139‑1 (MQ=255)
                                                                                                                    gcggcAAAACGCTGGTAATTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATc           >  1:1662858/1‑139 (MQ=255)
                                                                                                                      ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                      <  1:509676/65‑1 (MQ=255)
                                                                                                                      ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                      >  2:509676/1‑65 (MQ=255)
                                                                                                                      ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGt                                  <  1:514857/115‑1 (MQ=255)
                                                                                                                      ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGt                                  >  2:514857/1‑115 (MQ=255)
                                                                                                                             cGCTGGTAATTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGgaga  <  2:1662858/139‑1 (MQ=255)
                                                                                                                              gCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCTTTAAGATATTCCTTATATCTTCAGgaga  <  2:1113730/139‑1 (MQ=255)
                                                                                                                              gCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGgaga  <  1:582082/139‑1 (MQ=255)
                                                                                                                              gCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGgaga  <  2:1537795/139‑1 (MQ=255)
                                                                                                                                  gTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGt                                  >  1:1769375/1‑103 (MQ=255)
                                                                                                                                     |                                                                                                                                     
CGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGA‑‑TGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGGAGA  >  NC_000913/819614‑819878

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: