Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,497 T→G 6.7% intergenic (+16/‑116) aroH → / → ydiE 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible/hemin uptake protein HemP homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,789,4970TG 6.7% 271.8 / 8.9 104intergenic (+16/‑116)aroH/ydiE3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible/hemin uptake protein HemP homolog
Reads supporting (aligned to +/- strand):  ref base T (45/52);  new base G (4/3);  total (49/55)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

Only 100 of 105 total aligned reads displayed.

GCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCT  >  NC_000913/1789348‑1789645
                                                                                                                                                     |                                                                                                                                                    
gcgAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCt                                                                                                                                                      >  1:839778/1‑150 (MQ=255)
   aaGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTTCTGTCAATCCATTACCGAACCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGAc                                                                                                                                                   <  1:1587842/150‑1 (MQ=255)
    aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGAcc                                                                                                                                                  >  2:1353509/1‑149 (MQ=255)
     ggAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACgg                                                                                                                                                 <  2:1326212/150‑1 (MQ=255)
     ggAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACgg                                                                                                                                                 >  1:1619499/1‑150 (MQ=255)
        aCGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCGTCTGAATGCGTGCCCATTCGTGACGGTAt                                                                                                                                              >  2:177556/1‑150 (MQ=255)
          gCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATgg                                                                                                                                            >  2:1468219/1‑150 (MQ=255)
           cAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATggg                                                                                                                                           >  1:530904/1‑150 (MQ=255)
            aaaaacTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCTCTTCTGAATGCGTGCCCATTCCGTCAGGAATGGGc                                                                                                                                          <  2:934653/150‑1 (MQ=255)
            aaaaaaTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGc                                                                                                                                          <  1:560610/150‑1 (MQ=255)
             aaaaaTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCGGGGAGGATCCCGAACGCCTGGCCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCGTCAGGAATGGGCa                                                                                                                                         <  2:122244/150‑1 (MQ=255)
               aaaTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCGTCAGGAATGGGCacg                                                                                                                                       <  2:1062033/150‑3 (MQ=255)
                 aTCGTCGGCAGTCAGCCGCTCACTTACTGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTc                                                                                                                                     <  2:1689394/150‑1 (MQ=255)
                   cgtcgGCAGTCAGCCTCTCACATACGGTCAATACATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCGGGTCGAAATACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTg                                                                                                                                   <  1:1722295/150‑1 (MQ=255)
                   cgtcgGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGCCCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTg                                                                                                                                   >  2:1581377/1‑150 (MQ=255)
                   cgtcgGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGCCGGAATGGGCATTTCTg                                                                                                                                   >  2:73234/1‑150 (MQ=255)
                      cgGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTgcgc                                                                                                                                >  1:1199032/1‑150 (MQ=255)
                       ggCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCa                                                                                                                               <  2:621042/150‑1 (MQ=255)
                        gCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCaa                                                                                                                              >  1:1531250/1‑150 (MQ=255)
                           gTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACtt                                                                                                                           >  2:1014664/1‑150 (MQ=255)
                           gTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACtt                                                                                                                           >  1:1283447/1‑150 (MQ=255)
                             cAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACttgt                                                                                                                         >  2:631466/1‑150 (MQ=255)
                               gCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACttgttt                                                                                                                       >  2:1516435/1‑149 (MQ=255)
                               gCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACttgttg                                                                                                                       <  1:1004499/150‑1 (MQ=255)
                               gCCGCTCACTTACGGGCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACttgttg                                                                                                                       <  2:759449/150‑1 (MQ=255)
                                    gcACTTACGGTCAATCCATAACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtcttc                                                                                                                  <  2:1499852/149‑1 (MQ=255)
                                      aCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTtctc                                                                                                                >  1:366510/1‑150 (MQ=255)
                                       cTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGACCGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTGCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTATGCGCAACTTGGTGTCTTCTca                                                                                                               >  1:413559/1‑150 (MQ=255)
                                          aCGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTcaaca                                                                                                            <  1:1153119/150‑1 (MQ=255)
                                           cGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTcaacaa                                                                                                           >  2:1061684/1‑150 (MQ=255)
                                             gTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAACACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAAt                                                                                                         >  2:1719854/1‑150 (MQ=255)
                                              tCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAAtt                                                                                                        >  2:679708/1‑150 (MQ=255)
                                              tCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAAtt                                                                                                        >  1:846224/1‑150 (MQ=255)
                                               cAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTa                                                                                                       <  1:1675999/150‑1 (MQ=255)
                                               cAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAAGGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTa                                                                                                       <  1:1175207/150‑1 (MQ=255)
                                                 aTCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGATTTCTTCTCAACAAATTACt                                                                                                     >  1:778340/1‑150 (MQ=255)
                                                 aTCCATTACCGACCCGTGTCTGGGCTGGGAGGATACAGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATCCCCGCTGCTGAATGCGTGCCCATTCATGACGGAATGGGCATTTCTGCGCAACTTGTTTGCTTCTCAACAAAATACt                                                                                                     >  1:158130/1‑150 (MQ=255)
                                                    cATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACtgct                                                                                                  >  2:133954/1‑150 (MQ=255)
                                                    cATTACCGACCCGTGGCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACtgct                                                                                                  <  2:369512/150‑1 (MQ=255)
                                                       tACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCCACTTGTTGTCTTCTCAACAAATTACtgcttgc                                                                                               >  1:930174/1‑150 (MQ=255)
                                                        aCCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTAGTCTTCTCAACAAATTACtgcttgct                                                                                              >  2:544874/1‑150 (MQ=255)
                                                        aCCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTACCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACtgcttgct                                                                                              >  2:1669048/1‑150 (MQ=255)
                                                            aCCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTgg                                                                                          <  2:607810/150‑1 (MQ=255)
                                                            aCCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTgg                                                                                          <  2:372118/150‑1 (MQ=255)
                                                             cccGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGt                                                                                         >  2:1669355/1‑150 (MQ=255)
                                                                  gtCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGcc                                                                                    <  2:773699/150‑1 (MQ=255)
                                                                  gtCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGcc                                                                                    <  2:812658/150‑1 (MQ=255)
                                                                     tgggctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCata                                                                                 <  2:94260/150‑1 (MQ=255)
                                                                      gggctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCataa                                                                                <  1:1822582/150‑1 (MQ=255)
                                                                         gtgggAGTATACCTAACGCATGTTCGAAAAACTCGCGTCTCCTGTAGATACCGGCTTCTGAATGCGTGCCCATTCCTGACGGATTGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAatat                                                                             <  1:1398958/149‑1 (MQ=255)
                                                                          tgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATAtt                                                                            <  2:804086/150‑1 (MQ=255)
                                                                            ggaggaTACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATTTTGa                                                                          <  1:150416/150‑1 (MQ=255)
                                                                               ggaTACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATTATATTGataa                                                                       >  1:1115724/1‑150 (MQ=255)
                                                                                  tACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGAtaataa                                                                    <  1:1156816/150‑1 (MQ=255)
                                                                                    ccGAACGCCAGGTCGAAAAACTCGCCTCTGAGGTAGATACCCGCTACTGAATGCGAGCCCATTCCTGGCGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGa                                                                  >  1:956783/1‑150 (MQ=255)
                                                                                     cGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGaa                                                                 <  2:839778/150‑1 (MQ=255)
                                                                                      gAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACACGCTTCTGAATGCGTGACCATTCCTGACTGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATTATATTGATAATAAGAAt                                                                <  1:1778878/150‑1 (MQ=255)
                                                                                          gCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCAtt                                                            <  1:934742/150‑1 (MQ=255)
                                                                                            cTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAAGCATTGt                                                          <  1:1070740/150‑1 (MQ=255)
                                                                                             tgttCGAAAAACTCGCCTATGCGGTTGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCTCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGGCAGCCATAATATTGATAATAAGATTCATTGtt                                                         <  2:409548/147‑1 (MQ=255)
                                                                                              ggTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTta                                                        <  1:1418722/150‑1 (MQ=255)
                                                                                               gTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTtat                                                       <  2:1561824/150‑1 (MQ=255)
                                                                                                tCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCTTAATATTGATAATAAGAATCATTGTtata                                                      >  1:345477/1‑150 (MQ=255)
                                                                                                  gAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATc                                                    >  1:1523323/1‑150 (MQ=255)
                                                                                                   aaaaaCTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCa                                                   <  1:1468219/150‑1 (MQ=255)
                                                                                                        cTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAattat                                              <  2:1189257/150‑1 (MQ=255)
                                                                                                        cTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAattat                                              >  1:6508/1‑150 (MQ=255)
                                                                                                        cTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAattat                                              >  2:1470831/1‑150 (MQ=255)
                                                                                                           gCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTTACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAattattat                                           <  1:1353509/150‑1 (MQ=255)
                                                                                                             ctctGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAAttattatta                                         <  2:1574168/150‑1 (MQ=255)
                                                                                                                  cGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAAtttt                                    <  2:530904/150‑1 (MQ=255)
                                                                                                                  cGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAAtttt                                    <  2:1283447/150‑1 (MQ=255)
                                                                                                                    gTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtcttc                                                                                                                  <  1:1159811/70‑1 (MQ=255)
                                                                                                                    gTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtcttc                                                                                                                  >  2:1159811/1‑70 (MQ=255)
                                                                                                                      aGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGTAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATg                                <  1:1330530/150‑1 (MQ=255)
                                                                                                                       gATACCCGCTTCTGAATTCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGc                               <  2:1352081/150‑1 (MQ=255)
                                                                                                                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGt                             >  2:536328/1‑150 (MQ=255)
                                                                                                                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGt                             >  2:514307/1‑150 (MQ=255)
                                                                                                                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCAGGTCAGCCATAAAATTGATAATAAGAATCATTGTTATATAAAGTATTATTAATTTTTATGCGt                             >  1:925049/1‑150 (MQ=255)
                                                                                                                          aCCAGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGtt                            <  1:1736388/150‑1 (MQ=255)
                                                                                                                              gCTTCTGAATGCGTGCCCATTCCTGACTGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTTCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATAc                        <  1:902638/150‑1 (MQ=255)
                                                                                                                                   tGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATa                   >  1:1008995/1‑150 (MQ=255)
                                                                                                                                   tGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATa                   <  2:1531250/150‑1 (MQ=255)
                                                                                                                                   tGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATa                   <  2:1551272/150‑1 (MQ=255)
                                                                                                                                       tGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAg               >  2:1064226/1‑150 (MQ=255)
                                                                                                                                        gCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGa              <  1:365729/150‑1 (MQ=255)
                                                                                                                                         cGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACATATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGaa             >  1:1184005/1‑150 (MQ=255)
                                                                                                                                         aatGCCCATTCCGTCAGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGaa             >  1:1246121/3‑150 (MQ=255)
                                                                                                                                            gCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc          <  2:1523323/150‑1 (MQ=255)
                                                                                                                                            gCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc          >  1:874281/1‑150 (MQ=255)
                                                                                                                                            gCCCATTCCGTCAGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc          >  1:1155162/1‑150 (MQ=255)
                                                                                                                                            gCCCATTCCGTCAGGAATGGGCATTTCTGCGCAACTTGATGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAAAATTGATAATAAGAATCATAGTTATATCAATTATAATTAATTTTTATGCGTTATACGGATAGCAGAAAAc          >  1:713294/1‑150 (MQ=255)
                                                                                                                                            gCCCATTCCGTCAGGAATGGGCATTTCTGCGAAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTTTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc          >  1:905387/1‑150 (MQ=255)
                                                                                                                                              ccATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTc        <  2:1011370/150‑1 (MQ=255)
                                                                                                                                                aTTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCAc      <  1:1669048/150‑1 (MQ=255)
                                                                                                                                                  tCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGc    >  1:333822/1‑150 (MQ=255)
                                                                                                                                                  tCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGc    >  1:1288559/1‑150 (MQ=255)
                                                                                                                                                   ccTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGcc   <  2:930174/150‑1 (MQ=255)
                                                                                                                                                    cTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCt  >  1:1395509/1‑150 (MQ=255)
                                                                                                                                                    cTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCt  >  1:1544807/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCT  >  NC_000913/1789348‑1789645

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