Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 18,683 C→T 6.7% intergenic (+28/‑32) nhaA → / → nhaR Na(+):H(+) antiporter NhaA/DNA‑binding transcriptional activator NhaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091318,6830CT 6.7% 275.8 / 6.8 90intergenic (+28/‑32)nhaA/nhaRNa(+):H(+) antiporter NhaA/DNA‑binding transcriptional activator NhaR
Reads supporting (aligned to +/- strand):  ref base C (44/38);  new base T (2/4);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

CCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGAC  >  NC_000913/18534‑18832
                                                                                                                                                     |                                                                                                                                                     
ccTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAAc                                                                                                                                                       >  1:538542/1‑150 (MQ=255)
      gTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTcc                                                                                                                                                 >  2:599578/1‑150 (MQ=255)
         gCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCGTTTATGGCTccccg                                                                                                                                              <  2:400495/150‑3 (MQ=255)
         gCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGGTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCGTTTATGGCTccccg                                                                                                                                              <  1:1696072/150‑3 (MQ=255)
          cGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCGTTTATGGCTccccgg                                                                                                                                             <  1:1752751/150‑4 (MQ=255)
          cGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCttt                                                                                                                                             <  2:1486699/150‑1 (MQ=255)
            tAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTc                                                                                                                                           <  1:1930597/150‑1 (MQ=255)
            tAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTc                                                                                                                                           >  2:1669352/1‑150 (MQ=255)
                tCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTg                                                                                                                                       >  1:1480950/1‑150 (MQ=255)
                 ccAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTAACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGt                                                                                                                                      <  2:1476958/150‑1 (MQ=255)
                  cAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGtt                                                                                                                                     >  2:1334998/1‑150 (MQ=255)
                      aCTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCa                                                                                                                                 >  1:975942/1‑150 (MQ=255)
                       cTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCGTTTATGGCTccccg                                                                                                                                              >  1:1555510/1‑134 (MQ=255)
                       cTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCGTTTATGGCTccccg                                                                                                                                              <  2:1555510/136‑3 (MQ=255)
                       cTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAg                                                                                                                                >  2:1435768/1‑150 (MQ=255)
                        tGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAgg                                                                                                                               <  2:779469/150‑1 (MQ=255)
                              aCTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGgagaga                                                                                                                         <  1:386725/150‑1 (MQ=255)
                              aCTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGgagaga                                                                                                                         >  1:1825128/1‑150 (MQ=255)
                               cTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGaa                                                                                                                        >  2:313016/1‑150 (MQ=255)
                                 gggCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAAt                                                                                                                      >  1:1430275/1‑150 (MQ=255)
                                  ggCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTGGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATg                                                                                                                     >  2:1382072/1‑150 (MQ=255)
                                    cGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAg                                                                                                                   >  1:607221/1‑150 (MQ=255)
                                    cGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAg                                                                                                                   >  2:1081350/1‑150 (MQ=255)
                                      aaaCTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGcc                                                                                                                 <  2:642404/150‑2 (MQ=255)
                                        aCTGGGTATCGTGGTCGGTTCTATCTCTTCGGCGGTAATTGGTTACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATg                                                                                                               <  2:448739/150‑1 (MQ=255)
                                        aCTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGt                                                                                                                                      <  1:191093/127‑1 (MQ=255)
                                        aCTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGt                                                                                                                                      >  2:191093/1‑127 (MQ=255)
                                         cTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGt                                                                                                              >  2:768098/1‑150 (MQ=255)
                                          tCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGtc                                                                                                             <  2:1225492/150‑1 (MQ=255)
                                           cGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGtct                                                                                                            >  2:220851/1‑150 (MQ=255)
                                           cGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGtct                                                                                                            <  1:1669352/150‑1 (MQ=255)
                                           cGGTATCCTGGGCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGGTTGCGTCGATCAGGTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTGTTCATTGTTATCAGGGAGAGAAATGAGCATGtct                                                                                                            >  2:302290/1‑150 (MQ=255)
                                             gTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCa                                                                                                          <  1:1103065/150‑1 (MQ=255)
                                                   tGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACTGCTCCCTTTTCATTGTTATCAGGGATAGAAATGAGCATGTCTCATATcta                                                                                                    <  2:1193251/150‑3 (MQ=255)
                                                       cGGTTCTATCTCTTCGGCGGTAATGGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTAc                                                                                                <  2:23968/150‑1 (MQ=255)
                                                            cTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCa                                                                                           >  1:1900264/1‑150 (MQ=255)
                                                                ctctTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTg                                                                                       <  2:1824034/150‑1 (MQ=255)
                                                                 tctTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGt                                                                                      <  1:1334998/150‑1 (MQ=255)
                                                                    tCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTAtt                                                                                   <  2:550508/150‑1 (MQ=255)
                                                                      ggcggTACTTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGCGTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCGCATATCAATTACAACCACTTGTATTAc                                                                                 <  2:913966/150‑1 (MQ=255)
                                                                           tAATTGGATACAGCTGGTTACGCTTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACTGCTCCCTTTTCATTGTTATCAGTGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTg                                                                            <  1:250594/150‑1 (MQ=255)
                                                                           tAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTAATTGTTATAAGGGAGAGAAATGAGCATTTCTCATATAAATTACAACCACTTGTATTACTTCtt                                                                            >  1:1659192/1‑149 (MQ=255)
                                                                             aTTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGc                                                                          >  2:212846/1‑150 (MQ=255)
                                                                             aTTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGc                                                                          >  2:1405189/1‑150 (MQ=255)
                                                                                ggATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATg                                                                       <  2:607221/150‑1 (MQ=255)
                                                                                ggATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATg                                                                       >  2:1212469/1‑150 (MQ=255)
                                                                                ggATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCGTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATg                                                                       <  2:538542/150‑1 (MQ=255)
                                                                                  aTACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTc                                                                     >  1:991811/1‑150 (MQ=255)
                                                                                   tACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCt                                                                    >  1:854866/1‑150 (MQ=255)
                                                                                     cAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCtat                                                                  >  1:592127/1‑150 (MQ=255)
                                                                                     cAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCtat                                                                  <  2:991811/150‑1 (MQ=255)
                                                                                      aGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCtata                                                                 <  1:212846/150‑1 (MQ=255)
                                                                                       gCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATaa                                                                >  1:1047506/1‑150 (MQ=255)
                                                                                       gCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTACATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTGCAACCACTTGTATTACTTCTGGCATGTCTATaa                                                                >  1:1991190/1‑150 (MQ=255)
                                                                                         tGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAaag                                                              <  1:644436/150‑1 (MQ=255)
                                                                                             tACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAgg                                                          >  1:1649493/1‑150 (MQ=255)
                                                                                              aCGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGt                                                         >  1:1287381/1‑150 (MQ=255)
                                                                                               cgcgTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGtt                                                        >  2:1587881/1‑150 (MQ=255)
                                                                                                 cgttcgttCGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTcc                                                      >  1:998770/1‑150 (MQ=255)
                                                                                                   ttcgttTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGt                                                    >  2:1165026/1‑150 (MQ=255)
                                                                                                         tGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTgg                                              <  1:599578/150‑1 (MQ=255)
                                                                                                         tGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAAAGGCTCCCTTTTCATTGTTATAAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTgg                                              >  2:1367832/1‑150 (MQ=255)
                                                                                                          gCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGgc                                             >  1:837275/1‑150 (MQ=255)
                                                                                                          gCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGgc                                             <  1:1789451/150‑1 (MQ=255)
                                                                                                               cATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGc                                        >  1:1042951/1‑150 (MQ=255)
                                                                                                               cATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGc                                        >  1:1583713/1‑150 (MQ=255)
                                                                                                                aTCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCg                                       <  2:1825128/150‑1 (MQ=255)
                                                                                                                 cccgTTTGAGAGCTCTGTTTACCTGGGAGCCTTAAACGGCTCCCGGTTCATTGTTATCCGGGTCACAAATGTGCATGTGTGTTTTCAATTCCTTCCACTTGTATTTCCTCTGGCCTGTCTATAAAGACGGTTCCGTGGTTGGCGCAGCgg                                      <  2:1899012/147‑1 (MQ=255)
                                                                                                                    ttttGACAGGAGGGTTTACCGGGGAGCCATAAAAGGCTCCCTTTTCATTGTTATCAGTGTGAGAAATGAGCATGTCTCATATCAATTGCAACCACTTGTATTACTTCTGGCTTGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAgg                                   <  2:1991190/149‑1 (MQ=255)
                                                                                                                          cAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCttt                             >  1:424888/1‑150 (MQ=255)
                                                                                                                          cAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCttt                             <  2:475032/150‑1 (MQ=255)
                                                                                                                            ggACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttat                           <  2:746966/150‑1 (MQ=255)
                                                                                                                            ggACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCcttat                           >  2:610544/1‑150 (MQ=255)
                                                                                                                              aCGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttattt                         <  1:684540/150‑1 (MQ=255)
                                                                                                                                 gTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAAc                      <  2:1480950/150‑1 (MQ=255)
                                                                                                                                  tttACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACt                     <  2:1047506/150‑1 (MQ=255)
                                                                                                                                    tACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTcc                   >  1:1397187/1‑150 (MQ=255)
                                                                                                                                     aaaagggAGCCGTTTATGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCca                  >  1:574483/5‑150 (MQ=255)
                                                                                                                                     aCCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCca                  >  1:1746142/1‑150 (MQ=255)
                                                                                                                                     aCCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCca                  >  1:920171/1‑150 (MQ=255)
                                                                                                                                     aCCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCca                  <  2:1430275/150‑1 (MQ=255)
                                                                                                                                     aCCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCca                  >  2:861902/1‑150 (MQ=255)
                                                                                                                                      ccGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCcac                 >  1:1846914/1‑150 (MQ=255)
                                                                                                                                       cGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCcaca                >  1:444423/1‑150 (MQ=255)
                                                                                                                                       cGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCcaca                <  1:1382072/150‑1 (MQ=255)
                                                                                                                                          ggAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGTCTCAGCGGAGGCGCTTTATTTAACTCCACAAAc             <  2:1649493/150‑1 (MQ=255)
                                                                                                                                               cATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTa        >  1:546675/1‑150 (MQ=255)
                                                                                                                                                aTAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGACGCGCTTTATTTAACTCAACAAACCATTaa       >  2:1735802/1‑149 (MQ=255)
                                                                                                                                                 tAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTAcc      <  1:220851/150‑1 (MQ=255)
                                                                                                                                                  aaaCGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCg     <  1:1587881/150‑1 (MQ=255)
                                                                                                                                                    aCGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGa   <  1:1081350/150‑1 (MQ=255)
                                                                                                                                                     cGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGAc  <  1:1212469/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGAC  >  NC_000913/18534‑18832

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