Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,206 G→T 7.8% intergenic (+150/+48) yphH → / ← glyA putative DNA‑binding transcriptional regulator, NAGC‑like/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,2060GT 7.8% 252.0 / 8.8 77intergenic (+150/+48)yphH/glyAputative DNA‑binding transcriptional regulator, NAGC‑like/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (36/35);  new base T (2/4);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

TGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTCACACATCC  >  NC_000913/2684057‑2684354
                                                                                                                                                     |                                                                                                                                                    
tGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                                      >  2:1043992/1‑150 (MQ=255)
 ggCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                                     <  2:444474/150‑1 (MQ=255)
     cgATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                 <  1:662395/150‑1 (MQ=37)
       aTAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                               >  1:6647/1‑150 (MQ=255)
         aaCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                             <  2:1324149/150‑2 (MQ=25)
         aaCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                             <  1:1385105/150‑2 (MQ=25)
          aCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                            <  1:715145/150‑1 (MQ=255)
           cGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                                           >  1:673804/1‑150 (MQ=255)
           cGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                                           <  2:452471/150‑1 (MQ=255)
            gTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                                                          <  1:95587/150‑5 (MQ=25)
             catcaaGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                                                                                                                                         >  2:379132/5‑150 (MQ=255)
               gAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaaca                                                                                                                                       <  1:105102/150‑1 (MQ=255)
                caaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                                                      >  1:1290338/2‑150 (MQ=255)
                    gCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAg                                                                                                                                  >  1:217048/1‑150 (MQ=255)
                      ttCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca                                                                                                                                <  2:6647/150‑1 (MQ=255)
                      ttCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGAAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca                                                                                                                                <  1:378106/150‑1 (MQ=255)
                        cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca                                                                                                                              >  2:1004482/1‑150 (MQ=255)
                            aaggaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa                                                                                                                          <  1:1419007/150‑1 (MQ=255)
                             aggaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGAc                                                                                                                         <  2:1369130/150‑1 (MQ=255)
                                aagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGc                                                                                                                      <  1:532796/150‑1 (MQ=255)
                                 agCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCa                                                                                                                     <  1:1064590/150‑1 (MQ=255)
                                   ccTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCaaa                                                                                                                   >  1:544673/1‑150 (MQ=255)
                                         tgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCg                                                                                                             >  2:108736/1‑150 (MQ=255)
                                            tCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGttt                                                                                                          >  1:1042024/1‑150 (MQ=255)
                                                aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCt                                                                                                      <  1:1411128/150‑1 (MQ=255)
                                                 aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATATGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATAACCGTTTCGCtt                                                                                                     >  1:1149388/1‑150 (MQ=255)
                                                 aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCtt                                                                                                     >  2:739022/1‑150 (MQ=255)
                                                 aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCtt                                                                                                     >  2:173607/1‑150 (MQ=255)
                                                  aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTa                                                                                                    <  2:673804/150‑1 (MQ=255)
                                                   gaaaaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTAt                                                                                                   <  2:917263/148‑1 (MQ=255)
                                                   cGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTAt                                                                                                   >  1:611377/1‑150 (MQ=255)
                                                   cGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTAt                                                                                                   >  1:725575/1‑150 (MQ=255)
                                                   cGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTAt                                                                                                   >  1:142546/1‑150 (MQ=255)
                                                     aaaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGc                                                                                                 <  1:1004482/150‑1 (MQ=255)
                                                      aaaaTTGCCTGATGCGCTACGCTTATCAGTCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCg                                                                                                <  1:379132/150‑1 (MQ=255)
                                                      aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCg                                                                                                <  2:838403/150‑1 (MQ=255)
                                                               tGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCggg                                                                                       <  1:708545/150‑1 (MQ=255)
                                                                 aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTa                                                                                     <  2:175648/150‑1 (MQ=255)
                                                                      cTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTg                                                                                <  2:217048/150‑1 (MQ=255)
                                                                       tACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTgc                                                                               <  1:739022/150‑1 (MQ=255)
                                                                         cGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTgcgc                                                                             >  1:852815/1‑150 (MQ=255)
                                                                          gCTTATCAGGCCTACACGGGGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCa                                                                            >  2:1251867/1‑150 (MQ=255)
                                                                           cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAg                                                                           >  2:371819/1‑150 (MQ=255)
                                                                           cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAg                                                                           >  2:994453/1‑150 (MQ=255)
                                                                                     cTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATg                                                                        >  2:440094/1‑143 (MQ=255)
                                                                                     cTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATg                                                                        <  1:440094/143‑1 (MQ=255)
                                                                                      tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCAGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAAc                                                                >  2:67131/1‑150 (MQ=255)
                                                                                          cGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGTACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTa                                                            <  1:473158/150‑1 (MQ=255)
                                                                                          cGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTa                                                            >  2:415128/1‑150 (MQ=255)
                                                                                            cagaatccAATATATTGAAATTCCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGCCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTAcc                                                          >  2:775496/8‑150 (MQ=38)
                                                                                                 tGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGa                                                     <  1:173607/150‑1 (MQ=255)
                                                                                                  gCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                                    >  2:1360460/1‑150 (MQ=255)
                                                                                                   cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATg                                                   <  1:1366712/150‑1 (MQ=255)
                                                                                                    aaTATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATgc                                                  <  1:1134912/150‑1 (MQ=255)
                                                                                                        tatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTc                                              <  2:897322/150‑1 (MQ=255)
                                                                                                           tGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGAt                                           <  1:1043992/150‑1 (MQ=255)
                                                                                                               tttGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACg                                       <  1:763401/150‑1 (MQ=255)
                                                                                                               tttGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACg                                       <  1:758297/150‑1 (MQ=255)
                                                                                                                 tGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGc                                     >  2:1149927/1‑150 (MQ=255)
                                                                                                                  gCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCt                                    >  1:846333/1‑150 (MQ=255)
                                                                                                                  gCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCt                                    >  1:1373254/1‑150 (MQ=255)
                                                                                                                  gCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCt                                    >  1:498305/1‑150 (MQ=255)
                                                                                                                       gTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGTGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtcatc                               >  2:53303/1‑150 (MQ=255)
                                                                                                                       gTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtcatc                               >  2:235139/1‑150 (MQ=255)
                                                                                                                       gTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGATtcatc                               >  1:85149/1‑150 (MQ=255)
                                                                                                                        ttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtcatca                              <  2:725575/150‑1 (MQ=255)
                                                                                                                         tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTcatcat                             <  1:974608/150‑1 (MQ=255)
                                                                                                                          ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                            <  1:1282737/150‑1 (MQ=255)
                                                                                                                            gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACTGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGa                          >  2:1350536/1‑150 (MQ=255)
                                                                                                                             tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGAt                         <  2:498305/150‑1 (MQ=255)
                                                                                                                              atgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATg                        >  1:691402/3‑150 (MQ=39)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGc                       <  1:1360460/150‑1 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGc                       <  2:852815/150‑1 (MQ=255)
                                                                                                                                 ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTg                     <  1:371819/150‑1 (MQ=255)
                                                                                                                                        cggCGTGAACGCCTTATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCa              >  2:685727/2‑150 (MQ=255)
                                                                                                                                            gTTCATGCCGGATGCGGCGTGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTc          >  1:282045/1‑150 (MQ=255)
                                                                                                                                            gTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTc          >  1:395058/1‑150 (MQ=255)
                                                                                                                                               cACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTcaca       >  1:414800/1‑150 (MQ=255)
                                                                                                                                               cACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTcaca       >  1:425384/1‑150 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTcacaca     >  2:495493/1‑150 (MQ=255)
                                                                                                                                                   ccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTCACACATc   <  1:67131/150‑1 (MQ=255)
                                                                                                                                                    cgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGATAACTTTTCCTTTGATGCGCTCGATAACGTATTCATCATTTATGCTGTCCAGCACGTGACTCAttc  <  2:142546/150‑3 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTCACACATCC  >  NC_000913/2684057‑2684354

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