Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,497 T→G 5.6% 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 5.6% 286.4 / 5.9 101intergenic (+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 (47/48);  new base G (3/3);  total (50/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

Only 100 of 103 total aligned reads displayed.

GCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTG  >  NC_000913/1789348‑1789646
                                                                                                                                                     |                                                                                                                                                     
gcgAAGGAACTCAAAAAATCGTCGGCAGTCATCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGCTACCCGCTTCTGAATGCGTGCCCATTCCt                                                                                                                                                       >  2:1557602/1‑150 (MQ=255)
  gAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGa                                                                                                                                                     >  1:1360035/1‑150 (MQ=255)
         cGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATg                                                                                                                                              >  1:17972/1‑150 (MQ=255)
           cAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGGCAATCCATTACCGACCCGTGGCTGGGCTGGGCGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATggg                                                                                                                                            <  1:1362182/150‑1 (MQ=255)
            aaaaaaTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGc                                                                                                                                           <  1:428932/150‑1 (MQ=255)
            aaaaaaTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCGGCAGGAATGGGc                                                                                                                                           <  2:1431779/150‑1 (MQ=255)
            aaaaaaTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCGGACGGAATGGGc                                                                                                                                           >  2:82050/1‑150 (MQ=255)
             aaaaaTCGTCGGCAGGCAGCCGCTCACCTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCGTCAGGAATGGGCa                                                                                                                                          <  1:610274/150‑1 (MQ=255)
              aaaaTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCAt                                                                                                                                         >  1:1168295/1‑150 (MQ=255)
              aaaaTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCGTCAGGAATGGGCac                                                                                                                                         <  2:514758/150‑2 (MQ=255)
               aaaTCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGAAACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCaat                                                                                                                                        >  1:1610049/1‑148 (MQ=255)
                    gtcgGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTgc                                                                                                                                   <  1:616946/150‑1 (MQ=255)
                    gtcgGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTgc                                                                                                                                   <  2:1047061/150‑1 (MQ=255)
                     tcgGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTgcg                                                                                                                                  <  1:693910/150‑1 (MQ=255)
                       ggCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCa                                                                                                                                >  2:1332639/1‑150 (MQ=255)
                        gCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGCAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTCACGGAATGGGCATTTCTGCGCaa                                                                                                                               >  1:57271/1‑150 (MQ=255)
                          aGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACt                                                                                                                             >  1:624555/1‑150 (MQ=255)
                           gTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACtt                                                                                                                            >  1:1360082/1‑150 (MQ=255)
                               gCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACttgttg                                                                                                                        <  2:1168295/150‑1 (MQ=255)
                                  gCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGACTCTGCGGTAGATACCCTCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTTTTGtct                                                                                                                     >  1:127291/1‑150 (MQ=255)
                                     cACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGTATTTCTGCGCAACTTGTTGtcttct                                                                                                                  >  1:1717044/1‑150 (MQ=255)
                                            ggTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACaaa                                                                                                           <  1:188410/150‑1 (MQ=255)
                                              tCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAAtt                                                                                                         <  1:1395624/150‑1 (MQ=255)
                                                aaTCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTAc                                                                                                       <  1:1478997/150‑1 (MQ=255)
                                                 aTCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACt                                                                                                      >  2:1625790/1‑150 (MQ=255)
                                                 aTCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACt                                                                                                      >  2:207623/1‑150 (MQ=255)
                                                  tCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACtg                                                                                                     >  2:3034/1‑150 (MQ=255)
                                                   ccATTACCGTCCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTATGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACtgc                                                                                                    <  2:1597156/150‑1 (MQ=255)
                                                   ccATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACtgc                                                                                                    <  1:1561754/150‑1 (MQ=255)
                                                       tACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACtgcttgc                                                                                                >  2:1315516/1‑150 (MQ=255)
                                                        aCCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACACGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACtgcttgct                                                                                               <  2:924861/150‑1 (MQ=255)
                                                          cGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAACCTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGCCGGAATGGGCATTTCTGCGCCACTTGTTGTCTTCTCAACAAATTACTGCTTGctct                                                                                             >  1:1595804/1‑150 (MQ=255)
                                                            aCCCTTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTgg                                                                                           <  1:143020/150‑1 (MQ=255)
                                                            aCCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTgg                                                                                           <  1:1612603/150‑1 (MQ=255)
                                                             cccGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGt                                                                                          <  1:999040/150‑1 (MQ=255)
                                                             cccGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGt                                                                                          >  1:838599/1‑150 (MQ=255)
                                                             cccGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGt                                                                                          >  1:620146/1‑150 (MQ=255)
                                                             cccGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGt                                                                                          <  2:17972/150‑1 (MQ=255)
                                                               cGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCa                                                                                        >  1:1339885/1‑150 (MQ=255)
                                                                  gtCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGCTACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTTCTCTGCTCAGcc                                                                                     >  2:1344010/1‑150 (MQ=255)
                                                                  gtCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGcc                                                                                     >  1:1594544/1‑150 (MQ=255)
                                                                     tgggctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCata                                                                                  >  2:1645825/1‑150 (MQ=255)
                                                                       ggctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCataat                                                                                <  1:1353628/150‑1 (MQ=255)
                                                                        gctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCataata                                                                               <  1:777811/150‑1 (MQ=255)
                                                                        gcggggAGGATATCGAACGCATGGGCGAAAAACTCGCCGCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGGATTTCTGCGCAAGTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCataata                                                                               <  2:433708/147‑1 (MQ=255)
                                                                         ctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAatat                                                                              >  1:52276/1‑150 (MQ=255)
                                                                         ctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTGGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAatat                                                                              >  2:983448/1‑150 (MQ=255)
                                                                          tgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATAtt                                                                             <  1:803183/150‑1 (MQ=255)
                                                                          tgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATAtt                                                                             <  1:983766/150‑1 (MQ=255)
                                                                           gggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTg                                                                            <  2:1188706/150‑1 (MQ=255)
                                                                             gaggaTACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGACCATTCCTGACGGAATGGGCAATTCTGCGCAACTTGAGGTCTTCTCAACAAATTACTGCTTGCTCTGTGCAGCCATAATATTGat                                                                          >  1:334408/1‑150 (MQ=255)
                                                                               ggaTACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTACTGCGCAACTTCTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGataa                                                                        >  1:701862/1‑150 (MQ=255)
                                                                                  tACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGCCGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGAtaataa                                                                     >  2:300031/1‑150 (MQ=255)
                                                                                    ccGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGa                                                                   >  2:1189067/1‑150 (MQ=255)
                                                                                      gAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAAt                                                                 >  2:1141562/1‑150 (MQ=255)
                                                                                         cGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGATGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCAt                                                              <  2:1190941/150‑1 (MQ=255)
                                                                                          gCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCAtt                                                             >  1:884449/1‑150 (MQ=255)
                                                                                           ccTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTg                                                            >  2:1395641/1‑150 (MQ=255)
                                                                                            cTGGTCGAAAACCTGGCGTATGCGGTAGATACCCGGTTTTGTATGGGTGCCCTTTCCTGACGGATTGGGCATTTCTGCGCAACTTGTTGTCTTCTCATCAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGt                                                           <  1:1557602/150‑1 (MQ=255)
                                                                                            cTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGt                                                           >  1:942554/1‑150 (MQ=255)
                                                                                               gTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTtat                                                        <  2:842691/150‑1 (MQ=255)
                                                                                                 cGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGCCGGATTGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATATGATTCATTGTTatat                                                      >  1:1199496/1‑150 (MQ=255)
                                                                                                 cGAAAAACGCGCCTCTGCGGGAGATACCCGCTTCTGAATGCGTGCCCATTCGTGACGGAATGGGCATTTCTGGGCAACTTGTTGTCTTCTCAAGAAATTACTGCTTGGTCTGGTCAGCCATAATATTGATAAGAAGAATCTTTGTTatat                                                      >  2:527588/1‑150 (MQ=255)
                                                                                                  gAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGTTAATAAGAATCATTGTTatatg                                                     >  1:97342/1‑149 (MQ=255)
                                                                                                   aaaaaCTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAATAATCATTGTTATATCa                                                    <  2:640717/150‑1 (MQ=255)
                                                                                                        cTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCAGAATATTGATAATAAGAATCATTGTTATATCAattat                                               >  1:1263345/1‑150 (MQ=255)
                                                                                                           gCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAattattat                                            <  2:624555/150‑1 (MQ=255)
                                                                                                            ccTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAattattatt                                           <  1:790942/150‑1 (MQ=255)
                                                                                                             ctctGCGGTAGATACCCGGTTATGAATGCGTGCCCATAACTGACGGAATGGGCATTTCTTCGCAACTTGTTGTCTTCTCAAAAAATTACTGCTTGCTCTGGTCAGCCATAATATTTATAATAAGAATCATTGTTATATCAAttattatta                                          <  1:308713/150‑1 (MQ=255)
                                                                                                             ctctGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAAttattatta                                          <  1:295159/150‑1 (MQ=255)
                                                                                                             ctctGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAAttattatta                                          <  2:924415/150‑1 (MQ=255)
                                                                                                             ctctGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAAttattatta                                          >  1:866785/1‑150 (MQ=255)
                                                                                                                     tAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCTTAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTAt                                  <  1:1361371/150‑1 (MQ=255)
                                                                                                                      aGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATg                                 <  2:460720/150‑1 (MQ=255)
                                                                                                                       gATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATAATGATAATTAGAATCATTGTTATATAAATTATTATTAATTTTTATGc                                >  1:111638/1‑150 (MQ=255)
                                                                                                                        aTACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCg                               <  2:491298/150‑1 (MQ=255)
                                                                                                                        aTACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCg                               <  2:1360035/150‑1 (MQ=255)
                                                                                                                        aTACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCg                               <  2:1339885/150‑1 (MQ=255)
                                                                                                                        aTACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCg                               >  2:1668647/1‑150 (MQ=255)
                                                                                                                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGt                              <  1:1197553/150‑1 (MQ=255)
                                                                                                                          aCCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGTATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGtt                             <  2:1717044/150‑1 (MQ=255)
                                                                                                                          aCCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGtt                             >  2:18486/1‑150 (MQ=255)
                                                                                                                           cccGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTta                            <  2:744922/150‑1 (MQ=255)
                                                                                                                           cccGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTta                            >  1:1154588/1‑150 (MQ=255)
                                                                                                                           cccGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTta                            >  2:1561783/1‑150 (MQ=255)
                                                                                                                            ccGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTtat                           <  1:1645825/150‑1 (MQ=255)
                                                                                                                             cGCTTCTGAATGCGTGCCCATTCCTTACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTtata                          <  2:1149079/150‑1 (MQ=255)
                                                                                                                             cGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTtata                          <  2:555162/150‑1 (MQ=255)
                                                                                                                                  cTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATTATAAGAATCAGTGTTTTATCAATTTTTATTAATTTTTATGCGTTATACGGAt                     >  2:1432536/1‑150 (MQ=255)
                                                                                                                                    gAATGTGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAg                   >  1:506460/1‑150 (MQ=255)
                                                                                                                                      aTGCGTGCCCATTCCTGATGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCa                 <  1:1332639/150‑1 (MQ=255)
                                                                                                                                      aTGCGTGCCCATTCCTGACGGATTGGGCATTTCTGCGCCACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGGCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCa                 <  2:870664/150‑1 (MQ=255)
                                                                                                                                         cGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGaa              >  1:1383391/1‑150 (MQ=255)
                                                                                                                                         cGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGGCAGCCATTATATTGATATTTAGAATCATTGTTATATCAATTATTATGAATTTTTATGCGTTATAGGGATTGCAGaa              >  1:235941/1‑150 (MQ=255)
                                                                                                                                            gCCCATTCCTGACTGAATTGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc           <  1:491398/150‑1 (MQ=255)
                                                                                                                                            gCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc           <  2:591748/150‑1 (MQ=255)
                                                                                                                                            gCCCATTCCGTCAGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc           >  2:1259351/1‑150 (MQ=255)
                                                                                                                                             cccATTCCTTACGGAATGGGCATTTCTGCGCAACTTGTTGTATTCTCAACAAATTACTGCTTGCTCTGGTCAGACATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACt          <  1:345713/150‑1 (MQ=255)
                                                                                                                                               cATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACTAATTACTGCTTGCTCTGGTCAGCCATAATATTGTTAATAAGAATCAATGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCATAAAACTCa        >  1:1234826/1‑150 (MQ=255)
                                                                                                                                                     tGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTg  <  2:990296/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTG  >  NC_000913/1789348‑1789646

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: