Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,500 C→A 5.1% 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 5.1% 286.6 / 2.8 101intergenic (+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 (47/47);  new base A (2/3);  total (49/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

Only 100 of 101 total aligned reads displayed.

CGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTG  >  NC_000913/1789357‑1789646
                                                                                                                                               |                                                                                                                                                  
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                                          <  2:924415/150‑1 (MQ=255)
                                                                                                    ctctGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAAttattatta                                          <  1:295159/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: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)
                                                                                                                                               |                                                                                                                                                  
CGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTG  >  NC_000913/1789357‑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: