Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,499 A→C 7.0% intergenic (+18/‑114) 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,4990AC 7.0% 260.6 / 6.9 101intergenic (+18/‑114)aroH/ydiE3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible/hemin uptake protein HemP homolog
Reads supporting (aligned to +/- strand):  ref base A (44/49);  new base C (5/3);  total (49/52)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 102 total aligned reads displayed.

GAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTG  >  NC_000913/1789350‑1789646
                                                                                                                                                     |                                                                                                                                                   
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                                                                                                                                   <  2:1047061/150‑1 (MQ=255)
                  gtcgGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTgc                                                                                                                                   <  1:616946/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:207623/1‑150 (MQ=255)
                                               aTCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACt                                                                                                      >  2:1625790/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                                                                                          <  2:17972/150‑1 (MQ=255)
                                                           cccGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGt                                                                                          >  1:620146/1‑150 (MQ=255)
                                                           cccGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGt                                                                                          >  1:838599/1‑150 (MQ=255)
                                                           cccGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGt                                                                                          <  1:999040/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)
                                                                                                                   tAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCTTAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTAt                                  <  1:1361371/150‑1 (MQ=255)
                                                                                                                    aGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATg                                 <  2:460720/150‑1 (MQ=255)
                                                                                                                     gATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATAATGATAATTAGAATCATTGTTATATAAATTATTATTAATTTTTATGc                                >  1:111638/1‑150 (MQ=255)
                                                                                                                      aTACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCg                               <  2:1360035/150‑1 (MQ=255)
                                                                                                                      aTACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCg                               <  2:1339885/150‑1 (MQ=255)
                                                                                                                      aTACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCg                               <  2:491298/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)
                                                                                                                                        gTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGaaa             <  2:1345287/150‑1 (MQ=255)
                                                                                                                                          gCCCATTCCTGACTGAATTGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc           <  1:491398/150‑1 (MQ=255)
                                                                                                                                          gCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc           <  2:591748/150‑1 (MQ=255)
                                                                                                                                          gCCCATTCCGTCAGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc           >  2:1259351/1‑150 (MQ=255)
                                                                                                                                          gCCCATTCCGTCAGGAATGGGCATTTATGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAAc           >  1:419868/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)
                                                                                                                                                     |                                                                                                                                                   
GAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTG  >  NC_000913/1789350‑1789646

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: