Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,434,437 T→A 23.8% W109R (TGG→AGG)  ytfH → DUF24 family HxlR‑type putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,434,4370TA23.8% 297.0 / 97.6 155W109R (TGG→AGG) ytfHDUF24 family HxlR‑type putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (55/62);  new base A (20/17);  total (75/80)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 158 total aligned reads displayed.

CTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTT  >  NC_000913/4434338‑4434536
                                                                                                   |                                                                                                   
cTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACAGGCGGACt                                                                                                     >  1:3021680/1‑100 (MQ=255)
 ttAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTg                                                                                                    <  2:1721803/100‑1 (MQ=255)
   aaCCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAgg                                                                                                   >  1:410954/1‑99 (MQ=255)
    gatctATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGAt                                                                                                 <  2:3134297/96‑1 (MQ=255)
    aCCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGAt                                                                                                 <  1:603917/100‑1 (MQ=255)
    aCCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGAt                                                                                                 >  2:1779002/1‑100 (MQ=255)
     ccGTATCGCGTATCCGGTGGTGCGGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGAtt                                                                                                <  1:186489/100‑1 (MQ=255)
     ccGTAGCGTGTATCCGGTGGGGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGAt                                                                                                 <  1:1002930/99‑1 (MQ=255)
      cGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTg                                                                                               <  2:236820/100‑1 (MQ=255)
       gTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGa                                                                                              <  2:37609/100‑1 (MQ=255)
         aTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGt                                                                                            >  2:549582/1‑100 (MQ=255)
         aTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGt                                                                                            <  1:3161481/100‑1 (MQ=255)
         aTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGt                                                                                            <  2:2441828/100‑1 (MQ=255)
         aTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGAtagatc                                                                                            >  1:3134270/1‑95 (MQ=255)
         aTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGt                                                                                            <  2:2483489/100‑1 (MQ=255)
          ncgcgTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGtt                                                                                           >  2:3268169/2‑100 (MQ=255)
           cgcgTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTg                                                                                          <  1:1508028/100‑1 (MQ=255)
           cgcgTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTg                                                                                          <  2:3195045/100‑1 (MQ=255)
            gcgTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGa                                                                                         >  1:3331123/1‑100 (MQ=255)
              gTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAAt                                                                                       >  2:2813803/1‑100 (MQ=255)
               tATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAAtt                                                                                      <  1:926607/100‑1 (MQ=255)
               tATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAAtt                                                                                      <  1:2249518/100‑1 (MQ=255)
                 tCCGGTGGTGCCGCCGCATGTGGAATATGGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTg                                                                                    <  2:1188608/100‑1 (MQ=255)
                 tCCGGTGGTGCCGCCGCATGTGGAATATAGCCGCACGCCGCGGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTg                                                                                    <  1:594532/100‑1 (MQ=255)
                  ccGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGc                                                                                   <  1:48439/100‑1 (MQ=255)
                   cGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGcc                                                                                  <  1:1936420/100‑1 (MQ=255)
                       ggtgCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCTGAg                                                                              >  2:1623221/1‑100 (MQ=255)
                       ggtgCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAgg                                                                             >  1:2178613/1‑100 (MQ=255)
                        gtgCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAgg                                                                             <  1:3083425/100‑1 (MQ=255)
                         ngccgccGCATGTGGAATATAGCCTCACGCCGCTGGGCTAGCAGGTTAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGgt                                                                            >  1:155927/2‑100 (MQ=255)
                          gccgccGCATGTGGAATATAGCCTCACGCCGCTAGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGgtg                                                                           >  1:2042327/1‑100 (MQ=255)
                            cgccgcATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGtt                                                                         >  1:994596/1‑100 (MQ=255)
                                  aTGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGt                                                                   <  1:224821/100‑1 (MQ=255)
                                   tgtgGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTg                                                                  >  1:2008165/1‑100 (MQ=255)
                                   tgtgGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTg                                                                  <  1:2508893/100‑1 (MQ=255)
                                   tgtgGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTg                                                                  <  1:1162907/100‑1 (MQ=255)
                                    gtgGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCgg                                                                    <  2:2359585/97‑1 (MQ=255)
                                      ggAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGTTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCgg                                                               <  2:886911/100‑1 (MQ=255)
                                        aaTATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGa                                                             <  2:2646921/100‑1 (MQ=255)
                                        aaTATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGa                                                             <  2:972110/100‑1 (MQ=255)
                                         atatAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGAt                                                            <  2:2143980/100‑1 (MQ=255)
                                          tataGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATg                                                           >  2:1879906/1‑100 (MQ=255)
                                           ataGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGa                                                          >  2:2022913/1‑100 (MQ=255)
                                            taGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGaa                                                         <  2:1060168/100‑1 (MQ=255)
                                             atccTCACGCCGCAGGGCGAGCAGGTGAGCGGAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAAc                                                        >  2:2280422/3‑100 (MQ=255)
                                             aGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAAc                                                        >  1:709753/1‑100 (MQ=255)
                                              gCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACg                                                       >  2:1578211/1‑100 (MQ=255)
                                              gCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACg                                                       >  2:1337374/1‑100 (MQ=255)
                                              gCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACg                                                       <  2:2650728/100‑1 (MQ=255)
                                               ccTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGAGTTGGCGGTGCGGGATGAACGt                                                      >  2:1418125/1‑100 (MQ=255)
                                                cTCACGCCGCTGGGCGCGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGt                                                      >  1:3124876/1‑99 (MQ=255)
                                                cTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTg                                                     <  2:560269/100‑1 (MQ=255)
                                                 tCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGc                                                    <  1:1060092/100‑1 (MQ=255)
                                                   aCGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCgg                                                  <  2:466291/100‑1 (MQ=255)
                                                    cgccgcTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACTTGCGGc                                                 <  1:75309/100‑1 (MQ=255)
                                                    cgccgcTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGc                                                 <  2:620498/100‑1 (MQ=255)
                                                        gcTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATaa                                             >  2:1147132/1‑100 (MQ=255)
                                                        gcTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATaa                                             >  2:2614696/1‑100 (MQ=255)
                                                         cTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAAc                                            <  2:527525/100‑1 (MQ=255)
                                                            ggCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGc                                         >  2:1275499/1‑100 (MQ=255)
                                                              cGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGGGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCg                                       <  1:277308/100‑1 (MQ=255)
                                                              cGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCg                                       >  1:2915108/1‑100 (MQ=255)
                                                              cGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCg                                       <  1:277202/100‑1 (MQ=255)
                                                              cGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCg                                       <  1:3291778/100‑1 (MQ=255)
                                                              cGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAATGGATTGAGTTGTATTTGCCCGAGGTGTTGGGGGTGCGGGCTGAACGTGCGGCATAACTGcc                                        >  2:828894/1‑99 (MQ=255)
                                                                aGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGa                                     <  1:385674/100‑1 (MQ=255)
                                                                aGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGACCGTGCGGAATAACTGCCGGa                                     <  1:2197046/100‑1 (MQ=255)
                                                                 gCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGAt                                    >  1:3244375/1‑100 (MQ=255)
                                                                  cAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATg                                   >  1:441957/1‑100 (MQ=255)
                                                                   aGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATgc                                  <  2:2309430/100‑1 (MQ=255)
                                                                    ggTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATgcg                                 >  2:40716/1‑100 (MQ=255)
                                                                    ggTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATgcg                                 >  1:1389559/1‑100 (MQ=255)
                                                                    ggTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATgcg                                 >  1:2666674/1‑100 (MQ=255)
                                                                    ggTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCACGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATgcg                                 <  1:1674177/100‑1 (MQ=255)
                                                                     gTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATgcgc                                <  1:2489826/100‑1 (MQ=255)
                                                                     gTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATgcgc                                >  1:427998/1‑100 (MQ=255)
                                                                     gTGAGCGAAAAGGTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATgcgc                                <  2:2075596/100‑1 (MQ=255)
                                                                          cGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGc                           <  2:1914828/100‑1 (MQ=255)
                                                                          cGAAAAGGTTGCTGCACTGGCGGACGGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGc                           <  1:3160688/100‑1 (MQ=255)
                                                                           gAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCt                          <  2:387400/100‑1 (MQ=255)
                                                                                ggTTGCTGCACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCt                     <  2:205856/100‑1 (MQ=255)
                                                                                  ttGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTACCTgg                   >  1:3349007/1‑100 (MQ=255)
                                                                                    gctgcACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTc                 >  1:2769968/1‑100 (MQ=255)
                                                                                    gctgcACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTc                 >  2:709881/1‑100 (MQ=255)
                                                                                     ctgcACTGGCGGACTGGGTTGAGTTGAAGTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCt                <  1:1875669/100‑1 (MQ=255)
                                                                                     ctgcACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCt                >  1:488864/1‑100 (MQ=255)
                                                                                     ctgcACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCt                >  2:617179/1‑100 (MQ=255)
                                                                                     ctgcACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCt                >  1:2472556/1‑100 (MQ=255)
                                                                                     ctgcACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCt                >  2:2199067/1‑100 (MQ=255)
                                                                                      tgcACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTa               >  2:854905/1‑100 (MQ=255)
                                                                                       gcACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTAc              >  2:837520/1‑100 (MQ=255)
                                                                                       gcACTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTAc              >  2:2248705/1‑100 (MQ=255)
                                                                                         tcTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACgg            <  2:1355793/99‑1 (MQ=255)
                                                                                         aCTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACgg            <  1:1114325/100‑1 (MQ=255)
                                                                                          cTGGCGGACAGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGt           >  1:77437/1‑100 (MQ=255)
                                                                                            ggCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGc         >  2:317603/1‑100 (MQ=255)
                                                                                            ggCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGc         >  2:163439/1‑100 (MQ=255)
                                                                                             gCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCa        <  2:2808452/100‑1 (MQ=255)
                                                                                                 aCTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGc    <  1:1701428/100‑1 (MQ=255)
                                                                                                 aCTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGAGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATTTGGTCTACGGTGCATGGc    >  2:3329953/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTT  >  NC_000913/4434338‑4434536

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