Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,980,746 C→T 16.3% intergenic (‑111/+18) 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,7460CT16.3% 362.3 / 50.9 151intergenic (‑111/+18)ygeA/araEAsp/Glu_racemase family protein/arabinose transporter
Reads supporting (aligned to +/- strand):  ref base C (69/57);  new base T (13/12);  total (82/69)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 156 total aligned reads displayed.

GATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGA  >  NC_000913/2980647‑2980845
                                                                                                   |                                                                                                   
gATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGc                                                                                                     <  2:427466/100‑1 (MQ=255)
gATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAAAGAAAAAAGGGCCGGATGTACAGc                                                                                                     <  1:1867282/100‑1 (MQ=255)
      tgtTACGGCAACAAGGCTGAATAAGGAACGAGACAAGCGCCTAAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATcc                                                                                               <  2:3606961/100‑1 (MQ=255)
      tgtTACCGCAACAAGGCTGAATAAGGACCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGCATGTACAGCACATcc                                                                                               <  2:1741351/100‑1 (MQ=255)
      tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                               >  2:1442003/1‑100 (MQ=255)
      tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                               <  1:2838730/100‑1 (MQ=255)
       gtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCg                                                                                              >  1:2365503/1‑100 (MQ=255)
       gtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCg                                                                                              <  2:7558/100‑1 (MQ=255)
         tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCCACGGTCAAGTGCACAATCGCTATTAACGAAAAAAGGTCCGGATGTACAGcccacac                                                                                               >  2:83395/1‑94 (MQ=255)
         tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                            <  1:1497610/100‑1 (MQ=255)
         tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                            <  2:429982/100‑1 (MQ=255)
         tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGc                                                                                            >  2:1283690/1‑100 (MQ=255)
         tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGc                                                                                            <  2:101225/100‑1 (MQ=255)
          acaagaacaaGACTGAATAAGGAACGAGACGACCGCCTCAACGGCCAAGTGCCCAATAAAAAAAAAAAAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                            <  2:3260335/94‑1 (MQ=12)
          aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                           <  1:802264/100‑1 (MQ=39)
          aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGcc                                                                                           <  1:111261/100‑1 (MQ=255)
           ccGCAAGAAGTCTCAATAAGGAACTAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGccc                                                                                          <  1:1398805/100‑1 (MQ=37)
           ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGccc                                                                                          <  2:846077/100‑1 (MQ=39)
            cGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCg                                                                                         >  1:3288815/1‑100 (MQ=255)
             gCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGt                                                                                        <  1:318243/100‑1 (MQ=255)
              caacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTg                                                                                       >  2:1246756/1‑100 (MQ=255)
               aacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGa                                                                                      <  2:1072895/100‑1 (MQ=255)
               aacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGa                                                                                      <  1:3546512/100‑1 (MQ=255)
                 caaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCTCGTGaaa                                                                                    <  2:3535508/100‑1 (MQ=255)
                  aaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGcccttttttc                                                                                   <  2:1527761/100‑8 (MQ=37)
                   aGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATc                                                                                  >  2:2642441/1‑100 (MQ=255)
                    ggCTGAATAAGGAACGAGACAAACGACTCAACGGCCAAGTGCACAATCTCTATTAACGAAAAAATGGCCGGATGTACAGCACATCCGGCCCGTGAAATc                                                                                  <  2:2570969/99‑1 (MQ=255)
                      cTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGa                                                                               >  1:989513/1‑100 (MQ=255)
                        gAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACg                                                                             <  1:2529722/100‑1 (MQ=255)
                        gAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACg                                                                             >  2:1806807/1‑100 (MQ=255)
                             aGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGat                                                                        >  1:380397/1‑100 (MQ=255)
                              ggAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGata                                                                       >  2:799102/1‑100 (MQ=255)
                              ggAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGata                                                                       >  1:3483076/1‑100 (MQ=255)
                               gAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGatat                                                                      >  2:3755856/1‑100 (MQ=255)
                               gAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGatat                                                                      >  2:2305809/1‑100 (MQ=255)
                               gAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGatat                                                                      >  1:3546132/1‑100 (MQ=255)
                                 aCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATAttt                                                                    >  1:3320787/1‑100 (MQ=255)
                                  cGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTATATAAAAAAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                                   <  2:1092530/100‑1 (MQ=255)
                                     gaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCa                                                                >  1:2249827/1‑100 (MQ=255)
                                      aCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCaa                                                               <  1:2435189/100‑1 (MQ=255)
                                       cgaaCGCCTCAACGGCCAAGTGCCCAATAAAAAAAAAAAAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAAc                                                              <  2:590018/98‑1 (MQ=255)
                                          aCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTtc                                                           <  2:2102489/100‑1 (MQ=255)
                                             ccTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                        >  1:2915006/1‑100 (MQ=255)
                                             ccTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                        >  1:1417097/1‑100 (MQ=255)
                                              cTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGc                                                       >  2:2913450/1‑100 (MQ=255)
                                               tCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGcc                                                      >  2:3116651/1‑100 (MQ=255)
                                               tCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGcc                                                      <  2:1794675/100‑1 (MQ=255)
                                                cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCt                                                     >  2:2255428/1‑100 (MQ=255)
                                                cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCt                                                     >  1:1102038/1‑100 (MQ=255)
                                                 aaCGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTg                                                    <  1:628015/100‑1 (MQ=255)
                                                 aaCGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTg                                                    <  1:2059574/100‑1 (MQ=255)
                                                  aCGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGc                                                   <  2:2850275/100‑1 (MQ=255)
                                                         aGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGt                                            <  2:2295581/100‑1 (MQ=255)
                                                          gTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGtt                                           <  1:3352523/100‑1 (MQ=255)
                                                           tGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGttt                                          <  2:1466379/100‑1 (MQ=255)
                                                           tGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGttt                                          >  2:948562/1‑100 (MQ=255)
                                                             cccAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGc                                        <  2:3528603/100‑1 (MQ=255)
                                                               cAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGt                                      <  1:1091390/100‑1 (MQ=255)
                                                                 ntctcTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTc                                    >  2:886153/2‑100 (MQ=255)
                                                                    tctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGat                                 >  1:1881162/1‑100 (MQ=255)
                                                                    agtgttAACGAAAAAAGGGCCGGATGGACAGCACATGCGGCCCGTGAAATCAGACGCGCATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGat                                 <  2:1584416/96‑1 (MQ=255)
                                                                     ctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGatn                                <  2:282848/100‑2 (MQ=255)
                                                                        ttAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGt                             <  1:1733518/100‑1 (MQ=255)
                                                                          aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                           >  2:3820013/1‑100 (MQ=255)
                                                                          aaCGAAAAAAGGGCCGGATGTACAGCACATCCAGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                           <  2:679436/100‑1 (MQ=255)
                                                                            cGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCa                         >  1:466040/1‑100 (MQ=255)
                                                                             gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                        >  1:2813341/1‑100 (MQ=255)
                                                                              aaaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCATCTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCa                         <  2:2351527/98‑1 (MQ=255)
                                                                               aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGatagatcggaaga                      >  1:1043862/1‑90 (MQ=255)
                                                                                aaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGt                     <  1:3755884/100‑1 (MQ=255)
                                                                                 aaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTg                    <  1:1524667/100‑1 (MQ=255)
                                                                                  aaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTg                    <  2:3462237/99‑1 (MQ=255)
                                                                                   cgggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGAc                  >  1:1623276/2‑100 (MQ=38)
                                                                                   cgggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGAc                  >  1:3076718/2‑100 (MQ=38)
                                                                                    gggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa                 >  1:2035417/1‑100 (MQ=39)
                                                                                    gggCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa                 >  2:1346001/1‑100 (MQ=255)
                                                                                     ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt                >  1:1550367/1‑100 (MQ=39)
                                                                                     ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt                >  2:1686870/1‑100 (MQ=39)
                                                                                     ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt                >  1:563258/1‑100 (MQ=39)
                                                                                     ggCCGGATGGGCTGTACATGCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa                 >  2:3129074/1‑99 (MQ=255)
                                                                                      gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt               >  1:2919537/1‑100 (MQ=255)
                                                                                      gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt               >  1:2024021/1‑100 (MQ=255)
                                                                                      gCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTCTGCGTTCGATATGTTCCAGCGTGACAtt               >  2:3222412/1‑100 (MQ=255)
                                                                                      gCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCACAAATTCTCGCCTGCCACCCGTCTGCGTTCGATATGTTCCAGCGTGAc                  >  2:2696726/1‑97 (MQ=255)
                                                                                        cGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtttt             >  1:2696486/1‑100 (MQ=255)
                                                                                         ggATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAttttt            >  1:1349868/1‑100 (MQ=255)
                                                                                         ggATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAttttt            >  1:3474887/1‑100 (MQ=255)
                                                                                         ggATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAttttt            <  1:793633/100‑1 (MQ=255)
                                                                                          gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCTACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTg           <  1:2960863/100‑1 (MQ=255)
                                                                                          gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTg           <  1:1761447/100‑1 (MQ=255)
                                                                                          gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTg           <  2:2864209/100‑1 (MQ=255)
                                                                                          gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTg           <  1:2987722/100‑1 (MQ=255)
                                                                                           aTGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTgg          <  2:2903317/100‑1 (MQ=255)
                                                                                             gTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGtt        >  2:1015809/1‑100 (MQ=255)
                                                                                               aCAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTc      >  2:3638502/1‑100 (MQ=255)
                                                                                               aCAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTc      >  2:240778/1‑100 (MQ=255)
                                                                                                cAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTcc     >  1:3301394/1‑100 (MQ=255)
                                                                                                cAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTcc     >  1:121032/1‑100 (MQ=255)
                                                                                                 aGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCg    <  1:3303899/100‑1 (MQ=255)
                                                                                                  gCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCgg   >  1:1202400/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGA  >  NC_000913/2980647‑2980845

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: