Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,500 C→A 6.9% intergenic (+19/‑113) 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,5000CA 6.9% 279.8 / 9.5 102intergenic (+19/‑113)aroH/ydiE3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible/hemin uptake protein HemP homolog
Reads supporting (aligned to +/- strand):  ref base C (44/51);  new base A (4/3);  total (48/54)
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.11e-01

Only 100 of 105 total aligned reads displayed.

AAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAA  >  NC_000913/1789351‑1789649
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                               >  1:1283447/1‑150 (MQ=255)
                        gTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACtt                                                                                                                               >  2:1014664/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                                                                                                            >  1:846224/1‑150 (MQ=255)
                                           tCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAAtt                                                                                                            >  2:679708/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)
                                                    tACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCCACTTGTTGTCTTCTCAACAAATTACtgcttgc                                                                                                   >  1:930174/1‑150 (MQ=255)
                                                     aCCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTAGTCTTCTCAACAAATTACtgcttgct                                                                                                  >  2:544874/1‑150 (MQ=255)
                                                     aCCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTACCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACtgcttgct                                                                                                  >  2:1669048/1‑150 (MQ=255)
                                                         aCCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTgg                                                                                              <  2:372118/150‑1 (MQ=255)
                                                         aCCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTgg                                                                                              <  2:607810/150‑1 (MQ=255)
                                                          cccGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGt                                                                                             >  2:1669355/1‑150 (MQ=255)
                                                               gtCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGcc                                                                                        <  2:812658/150‑1 (MQ=255)
                                                               gtCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGcc                                                                                        <  2:773699/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)
                                                                      ctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAatat                                                                                 >  2:1798465/1‑150 (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)
                                                                                aCCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAg                                                                       <  2:1232105/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                                                  >  1:6508/1‑150 (MQ=255)
                                                                                                     cTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAattat                                                  >  2:1470831/1‑150 (MQ=255)
                                                                                                     cTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAattat                                                  <  2:1189257/150‑1 (MQ=255)
                                                                                                        gCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTTACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAattattat                                               <  1:1353509/150‑1 (MQ=255)
                                                                                                          ctctGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAAttattatta                                             <  2:1574168/150‑1 (MQ=255)
                                                                                                               cGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAAtttt                                        <  2:530904/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:1551272/150‑1 (MQ=255)
                                                                                                                                tGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATa                       <  2:1531250/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:1288559/1‑150 (MQ=255)
                                                                                                                                               tCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGc        >  1:333822/1‑150 (MQ=255)
                                                                                                                                                ccTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGcc       <  2:930174/150‑1 (MQ=255)
                                                                                                                                                 cTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCt      >  1:1544807/1‑150 (MQ=255)
                                                                                                                                                 cTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCt      >  1:1395509/1‑150 (MQ=255)
                                                                                                                                                     cGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGTATAGCAGAAAACTCACGCCTGaaa  >  2:604279/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAA  >  NC_000913/1789351‑1789649

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: