Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,693,606 A→G intergenic (‑43/‑353) purL ← / → yfhD phosphoribosylformyl‑glycineamide synthetase/putative periplasmic binding transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,693,6060AG99.3% 472.3 / ‑4.5 150intergenic (‑43/‑353)purL/yfhDphosphoribosylformyl‑glycineamide synthetase/putative periplasmic binding transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (78/70);  total (78/72)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

Only 100 of 166 total aligned reads displayed.

GATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGC  >  NC_000913/2693516‑2693695
                                                                                           |                                                                                         
gATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGgaa                                                                                         <  2:3163346/93‑1 (MQ=255)
gATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGgaa                                                                                         <  1:4522538/93‑1 (MQ=255)
gATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGgaa                                                                                         <  1:4315602/93‑1 (MQ=255)
gATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGgaa                                                                                         <  1:6353716/93‑1 (MQ=255)
gATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGgaa                                                                                         <  2:5131333/93‑1 (MQ=255)
 aTTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAc                                                                                        <  2:626985/93‑1 (MQ=255)
 aTTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAc                                                                                        >  2:1841005/1‑93 (MQ=255)
  ttCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACg                                                                                       >  1:1363904/1‑93 (MQ=255)
   tCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTGGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACgg                                                                                      <  2:4821676/93‑1 (MQ=255)
     ggAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGc                                                                                    <  2:1080627/93‑1 (MQ=255)
      gAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCg                                                                                   <  1:2200992/93‑1 (MQ=255)
      gAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCg                                                                                   >  2:3745624/1‑93 (MQ=255)
      gAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCg                                                                                   >  2:6135680/1‑93 (MQ=255)
        aTGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTc                                                                                 <  2:4281912/93‑1 (MQ=255)
        aTGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTc                                                                                 >  1:317611/1‑93 (MQ=255)
        aTGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTc                                                                                 >  1:2399297/1‑93 (MQ=255)
         tGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCa                                                                                >  1:5446653/1‑93 (MQ=255)
         tGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCa                                                                                <  2:3370204/93‑1 (MQ=255)
          gCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCAt                                                                               <  2:4122107/93‑1 (MQ=255)
            cGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATta                                                                             <  1:4840467/93‑1 (MQ=255)
             gACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATtat                                                                            >  2:2727268/1‑93 (MQ=255)
             gACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATtat                                                                            >  2:3734548/1‑93 (MQ=255)
              aCAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATtata                                                                           <  1:5696640/93‑1 (MQ=255)
                aGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATaaa                                                                         >  1:1694102/1‑93 (MQ=255)
                 gTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAg                                                                        >  2:3426189/1‑93 (MQ=255)
                 gTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAg                                                                        >  1:3928392/1‑93 (MQ=255)
                 gTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAg                                                                        <  2:6614738/93‑1 (MQ=255)
                  tGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGa                                                                       >  1:5979718/1‑93 (MQ=255)
                  tGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGa                                                                       >  1:5719950/1‑93 (MQ=255)
                   gCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGaa                                                                      >  2:2852234/1‑93 (MQ=255)
                   gCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGaa                                                                      >  2:2210444/1‑93 (MQ=255)
                    cAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAAt                                                                     >  1:6617951/1‑93 (MQ=255)
                    cAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAAt                                                                     >  1:17144/1‑93 (MQ=255)
                      ggCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCt                                                                   >  2:3063539/1‑93 (MQ=255)
                       gCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTg                                                                  <  1:3277822/93‑1 (MQ=255)
                       gCGAACCACGCAGAATTTCCATCATAAGCCCCTGGTCTTCTAACCTTTCGGTGTCCCAAAGGGGGGGGAACGGCGTCATTATAAAGAATCTGa                                                                 <  2:872691/93‑1 (MQ=255)
                          aaCCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATg                                                               >  2:5341487/1‑93 (MQ=255)
                          aaCCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATg                                                               >  1:3947182/1‑93 (MQ=255)
                          aaCCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATg                                                               <  1:154324/93‑1 (MQ=255)
                           aCCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATgc                                                              >  2:5517723/1‑93 (MQ=255)
                           aCCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATgc                                                              >  1:2239715/1‑93 (MQ=255)
                              aCGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCt                                                           <  2:1689823/93‑1 (MQ=255)
                              aCGCAGAATTTCCATCATAAGTCTCTCGTCATCGAAGCTTTGGGTGCACCCAAGGGGGGTGAACGGGCGTCATTATAAAGAATCTGATGCGCt                                                           <  1:4935837/93‑1 (MQ=255)
                               cGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTg                                                          >  1:3146060/1‑93 (MQ=255)
                               cGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTg                                                          >  1:2970822/1‑93 (MQ=255)
                                   gAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGa                                                      >  2:1360553/1‑93 (MQ=255)
                                     aTTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGaaa                                                    <  1:937934/93‑1 (MQ=255)
                                       ttCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAAcc                                                  <  1:3759036/93‑1 (MQ=255)
                                       ttCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAAcc                                                  <  2:1455562/93‑1 (MQ=255)
                                        tCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCg                                                 <  2:5150298/93‑1 (MQ=255)
                                         ccATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGt                                                <  2:6224395/93‑1 (MQ=255)
                                         ccATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGt                                                <  1:6274167/93‑1 (MQ=255)
                                          catcatAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGACTCTGATGCGCTGACGAAACCgct                                               >  1:1774478/1‑91 (MQ=255)
                                          catcatAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGtt                                               <  2:4920850/93‑1 (MQ=255)
                                          catcatAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCCTCATTATAAAGAATCTGATGCGCTGACGAAACCGtt                                               <  2:2676733/93‑1 (MQ=255)
                                           atcatAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGttt                                              >  2:1398260/1‑93 (MQ=255)
                                            tcatAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTg                                             <  2:5099990/93‑1 (MQ=255)
                                            tcatAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTg                                             >  2:6006612/1‑93 (MQ=255)
                                             catAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGc                                            <  2:102381/93‑1 (MQ=255)
                                              atAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCg                                           >  1:1942227/1‑93 (MQ=255)
                                               tAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGt                                          >  2:5439845/1‑93 (MQ=255)
                                                aaGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTg                                         <  1:1604752/93‑1 (MQ=255)
                                                 aGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTgg                                        <  2:2542183/93‑1 (MQ=255)
                                                 aGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTgg                                        >  1:3669726/1‑93 (MQ=255)
                                                 aGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTgg                                        >  1:2528870/1‑93 (MQ=255)
                                                  gTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTTATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGa                                       >  1:1213443/1‑93 (MQ=255)
                                                  gTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGa                                       <  1:2922699/93‑1 (MQ=255)
                                                         gTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATaaaa                                >  1:5077607/1‑93 (MQ=255)
                                                         gTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATaaaa                                <  1:5044282/93‑1 (MQ=255)
                                                         gTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATaaaa                                >  2:4066680/1‑93 (MQ=255)
                                                         gTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGGGGAAATaaaa                                >  2:3867588/1‑93 (MQ=255)
                                                           cttctAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGTAATAAAATc                              >  1:712983/1‑93 (MQ=255)
                                                           cttctAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATc                              <  2:1944571/93‑1 (MQ=255)
                                                            ttctAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCTTTTCCGGGGAAAAAAAAtac                             >  1:2002781/1‑91 (MQ=255)
                                                            ttctAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCa                             >  1:2160258/1‑93 (MQ=255)
                                                            ttctAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCa                             >  2:240059/1‑93 (MQ=255)
                                                             tctAAGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCAc                            <  1:1669342/93‑1 (MQ=255)
                                                                aaGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCAt                         <  1:6460129/93‑1 (MQ=255)
                                                                aaGCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCAt                         <  2:2694343/93‑1 (MQ=255)
                                                                  gCTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCg                       >  1:2737379/1‑93 (MQ=255)
                                                                   cTTTCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGt                      >  1:3082806/1‑93 (MQ=255)
                                                                      tCGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGaa                   <  1:3336831/93‑1 (MQ=255)
                                                                       cGGTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAAt                  <  1:2369838/93‑1 (MQ=255)
                                                                        ggTGTACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAAtt                 <  1:5575280/93‑1 (MQ=255)
                                                                         gtgtACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTa                >  1:3809115/1‑93 (MQ=255)
                                                                          tgtACCCAAGGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTa                <  2:6012997/93‑1 (MQ=255)
                                                                          tgtACCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAg               >  1:805595/1‑93 (MQ=255)
                                                                             aCCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCaa            >  2:4582341/1‑93 (MQ=255)
                                                                             aCCCAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCaa            <  2:6694827/93‑1 (MQ=255)
                                                                              cccAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAAc           >  1:5320654/1‑93 (MQ=255)
                                                                              cccAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAAc           >  2:4556102/1‑93 (MQ=255)
                                                                              cccAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAAc           <  2:2972437/93‑1 (MQ=255)
                                                                               ccAAGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAAcg          >  1:2317261/1‑93 (MQ=255)
                                                                                 aaGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAAcgcg        >  1:2784622/1‑93 (MQ=255)
                                                                                 aaGGGGGGGGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAAcgcg        >  1:662789/1‑93 (MQ=255)
                                                                                    ggggggggAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTg      >  2:4794715/1‑93 (MQ=255)
                                                                                    ggggggggAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTg      <  2:6445454/93‑1 (MQ=255)
                                                                                      ggggggAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTgcc    <  2:5859689/93‑1 (MQ=255)
                                                                                      ggggggAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTgcc    <  2:3594709/93‑1 (MQ=255)
                                                                                       ggggGAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTgccg   <  1:5519527/93‑1 (MQ=255)
                                                                                           |                                                                                         
GATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGC  >  NC_000913/2693516‑2693695

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: