Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,980,742 A→G 15.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,7420AG15.6% 351.3 / 52.6 154intergenic (‑107/+22)ygeA/araEAsp/Glu_racemase family protein/arabinose transporter
Reads supporting (aligned to +/- strand):  ref base A (66/63);  new base G (13/11);  total (80/74)
Fisher's exact test for biased strand distribution p-value = 8.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

Only 100 of 159 total aligned reads displayed.

TCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTC  >  NC_000913/2980643‑2980841
                                                                                                   |                                                                                                   
tcCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTa                                                                                                     >  1:2245323/1‑100 (MQ=255)
tcCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTa                                                                                                     >  1:196487/1‑100 (MQ=255)
 ccTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                    >  2:1542563/1‑100 (MQ=255)
 ccTGATGTGTGTGACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                    <  2:361763/100‑1 (MQ=255)
  cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTaaa                                                                                                   <  1:2592393/100‑3 (MQ=255)
  cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                    <  2:1809004/99‑1 (MQ=255)
  cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                   >  1:2843303/1‑100 (MQ=255)
  cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                   <  1:588947/100‑1 (MQ=255)
   tGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAg                                                                                                  >  1:1048259/1‑100 (MQ=255)
    gATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAAAGAAAAAAGGGCCGGATGTACAGc                                                                                                 <  1:1867282/100‑1 (MQ=255)
     aTGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGca                                                                                                <  1:3300501/100‑1 (MQ=255)
     aTGTGTGTTACCGAAACAAGGCTGAATAAGGAACGAGAAAAACGACTCAACGGCCAAGTGCCCAATAAAAAAAAAAAAAAAAAGGGCCGGATGTACa                                                                                                   <  2:592790/97‑1 (MQ=255)
          tgtTACGGCAACAAGGCTGAATAAGGAACGAGACAAGCGCCTAAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATcc                                                                                           <  2:3606961/100‑1 (MQ=255)
          tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                           <  1:2838730/100‑1 (MQ=255)
           gtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCg                                                                                          >  1:1092885/1‑100 (MQ=255)
           gtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCg                                                                                          <  2:7558/100‑1 (MQ=255)
           gtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCg                                                                                          >  1:2365503/1‑100 (MQ=255)
             tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCCACGGTCAAGTGCACAATCGCTATTAACGAAAAAAGGTCCGGATGTACAGcccacac                                                                                           >  2:83395/1‑94 (MQ=255)
             tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                        <  1:964999/100‑1 (MQ=255)
             tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                        <  2:429982/100‑1 (MQ=255)
             tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGc                                                                                        <  2:101225/100‑1 (MQ=255)
               ccGCAAGAAGTCTCAATAAGGAACTAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGccc                                                                                      <  1:1398805/100‑1 (MQ=37)
               ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGccc                                                                                      <  2:846077/100‑1 (MQ=39)
               ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGccc                                                                                      >  2:1137239/1‑100 (MQ=255)
                cGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCg                                                                                     <  1:1463971/100‑1 (MQ=255)
                 gCAACAAGGCAGAATAAGGAACGAGACAAACGCCGCAGCGGCCAAGTGCCCAATCAAAAAAAAAGAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                       <  2:3444658/97‑1 (MQ=19)
                  caacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTg                                                                                   <  1:2745897/100‑1 (MQ=255)
                  caacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTg                                                                                   >  2:1246756/1‑100 (MQ=255)
                   aacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGa                                                                                  <  1:3546512/100‑1 (MQ=255)
                   aacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGa                                                                                  <  2:1072895/100‑1 (MQ=255)
                     caaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCTCGTGaaa                                                                                <  2:3535508/100‑1 (MQ=255)
                     caaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGaaa                                                                                >  1:2927439/1‑100 (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                                                                         >  2:1806807/1‑100 (MQ=255)
                                  ggAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGata                                                                   >  1:3483076/1‑100 (MQ=255)
                                  ggAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGata                                                                   >  2:799102/1‑100 (MQ=255)
                                   gAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGgta                                                                   >  2:649568/1‑99 (MQ=255)
                                   gAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGatat                                                                  >  2:3755856/1‑100 (MQ=255)
                                   gAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGatat                                                                  >  1:3546132/1‑100 (MQ=255)
                                   gAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGatat                                                                  >  2:2305809/1‑100 (MQ=255)
                                    aaCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATAtt                                                                 >  1:1530009/1‑100 (MQ=255)
                                       gagaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtct                                                              <  1:3022694/100‑1 (MQ=255)
                                          aCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCaa                                                           <  1:2435189/100‑1 (MQ=255)
                                             aaCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACtt                                                        >  2:3751680/1‑100 (MQ=255)
                                              aCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTtc                                                       <  2:2102489/100‑1 (MQ=255)
                                                 ccTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                    >  1:1417097/1‑100 (MQ=255)
                                                 ccTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                    >  1:2915006/1‑100 (MQ=255)
                                                  cacAACGGCCAAGTGCCAAATCTCTATCAACGAAAAAAGGGCCTGAAGTACAGCACAGCGCGCACGTGAAATTAGTTGCCGATATTTCTCaa                                                           <  2:3211038/90‑1 (MQ=255)
                                                   tCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGcc                                                  >  1:2110044/1‑100 (MQ=255)
                                                   tCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGcc                                                  >  2:3116651/1‑100 (MQ=255)
                                                    cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCt                                                 >  2:2775381/1‑100 (MQ=255)
                                                    cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCt                                                 >  1:1102038/1‑100 (MQ=255)
                                                    cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCt                                                 >  2:2255428/1‑100 (MQ=255)
                                                      aCGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGc                                               <  2:2850275/100‑1 (MQ=255)
                                                       nggCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGcc                                              >  2:3070628/2‑100 (MQ=255)
                                                       cGGCCAAGTGCCCAATCTCTATTAACGAACAAAAGGCTGGCTATAGAGTCCATCATGCCTGTGAGATCAGGCGACGATAGTTCTCAACTGCTCGCCTg                                                <  2:1757045/98‑1 (MQ=255)
                                                          ccAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATc                                           <  1:704742/100‑1 (MQ=255)
                                                             aGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGt                                        <  2:2295581/100‑1 (MQ=255)
                                                              gTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGtt                                       <  1:3352523/100‑1 (MQ=255)
                                                               tGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGttt                                      >  2:948562/1‑100 (MQ=255)
                                                               tGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGttt                                      <  2:1466379/100‑1 (MQ=255)
                                                                gCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTg                                     <  1:125892/100‑1 (MQ=255)
                                                                     ntctcTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTc                                >  2:886153/2‑100 (MQ=255)
                                                                        tctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGat                             <  1:3211890/100‑1 (MQ=255)
                                                                        tctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGat                             >  1:1881162/1‑100 (MQ=255)
                                                                         ctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGatn                            <  2:282848/100‑2 (MQ=255)
                                                                            ttAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGt                         <  1:1733518/100‑1 (MQ=255)
                                                                             tAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGtt                        >  1:3307460/1‑100 (MQ=255)
                                                                              aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                       >  1:2032672/1‑100 (MQ=255)
                                                                                cGAAAAACGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGCTCGATATGTTCCa                     >  1:3428827/1‑100 (MQ=255)
                                                                                cGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCa                     >  1:1150895/1‑100 (MQ=255)
                                                                                cGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCa                     >  1:466040/1‑100 (MQ=255)
                                                                                  aaaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCATCTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCa                     <  2:2351527/98‑1 (MQ=255)
                                                                                   aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGatagatcggaaga                  >  1:1043862/1‑90 (MQ=255)
                                                                                    aaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGt                 >  1:522237/1‑100 (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)
                                                                                       aGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGAc              <  2:40731/100‑1 (MQ=255)
                                                                                       aGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGAc              >  2:1463717/1‑100 (MQ=255)
                                                                                        gggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa             >  1:2035417/1‑100 (MQ=39)
                                                                                        gggCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa             >  2:1346001/1‑100 (MQ=255)
                                                                                         ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt            >  2:1686870/1‑100 (MQ=39)
                                                                                          gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt           >  2:434011/1‑100 (MQ=255)
                                                                                          gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt           >  2:1666899/1‑100 (MQ=255)
                                                                                          gCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTCTGCGTTCGATATGTTCCAGCGTGACAtt           >  2:3222412/1‑100 (MQ=255)
                                                                                            cGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtttt         >  1:2696486/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)
                                                                                               aTGTACAGCACATCCTGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTgg      <  1:3763261/100‑1 (MQ=255)
                                                                                               aTGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTgg      <  2:2903317/100‑1 (MQ=255)
                                                                                                  tACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTCTTTGttt   >  2:1056448/1‑100 (MQ=255)
                                                                                                   aCAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTc  >  2:240778/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTC  >  NC_000913/2980643‑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: