Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,207,576 A→G 5.6% intergenic (‑44/‑571) purH ← / → rrsE IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase/16S ribosomal RNA of rrnE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,207,5760AG 5.6% 373.3 / 9.2 124intergenic (‑44/‑571)purH/rrsEIMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase/16S ribosomal RNA of rrnE operon
Reads supporting (aligned to +/- strand):  ref base A (66/51);  new base G (2/5);  total (68/56)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 124 total aligned reads displayed.

CAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATAC  >  NC_000913/4207428‑4207723
                                                                                                                                                    |                                                                                                                                                   
cagcagCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAAt                                                                                                                                                    <  2:156366/150‑1 (MQ=255)
  gcagcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAAttt                                                                                                                                                  >  1:813404/1‑150 (MQ=255)
   cagcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGAGTTGACTATTACAGAGAGCGTTAGCTGAAtttt                                                                                                                                                 >  1:629433/1‑150 (MQ=255)
      cTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTcg                                                                                                                                              <  2:1794122/150‑1 (MQ=255)
       tCCACACCGCGGGCGGAAAGTGCCTGGGAGAATTCGACGATTCCGGCTTTGTCTGAAACACTGAGCAGTGCGCGGCGGACTGGACGACGTGGTGGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTcgc                                                                                                                                             <  1:626495/150‑1 (MQ=255)
        ccACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGGTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTcgcg                                                                                                                                            >  2:250548/1‑150 (MQ=255)
         cacaCCGCGTGCGGAAAGTGCCTGGGCGAATTCGATGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGa                                                                                                                                           >  1:986038/1‑150 (MQ=255)
         cacaCCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGa                                                                                                                                           >  1:1013054/1‑150 (MQ=255)
              cgcgTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAAc                                                                                                                                      <  2:1999701/150‑1 (MQ=255)
                cgTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTc                                                                                                                                    <  2:1795131/150‑1 (MQ=255)
                  tGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGCGCGTTAGCTGAATTTTTCGCGAAAAACTCAg                                                                                                                                  >  1:390047/1‑150 (MQ=255)
                   gCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGc                                                                                                                                 <  2:586907/150‑1 (MQ=255)
                    cGGAAAGTGCCGGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGCGCGCGGCGGACTGGACGACGTTGTTGCATTGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAATTCAGCt                                                                                                                                <  1:1308191/150‑1 (MQ=255)
                      gAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTaa                                                                                                                              >  1:512671/1‑150 (MQ=255)
                         aGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAATTCAGCTAACGc                                                                                                                           <  1:861926/150‑1 (MQ=255)
                          gTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGcc                                                                                                                          <  2:1405661/150‑1 (MQ=255)
                           ggCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAATTCAGCTAACGctc                                                                                                                         <  1:1853965/149‑3 (MQ=255)
                            gCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTTCAGAGAGCGTAAGCTGAATTTTTCGCGAAAAACTCAGCTAACGcccc                                                                                                                        <  2:932429/150‑1 (MQ=255)
                            gCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGcccc                                                                                                                        <  1:1073516/150‑1 (MQ=255)
                            gCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATCTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGcccc                                                                                                                        <  1:153665/150‑1 (MQ=255)
                             ccTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAATTCAGCTAACGctctc                                                                                                                       <  2:734204/150‑5 (MQ=255)
                              cTGGGCGAATTCGACGATCCCGGCTTTGTCAGCAACACTGCGCAGAGGGCGGCGGACGGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTa                                                                                                                      <  2:1632503/150‑1 (MQ=255)
                              cTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTa                                                                                                                      >  2:853231/1‑150 (MQ=255)
                                gggCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAAc                                                                                                                    >  2:494675/1‑150 (MQ=255)
                                   ggAATTCGAAGATACCGGCTTTGTCAGAAAAACTGAGCAGAGCTCGGCGGCCTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTTTTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTAATCTAACGCCCCTAACggg                                                                                                                 <  2:412643/149‑1 (MQ=255)
                                    gAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACgggg                                                                                                                <  2:749224/150‑1 (MQ=255)
                                    gAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACgggg                                                                                                                <  2:796608/150‑1 (MQ=255)
                                      aTTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCa                                                                                                              >  1:650433/1‑150 (MQ=255)
                                          gacgaTACCGGCGTTGTCAGAAACCCTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCt                                                                                                          <  1:669194/150‑1 (MQ=255)
                                            cgaTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTa                                                                                                        >  1:816818/1‑150 (MQ=255)
                                              aTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAtt                                                                                                      >  2:907572/1‑150 (MQ=255)
                                                 cctgCTTTGTCAGAAACACTGAGCAGAGCGCGGCTGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACAATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAttttt                                                                                                   <  1:1315224/147‑1 (MQ=255)
                                                 ccGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAttttt                                                                                                   >  2:913619/1‑150 (MQ=255)
                                                  cGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTc                                                                                                  >  2:789913/1‑150 (MQ=255)
                                                   ggCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCg                                                                                                 >  1:1268484/1‑150 (MQ=255)
                                                    gCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGc                                                                                                >  1:1916433/1‑150 (MQ=255)
                                                    gCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTCGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGc                                                                                                <  1:762645/150‑1 (MQ=255)
                                                       ttGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCg                                                                                             <  2:1691858/150‑1 (MQ=255)
                                                       ttGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCg                                                                                             >  1:131407/1‑150 (MQ=255)
                                                         gTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTACCGGGGCATCCTTATTTTTCGCCCGCa                                                                                           >  1:1015682/1‑150 (MQ=255)
                                                         gTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCa                                                                                           >  1:434993/1‑150 (MQ=255)
                                                          tCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTAGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATACATATTTTTCGACCGCAt                                                                                          >  1:1802760/1‑150 (MQ=255)
                                                           cAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCGTTATTTTTCGCCCGCAtt                                                                                         >  2:396760/1‑150 (MQ=255)
                                                            aGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTg                                                                                        >  2:1570596/1‑150 (MQ=255)
                                                            aGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTg                                                                                        <  2:1447560/150‑1 (MQ=255)
                                                                 cacTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCATAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACg                                                                                   >  1:487049/1‑150 (MQ=255)
                                                                  acTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGa                                                                                  >  1:598278/1‑150 (MQ=255)
                                                                  acTGAGCAGAGCGCGGCCGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGACCGCAAAGTAACGa                                                                                  >  1:1108312/1‑150 (MQ=255)
                                                                      aGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAAc                                                                              >  2:510136/1‑150 (MQ=255)
                                                                       gCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACg                                                                             >  1:358106/1‑150 (MQ=255)
                                                                         agagCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGtt                                                                           <  1:1182666/150‑1 (MQ=255)
                                                                           agCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACCAAAACGTTTg                                                                         <  2:1177179/150‑1 (MQ=255)
                                                                           agCGAGTCGGAATGTACGACGTTGTTGCATTGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTTTTTTTCGCCCGCATTGTAACGAAAAAGTTTg                                                                         <  2:1926261/150‑1 (MQ=255)
                                                                             cgcgGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTgcg                                                                       <  2:1682121/150‑1 (MQ=255)
                                                                                    gACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCt                                                                <  2:512671/150‑1 (MQ=255)
                                                                                              cgTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGCGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCACGCATTGTAACGAAAACGTTTGCGCACCGCGCGCGAAtttt                                                      >  2:143776/1‑150 (MQ=255)
                                                                                              cgTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAAtttt                                                      >  1:1866235/1‑150 (MQ=255)
                                                                                                ttgttgCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTCAATTTTTCGCGAAAAACTCAGCTAACGCCCCAAAAGGGGCATCCTTATTTTTCGCCCGCATTGTAAAGAAAACGTTTGCTCAACGCTCGCGAAttttta                                                    >  2:1924132/1‑149 (MQ=255)
                                                                                                    tgCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCtt                                                >  1:421827/1‑150 (MQ=255)
                                                                                                     gCATGGTAACACCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGCGGCATCCATATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCttt                                               <  1:1244639/150‑1 (MQ=255)
                                                                                                      cATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTc                                              >  2:947548/1‑150 (MQ=255)
                                                                                                      cATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTc                                              >  2:1403202/1‑150 (MQ=255)
                                                                                                       aTGGTAAATCCCCTGGATTTGACTATTTCAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGATTTTTTCTCGTTCa                                             >  1:1673476/1‑150 (MQ=255)
                                                                                                       aTGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCa                                             >  2:473532/1‑150 (MQ=255)
                                                                                                       aTGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCa                                             <  1:907044/150‑1 (MQ=255)
                                                                                                       aTGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCa                                             <  2:1013054/150‑1 (MQ=255)
                                                                                                        tGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCaa                                            >  2:1303086/1‑150 (MQ=255)
                                                                                                          gTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATg                                          <  1:310204/150‑1 (MQ=255)
                                                                                                          gTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATg                                          <  1:454897/150‑1 (MQ=255)
                                                                                                              aTCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGa                                      >  1:751686/1‑150 (MQ=255)
                                                                                                              aTCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTATTCGCGAAAAACTCAGCTAACGCCCCAAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAACCGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGa                                      <  2:314386/150‑1 (MQ=255)
                                                                                                                ccccTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATc                                    <  1:1721623/150‑1 (MQ=255)
                                                                                                                 cccTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAACTAAGCTAACGCCCCTAACGGGGCATCATTATTTTTCGCCCGCATTGTAACGAAAACGATTGCGCAACGCTCGCGAATTTTTCTCGTTCAATGGTGAGca                                   >  2:684969/1‑150 (MQ=255)
                                                                                                                 cccTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATca                                   >  1:1100991/1‑150 (MQ=255)
                                                                                                                  ccTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTTTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATcac                                  >  1:1499026/1‑150 (MQ=255)
                                                                                                                  ccTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATcaa                                  >  1:659735/1‑149 (MQ=255)
                                                                                                                    tGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACaa                                <  1:1376988/150‑1 (MQ=255)
                                                                                                                     ggATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAAt                               <  2:47902/150‑1 (MQ=255)
                                                                                                                      gATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAAtt                              >  2:831011/1‑150 (MQ=255)
                                                                                                                      gATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAAtt                              <  2:813404/150‑1 (MQ=255)
                                                                                                                       aTTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCGTTCAATGGTGATCACAAttt                             >  2:1818022/1‑150 (MQ=255)
                                                                                                                         ttGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTg                           >  2:1096545/1‑150 (MQ=255)
                                                                                                                          tGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGa                          >  1:1190906/1‑150 (MQ=255)
                                                                                                                           gACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGAc                         >  1:908079/1‑150 (MQ=255)
                                                                                                                           gACAATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGAc                         <  2:1065383/150‑1 (MQ=255)
                                                                                                                            aCTATTACAGAGAGCGTAAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCATCGCTCGCGAATTTTTCTCTTTCAATGGTGATCACTATTTTGACt                        <  2:1233240/150‑1 (MQ=255)
                                                                                                                             cTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGTCtg                       >  2:710951/1‑150 (MQ=255)
                                                                                                                             cTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACtg                       >  1:713443/1‑150 (MQ=255)
                                                                                                                              tATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTATCAATCGTGATCACAATTTTGACtgt                      >  2:649148/1‑150 (MQ=255)
                                                                                                                                ttACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTgg                    <  2:1268484/150‑1 (MQ=255)
                                                                                                                                 tACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCAAACGGGGCATCCTTAATTTTCGCCCGCATTGTAACGAAAACGTTAGCGCAACGCTCGCCAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGt                   >  2:1783440/1‑150 (MQ=255)
                                                                                                                                  aCAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGtt                  >  2:1039301/1‑150 (MQ=255)
                                                                                                                                        ggCGTTAGCTGAGTTTTTCGCGAAAAATTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTg            >  2:1636214/2‑150 (MQ=255)
                                                                                                                                          cGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGTGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATTGTGATCACAATTTTGACTGTGGTTACCGTggg          <  1:474443/150‑1 (MQ=255)
                                                                                                                                          cGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTggg          <  2:830369/150‑1 (MQ=255)
                                                                                                                                          cGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTggg          <  2:1214291/150‑1 (MQ=255)
                                                                                                                                          cGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTggg          >  1:636681/1‑150 (MQ=255)
                                                                                                                                           gTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCCACGCTCGCGAATTTTTCTCTTTCAATGGGGATCACAATTTTGACTGTGGTTACCGGGGGc         >  2:377344/1‑150 (MQ=255)
                                                                                                                                            ttAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCa        <  1:907572/150‑1 (MQ=255)
                                                                                                                                                  gAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATAc  >  2:1229310/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATAC  >  NC_000913/4207428‑4207723

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