Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,980,742 A→G 18.6% intergenic (‑107/+22) ygeA ← / ← araE Asp/Glu_racemase family protein/arabinose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,980,7420AG18.6% 325.6 / 67.6 150intergenic (‑107/+22)ygeA/araEAsp/Glu_racemase family protein/arabinose transporter
Reads supporting (aligned to +/- strand):  ref base A (80/42);  new base G (19/9);  total (99/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

Only 100 of 153 total aligned reads displayed.

CCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTC  >  NC_000913/2980644‑2980841
                                                                                                  |                                                                                                   
 cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAATGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                   <  1:2562629/100‑1 (MQ=255)
   gATGTGTGTTACCGCAACAAGGCTGAATAAGGATCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGc                                                                                                 <  2:925036/100‑1 (MQ=255)
   gATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGc                                                                                                 >  1:71735/1‑100 (MQ=255)
   gATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTAAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGc                                                                                                 <  2:3464940/100‑1 (MQ=255)
    atgtatgtTACCGCAACAAGGCTGAATAAGGAACGAAACAAAAGACTCAACGGCCAAGTGCCCAATAAAAAAAAAAAAAAAAAGGGCCGGATGTACAGca                                                                                                <  2:798474/100‑1 (MQ=255)
    aTGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGca                                                                                                <  2:284550/100‑1 (MQ=255)
     tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcac                                                                                               >  1:3447010/1‑100 (MQ=255)
     tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcac                                                                                               >  1:3560040/1‑100 (MQ=255)
      gtgtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcaca                                                                                              >  2:3552138/1‑100 (MQ=255)
      gtgtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcaca                                                                                              >  1:1467101/1‑100 (MQ=255)
         tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                           >  2:3507072/1‑100 (MQ=255)
         tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                           >  1:265130/1‑100 (MQ=255)
         tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                           >  2:1592203/1‑100 (MQ=255)
          gtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCg                                                                                          >  1:2270524/1‑100 (MQ=255)
          gtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCg                                                                                          >  1:1675169/1‑100 (MQ=255)
           ttACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCgg                                                                                         <  1:1834403/100‑1 (MQ=255)
            tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                        <  1:1625137/100‑1 (MQ=255)
            tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                        <  1:2893887/100‑1 (MQ=255)
            tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGc                                                                                        >  1:1288892/1‑100 (MQ=255)
             aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                        <  2:57299/99‑1 (MQ=39)
             aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGcc                                                                                       >  1:1224264/1‑100 (MQ=255)
             aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTAGTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                       <  2:1950463/100‑1 (MQ=38)
              ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGccc                                                                                      >  1:552150/1‑100 (MQ=255)
                gCAACAAGCCTGAATAAGGAACAAGAAAAACGCCTCAACGGCCAAGTGCCCAATAAAAAAAAAAAAAAAAAGGGCCGGATGTACAGCACATCCGGCCCg                                                                                     <  2:1772126/99‑1 (MQ=255)
                 caacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTg                                                                                   >  1:1593777/1‑100 (MQ=255)
                 caacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTg                                                                                   <  1:3054535/100‑1 (MQ=255)
                  aacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATATCAATAAAAAAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGa                                                                                  <  2:2660402/100‑1 (MQ=255)
                      aGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGCAAAAAGGACCGGATGTACAGCACATCCGGCCCGTGAAATc                                                                              >  1:41396/1‑100 (MQ=255)
                        gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAg                                                                            <  1:47183/100‑1 (MQ=255)
                         cTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGa                                                                           <  1:2189516/100‑1 (MQ=255)
                           gAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACg                                                                         >  1:2287976/1‑100 (MQ=255)
                           gAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACg                                                                         >  2:3145792/1‑100 (MQ=255)
                             tctaGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGcc                                                                       <  2:693787/97‑1 (MQ=255)
                                aGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGccaga                                                                    >  1:693790/1‑97 (MQ=255)
                                 ggAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGata                                                                   >  1:3135024/1‑100 (MQ=255)
                                  gAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGatat                                                                  <  1:3398653/100‑1 (MQ=255)
                                   aaCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATAtt                                                                 >  1:1853291/1‑100 (MQ=255)
                                    aCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATAttt                                                                >  1:1885681/1‑100 (MQ=255)
                                    aCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATAttt                                                                <  1:2424559/100‑1 (MQ=255)
                                      gagaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                               <  2:2802447/99‑1 (MQ=255)
                                       agaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtctc                                                             >  2:1070103/1‑100 (MQ=255)
                                        gaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCa                                                            <  1:1485679/100‑1 (MQ=255)
                                          cAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAAc                                                          >  1:3645411/1‑100 (MQ=255)
                                           aaaCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACt                                                         >  2:144766/1‑100 (MQ=255)
                                           aaaCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACt                                                         >  2:2881471/1‑100 (MQ=255)
                                                ccTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                    >  1:2674181/1‑100 (MQ=255)
                                                  tCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCCCGcc                                                  <  1:3566787/100‑1 (MQ=255)
                                                   cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCt                                                 >  1:2927396/1‑100 (MQ=255)
                                                   cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCt                                                 >  2:2156321/1‑100 (MQ=255)
                                                      cGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGcc                                              >  2:2736875/1‑100 (MQ=255)
                                                      cGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGcc                                              >  2:2972961/1‑100 (MQ=255)
                                                      cGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGACCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGcc                                              >  2:3371181/1‑100 (MQ=255)
                                                        gCCAAGTGCCTTATCCCTATTAGCAAAAAAATGGCCGGATGTATAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCa                                             <  2:2129647/99‑1 (MQ=255)
                                                              tGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGttt                                      >  2:3574553/1‑100 (MQ=255)
                                                               gCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTg                                     <  1:974741/100‑1 (MQ=255)
                                                               gCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTg                                     >  2:849703/1‑100 (MQ=255)
                                                                cccAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAACTCAGACGCCGATATTCCTCAACTTCTCGCCTGCCATCAGTTTGc                                    >  2:766029/1‑100 (MQ=255)
                                                                 ccAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAAACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCg                                   >  2:1047106/1‑100 (MQ=255)
                                                                  cAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGt                                  <  1:85999/100‑1 (MQ=255)
                                                                  cAATCTCTATTAACAAAAAAAGGGCCCGATATACAGCACATCCGGCCCGTGAAATCAGCCGCCGATATTTCTCAACTTCTCGCCTGCCATCATTTTGCg                                   <  1:81929/99‑1 (MQ=255)
                                                                    aTCACTATTAAAGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTc                                <  2:1068763/100‑1 (MQ=255)
                                                                       tctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGa                              <  2:1367083/99‑1 (MQ=255)
                                                                        ctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCCACTTCTCGCCTGCCATCAGTTTGCGTTCGata                            >  2:889529/1‑100 (MQ=255)
                                                                        ctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGata                            >  2:2361941/1‑100 (MQ=255)
                                                                            tAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGtt                        >  1:3236356/1‑100 (MQ=255)
                                                                            tAACGAAAAAAGCGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGtt                        <  2:41396/100‑1 (MQ=255)
                                                                             ctgatttcacgggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                       >  1:559739/11‑100 (MQ=25)
                                                                             aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                       >  2:3334455/1‑100 (MQ=255)
                                                                              aCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCATCTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTcc                      <  2:1946524/100‑1 (MQ=255)
                                                                               cGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCa                     >  1:1500643/1‑100 (MQ=255)
                                                                               cGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCa                     >  1:161166/1‑100 (MQ=255)
                                                                               cGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCa                     >  1:986953/1‑100 (MQ=255)
                                                                                gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                    >  1:1449170/1‑100 (MQ=255)
                                                                                gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                    <  1:2159953/100‑1 (MQ=255)
                                                                                gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                    <  2:3236331/100‑1 (MQ=255)
                                                                                  naaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                  >  2:1883813/2‑100 (MQ=255)
                                                                                  aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                  >  2:1860270/1‑100 (MQ=255)
                                                                                    aaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTg                <  2:167396/100‑1 (MQ=255)
                                                                                      aGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGAc              <  1:3645396/100‑1 (MQ=255)
                                                                                       gggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa             >  1:3257368/1‑100 (MQ=39)
                                                                                       gggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa             >  1:447974/1‑100 (MQ=39)
                                                                                       gggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTAGATATGTTCCAGCGTGCCa             >  1:346373/1‑100 (MQ=38)
                                                                                       gggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGATTGCGTTCGATATGTTCCAGCGTGACa             >  2:630605/1‑100 (MQ=38)
                                                                                       gggCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa             <  2:2674163/100‑1 (MQ=255)
                                                                                        ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt            >  1:1193100/1‑100 (MQ=39)
                                                                                        ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt            >  1:2863087/1‑100 (MQ=39)
                                                                                        ggCCGGATCTGCTGTACATCCGGCCCGTGCAATCAGACCCCGAAATTTCTCAACTTCTCGCCTGCCACCAACTTGCGTTCGCTATGTTCCAGCGTGACa             >  2:3510538/1‑99 (MQ=255)
                                                                                         gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt           >  1:39578/1‑100 (MQ=255)
                                                                                         gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt           >  2:1054898/1‑100 (MQ=255)
                                                                                         gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt           >  2:1500989/1‑100 (MQ=255)
                                                                                         gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt            >  2:1305714/1‑99 (MQ=39)
                                                                                         gCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt           >  2:3463992/1‑100 (MQ=255)
                                                                                         gCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt           >  1:110544/1‑100 (MQ=255)
                                                                                           cGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtttt         >  2:1818722/1‑100 (MQ=255)
                                                                                             gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTg       <  1:1592211/100‑1 (MQ=255)
                                                                                               tGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGt     >  2:2488955/1‑100 (MQ=255)
                                                                                                 tACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGttt   >  1:847996/1‑100 (MQ=255)
                                                                                                 tACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGttt   >  2:1919682/1‑100 (MQ=255)
                                                                                                  aCAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTc  >  1:315948/1‑100 (MQ=255)
                                                                                                  aCAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTc  >  1:2056401/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTC  >  NC_000913/2980644‑2980841

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: