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.3% 1029.6 / NA 275G42V (GGC→GTC) fabRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (123/150);  total (125/150)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

AAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATG  >  NC_000913/4159173‑4159368
                                                                                                  |                                                                                                 
aaaacCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCa                                                                                                 <  2:2065140/101‑1 (MQ=255)
aaaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCa                                                                                                 <  1:1000720/101‑1 (MQ=255)
 aaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAt                                                                                                >  2:5332643/1‑101 (MQ=255)
  acccccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAtt                                                                                               <  2:1415720/98‑1 (MQ=255)
  aacccccccGTTCGCTGGTGGAAGCCGCATTTAGCCAATGAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAtt                                                                                               <  2:1808992/98‑1 (MQ=255)
  aaccccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAtt                                                                                               <  1:341901/98‑1 (MQ=255)
  aaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAtt                                                                                               <  2:5946720/101‑1 (MQ=255)
   aaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTg                                                                                              >  1:3256204/1‑101 (MQ=255)
     cccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCt                                                                                            <  1:4126931/101‑1 (MQ=255)
     cccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCt                                                                                            <  1:1594522/101‑1 (MQ=255)
      ccgccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTc                                                                                           <  2:2453242/101‑1 (MQ=255)
        gccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTccc                                                                                         <  2:4304504/101‑1 (MQ=255)
        gccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTccc                                                                                         >  2:1226997/1‑101 (MQ=255)
         ccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCa                                                                                        <  1:6859987/101‑1 (MQ=255)
         ccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCa                                                                                        >  1:5721168/1‑101 (MQ=255)
         ccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCa                                                                                        >  1:1156957/1‑101 (MQ=255)
           gTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAcc                                                                                      >  2:3650188/1‑101 (MQ=255)
           gTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAcc                                                                                      >  1:3296970/1‑101 (MQ=255)
              cGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACctct                                                                                   >  1:6608339/1‑101 (MQ=255)
               gCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtt                                                                                  >  1:2851346/1‑101 (MQ=255)
                cTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCttt                                                                                 >  1:5263586/1‑101 (MQ=255)
                 tggtggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCCCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtttt                                                                                <  2:7053542/101‑1 (MQ=255)
                   gtggAAGCCGCATTTAGCCAATTAAGTGCTGACCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTa                                                                              <  1:6287938/101‑1 (MQ=255)
                   gtggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTa                                                                              >  1:2570372/1‑101 (MQ=255)
                    tggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTAt                                                                             <  1:5959610/101‑1 (MQ=255)
                    tggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGGTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTAt                                                                             <  2:380942/101‑1 (MQ=255)
                     ggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATc                                                                            <  2:944483/101‑1 (MQ=255)
                     ggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATc                                                                            <  2:4757887/101‑1 (MQ=255)
                       aaGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCgg                                                                          <  2:2396286/101‑1 (MQ=255)
                        cgCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGc                                                                         <  2:4521343/100‑1 (MQ=255)
                        aGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGc                                                                         <  2:1564277/101‑1 (MQ=255)
                        aGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGc                                                                         >  2:4989176/1‑101 (MQ=255)
                           cGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAtt                                                                      <  2:1325101/101‑1 (MQ=255)
                            gCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAttt                                                                     <  2:5563661/101‑1 (MQ=255)
                              aTTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTcc                                                                   <  2:3141660/101‑1 (MQ=255)
                              aTTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTcc                                                                   >  2:2917516/1‑101 (MQ=255)
                               tttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                                                                  >  2:1771603/1‑101 (MQ=255)
                               tttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                                                                  <  2:1736784/101‑1 (MQ=255)
                               tttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                                                                  <  1:6516058/101‑1 (MQ=255)
                                 tAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcgcg                                                                <  1:5406165/101‑1 (MQ=255)
                                  aGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGa                                                               >  2:1269642/1‑101 (MQ=255)
                                     cAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt                                                            >  2:2168199/1‑101 (MQ=255)
                                       aTTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAg                                                          <  1:253087/101‑1 (MQ=255)
                                        ttAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGa                                                         <  1:2951328/101‑1 (MQ=255)
                                         tAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGAc                                                        >  1:1497636/1‑101 (MQ=255)
                                         tAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGAc                                                        <  1:363736/101‑1 (MQ=255)
                                             tGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAAc                                                    <  1:4660590/101‑1 (MQ=255)
                                             tGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAAc                                                    <  2:1344967/101‑1 (MQ=255)
                                              gCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTTTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                                   <  2:1052979/101‑1 (MQ=255)
                                              gCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAAcc                                                   >  2:3441719/1‑100 (MQ=255)
                                              gCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                                   <  2:6094848/101‑1 (MQ=255)
                                              gCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                                   <  1:1716613/101‑1 (MQ=255)
                                               cTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg                                                  <  1:441633/101‑1 (MQ=255)
                                                tGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGCAGCGGTCATCGCTCTCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACAgg                                                 >  2:6710471/1‑101 (MQ=255)
                                                tGAACGCAGCTTAACCAGGCTGAGTTTGCGTGAAGTGGCGCGTGGATCGGTCATTGAGCCCACCTCTTTTTATCGGAATTTCCGCGACGTATACTAACTgg                                                 <  2:6532467/101‑1 (MQ=255)
                                                 gAACGCAGCTCCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTggg                                                <  2:2752923/101‑1 (MQ=255)
                                                  aaCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGt                                               >  2:3065899/1‑101 (MQ=255)
                                                   aCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTc                                              <  1:805068/101‑1 (MQ=255)
                                                    cGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCCCCTCTTTTTATCGGCATTCCCGCGACGTAGACGAACTGGGCCt                                             <  2:3757959/101‑1 (MQ=255)
                                                    cGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCt                                             <  1:1669667/101‑1 (MQ=255)
                                                    cGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCt                                             <  1:6134852/101‑1 (MQ=255)
                                                      cAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGACTAGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGa                                           >  2:3340720/1‑101 (MQ=255)
                                                      cAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGa                                           <  1:224323/101‑1 (MQ=255)
                                                       aGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCCTTTCCGCGACGTAGACGAACTGGGTCTGAc                                          >  1:4868151/1‑101 (MQ=255)
                                                         cTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCCCCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCa                                        <  2:3945189/101‑1 (MQ=255)
                                                          ttCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCAt                                       <  1:3115355/101‑1 (MQ=255)
                                                          ttCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCAt                                       <  1:1994646/101‑1 (MQ=255)
                                                              ccAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGtt                                   >  2:4366575/1‑101 (MQ=255)
                                                              ccAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGtt                                   <  1:4999894/101‑1 (MQ=255)
                                                               cAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtg                                  >  2:6985683/1‑101 (MQ=255)
                                                                aGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtga                                 >  1:6783561/1‑101 (MQ=255)
                                                                  ccTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatg                               >  1:3524552/1‑101 (MQ=255)
                                                                     gAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgaga                            >  2:223846/1‑101 (MQ=255)
                                                                     gAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgaga                            <  1:4014995/101‑1 (MQ=255)
                                                                      aGTTTGCGTGGAGTGGCGCGTGAAGCGGTCATTGCTCCCCCCTCCTTTTATCGGCATTTCCGCTACGTCGGCGAACTGGGTCTTCCCATGGTTGATGagag                           <  2:4616557/101‑1 (MQ=255)
                                                                      aGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGagag                           >  1:1393290/1‑101 (MQ=255)
                                                                      aGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGagag                           >  1:1023411/1‑101 (MQ=255)
                                                                        tttGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCg                         >  1:902997/1‑101 (MQ=255)
                                                                          tGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGt                       <  1:4550005/101‑1 (MQ=255)
                                                                            agTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGGAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGttt                     <  2:3296970/100‑1 (MQ=255)
                                                                             gTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTa                    >  1:5763661/1‑101 (MQ=255)
                                                                               gAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAAt                  >  2:4817910/1‑101 (MQ=255)
                                                                                aaGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATg                 <  1:4317228/101‑1 (MQ=255)
                                                                                aaGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATg                 <  2:6608339/101‑1 (MQ=255)
                                                                                   tGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTa              >  1:2192960/1‑101 (MQ=255)
                                                                                   tGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTa              <  1:4078053/101‑1 (MQ=255)
                                                                                   tGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTa              <  2:179526/101‑1 (MQ=255)
                                                                                   tGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTa              <  1:2667759/101‑1 (MQ=255)
                                                                                       gcgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc          <  2:4143580/101‑1 (MQ=255)
                                                                                       gcgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc          <  2:2570372/101‑1 (MQ=255)
                                                                                        cgTAAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCa         <  1:2889408/101‑1 (MQ=255)
                                                                                          tGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAAc       >  1:646905/1‑101 (MQ=255)
                                                                                          tGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAAc       >  2:5232076/1‑101 (MQ=255)
                                                                                           gAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACt      >  2:3186926/1‑101 (MQ=255)
                                                                                            aaGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTc     >  2:6208169/1‑101 (MQ=255)
                                                                                            aaGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTc     <  2:5815832/101‑1 (MQ=255)
                                                                                             aGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCa    <  1:1226997/101‑1 (MQ=255)
                                                                                              gCGGTCATTGCACCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCAt   <  2:804853/101‑1 (MQ=255)
                                                                                               cGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATg  <  2:2851346/101‑1 (MQ=255)
                                                                                               cGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATg  >  2:376513/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
AAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATG  >  NC_000913/4159173‑4159368

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