Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,159,271 G→T G42V (GGC→GTC)  fabR → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,159,2710GT99.8% 1937.6 / NA 518G42V (GGC→GTC) fabRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (244/272);  total (245/273)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 539 total aligned reads displayed.

AAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATG  >  NC_000913/4159173‑4159368
                                                                                                  |                                                                                                 
aaaacCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCa                                                                                                 <  1:3785025/101‑1 (MQ=255)
 aaacCCCGCCGTTGACTGGTGGAAGCCGAATTAAACCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAt                                                                                                <  2:229723/101‑1 (MQ=255)
  ataCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGAGCTGAACACAGCTTCCCCAGCCTGAGTTTGCGTGAAGAGGAGCGAGAACCCGTCAtt                                                                                               >  2:5893523/3‑101 (MQ=255)
  aaccccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAtt                                                                                               <  1:11095598/98‑1 (MQ=255)
  aaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAtt                                                                                               <  1:12489376/101‑1 (MQ=255)
   nnnccgccgTTCGCTGGTGGAAGCCGAATTTAGCCAATTAACTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTg                                                                                              <  2:10549581/98‑1 (MQ=255)
    ccccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGc                                                                                             <  1:74567/100‑1 (MQ=255)
     cccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCt                                                                                            <  2:5104149/101‑1 (MQ=255)
      nnnccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGccc                                                                                           <  2:7789783/98‑3 (MQ=255)
      ccgccgTTCGCTGGTGTAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTc                                                                                           <  2:564203/101‑1 (MQ=255)
      ccgccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTc                                                                                           <  1:14038548/101‑1 (MQ=255)
       nnncgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCtac                                                                                          <  2:11355423/98‑3 (MQ=255)
       cgccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTcc                                                                                          <  1:14743801/101‑1 (MQ=255)
         ccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCa                                                                                        >  2:690778/1‑101 (MQ=255)
          nnntCGCGGGTGTAATGCGCATTTAGCCAATTAAGTGCTGCACGCAGCATCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAc                                                                                       <  2:11506598/98‑1 (MQ=255)
            tgcGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGACGCGGTCATTGCTCCCannn                                                                                     >  2:7599919/3‑98 (MQ=255)
              cGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACctct                                                                                   >  1:2470235/1‑101 (MQ=255)
              cGCTGGTGGAAGCCGCATTTAGCCAATAAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACctct                                                                                   <  1:321360/101‑1 (MQ=255)
                cTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCttt                                                                                 <  2:818872/101‑1 (MQ=255)
                cTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCttt                                                                                 <  2:3442013/101‑1 (MQ=255)
                 tggtggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtttt                                                                                >  2:3093543/1‑101 (MQ=255)
                      nnngCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCg                                                                           <  2:10143939/98‑1 (MQ=255)
                      gAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCg                                                                           <  1:9077606/101‑1 (MQ=255)
                      gAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCg                                                                           >  1:3072872/1‑101 (MQ=255)
                       aaGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCgg                                                                          <  1:7961151/101‑1 (MQ=255)
                       aaGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCgg                                                                          >  2:5354903/1‑101 (MQ=255)
                        aGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGc                                                                         <  1:5011081/101‑1 (MQ=255)
                          ccGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                                       >  1:4792700/1‑101 (MQ=255)
                          ccGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                                       <  1:13475331/101‑1 (MQ=255)
                           cGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAtt                                                                      >  2:1500740/1‑101 (MQ=255)
                            nnntttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAttt                                                                     <  2:8441246/98‑1 (MQ=255)
                            gCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCGttt                                                                     >  2:63571/1‑101 (MQ=255)
                            gCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAttt                                                                     >  2:424581/1‑101 (MQ=255)
                            gCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCGCTGTTTATCCGCAttt                                                                     >  1:7172183/1‑101 (MQ=255)
                            gCATTGAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAttt                                                                     <  1:6231388/101‑1 (MQ=255)
                             cATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                                    >  1:6046051/1‑101 (MQ=255)
                             cATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                                    <  2:281930/101‑1 (MQ=255)
                             cATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                                    >  1:4147883/1‑101 (MQ=255)
                             cATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                                    <  1:1172405/101‑1 (MQ=255)
                               tttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                                                                  >  1:13770530/1‑101 (MQ=255)
                                ttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcgc                                                                 <  2:3914143/101‑1 (MQ=255)
                                    ccAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACg                                                             <  1:6513666/101‑1 (MQ=255)
                                      nnntAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTa                                                           <  2:14231014/98‑1 (MQ=255)
                                        ttAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGa                                                         <  1:3646052/101‑1 (MQ=255)
                                         tAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGAc                                                        <  2:1776519/101‑1 (MQ=255)
                                          aaGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACg                                                       >  1:10981229/1‑101 (MQ=255)
                                              gCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGannn                                                   >  2:9209978/1‑98 (MQ=255)
                                                tGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACttg                                                 >  2:2749272/1‑99 (MQ=255)
                                                tGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTgg                                                 >  1:765154/1‑101 (MQ=255)
                                                  aaCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGt                                               <  1:1986476/101‑1 (MQ=255)
                                                  aaCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGt                                               >  1:8597283/1‑101 (MQ=255)
                                                  aaCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCGCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGTCGAACTgggc                                               >  2:4856142/1‑100 (MQ=255)
                                                   aCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTc                                              <  2:7176724/101‑1 (MQ=255)
                                                     gCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTg                                            >  2:441330/1‑101 (MQ=255)
                                                      cAGCTTCGCCAGCTTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                                   >  2:13849056/1‑93 (MQ=255)
                                                      cAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGa                                           >  2:6968246/1‑101 (MQ=255)
                                                       aGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCCTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCtnnn                                          >  2:13091583/1‑98 (MQ=255)
                                                         cTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCa                                        <  1:10005960/101‑1 (MQ=255)
                                                          ttCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCAt                                       <  1:10297498/101‑1 (MQ=255)
                                                          nnngCCAGCCTGAGTTTGCGTGAAGTGAAGCGTAAAGCGGTCATCGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCAt                                       <  2:8794288/98‑1 (MQ=255)
                                                           tCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATg                                      >  1:6679154/1‑101 (MQ=255)
                                                             gCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGt                                    <  2:241554/101‑1 (MQ=255)
                                                             gCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTATGACCATGGt                                    >  2:6724417/1‑101 (MQ=255)
                                                              nnngCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGtt                                   <  2:12479340/98‑1 (MQ=255)
                                                              ccAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATgnnn                                   >  2:11872392/1‑98 (MQ=255)
                                                              ccAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGtt                                   >  1:3412122/1‑101 (MQ=255)
                                                               cAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtg                                  >  1:1689922/1‑101 (MQ=255)
                                                               cAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtg                                  <  1:11352439/101‑1 (MQ=255)
                                                                 gCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgat                                >  1:6874048/1‑101 (MQ=255)
                                                                   nnnaGTGTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGTCGTAGACGAACTGGGTCTGACCATGGTtgatga                              <  2:14169555/98‑1 (MQ=255)
                                                                    taaGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGCCGAACTGGGTCTGACCATGGGtgatnnn                             >  2:11107775/3‑98 (MQ=255)
                                                                    tGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag                             >  1:9704614/1‑101 (MQ=255)
                                                                     gAGTTTGCGTGAAGTGGCGCGTGAAGCGGTTATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgaga                            >  2:4729156/1‑101 (MQ=255)
                                                                     gAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgagg                            >  1:1464807/1‑100 (MQ=255)
                                                                      cgTTTGCGAGAGGTGGCTCGTGATGGGGTGATTACTCCCCCCTCACTTTATCGGCCTTTCCGAGACGTAGACGAACTGGGTATGCCCATGGTTGATGagag                           <  2:1934330/100‑1 (MQ=255)
                                                                      aGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGagag                           >  1:8266162/1‑101 (MQ=255)
                                                                       gTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGc                          <  2:6189293/101‑1 (MQ=255)
                                                                         ttGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCgg                        <  1:12761757/101‑1 (MQ=255)
                                                                          tGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACCGGGTCTGACCATGGTTGATGAGAGCGGt                       <  1:14317614/101‑1 (MQ=255)
                                                                           gCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGtt                      <  1:5305382/101‑1 (MQ=255)
                                                                            nnngAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGttt                     <  2:12834655/98‑1 (MQ=255)
                                                                            cGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGttt                     <  1:12203689/101‑1 (MQ=255)
                                                                            cGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGttt                     <  1:7067213/101‑1 (MQ=255)
                                                                             gTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTa                    <  1:5523740/101‑1 (MQ=255)
                                                                              tGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTaa                   >  1:948992/1‑101 (MQ=255)
                                                                                  gTGGCGCTTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCt               <  1:8844498/101‑1 (MQ=255)
                                                                                  gTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCt               <  1:4517386/101‑1 (MQ=255)
                                                                                   tGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTa              >  2:5206854/1‑101 (MQ=255)
                                                                                    ggCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTAc             >  1:10768266/1‑101 (MQ=255)
                                                                                    ggCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTAc             >  2:971757/1‑101 (MQ=255)
                                                                                     nnncgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACg            <  2:7440530/98‑1 (MQ=255)
                                                                                     gcgcGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACg            >  1:2143437/1‑101 (MQ=255)
                                                                                        cgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGCGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCa         <  2:133393/101‑1 (MQ=255)
                                                                                         gTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCaa        >  1:1304596/1‑101 (MQ=255)
                                                                                          tGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAAc       >  2:1079129/1‑101 (MQ=255)
                                                                                          tGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAAc       <  1:13293244/101‑1 (MQ=255)
                                                                                           gAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACt      >  1:12006130/1‑101 (MQ=255)
                                                                                             aGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCa    >  1:1042268/1‑101 (MQ=255)
                                                                                               cGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGTGTCTGACCATGTTTGATGAGAGCGGTTTAATGCTACGCCAACTCAtc  >  1:5731313/1‑100 (MQ=255)
                                                                                               cGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATg  <  1:9953771/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
AAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATG  >  NC_000913/4159173‑4159368

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