Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,095 T→A 20.1% intergenic (+36/+20) purT → / ← eda phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,932,0950TA20.1% 255.3 / 55.4 114intergenic (+36/+20)purT/edaphosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase
Reads supporting (aligned to +/- strand):  ref base T (44/47);  new base A (11/12);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 116 total aligned reads displayed.

GTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCG  >  NC_000913/1931998‑1932194
                                                                                                 |                                                                                                   
gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCg                                                                                                   >  2:2008422/1‑100 (MQ=255)
  gttgttGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGAt                                                                                                 >  1:1999748/1‑100 (MQ=255)
  gttgttGACGCCATTGAACGCGCGAAACACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGAt                                                                                                 >  2:1774996/1‑100 (MQ=255)
    tgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCg                                                                                               >  1:1988292/1‑100 (MQ=255)
     gttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCgg                                                                                              <  2:2152387/100‑1 (MQ=255)
     gttgACGCCAGTGAACGCGCGAAGCACGCGGCCAGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCgg                                                                                              <  1:1496138/100‑1 (MQ=255)
      ttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCggg                                                                                             <  1:1329254/100‑1 (MQ=255)
       tgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGc                                                                                            >  2:2574693/1‑100 (MQ=255)
       tgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGc                                                                                            <  2:1049560/100‑1 (MQ=255)
        gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                           <  2:2236136/100‑1 (MQ=255)
        gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                                                           >  2:1189849/1‑100 (MQ=255)
        gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGGAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGc                                                                                            <  1:387997/99‑1 (MQ=255)
         aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAt                                                                                          <  2:2333794/100‑1 (MQ=39)
         aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAt                                                                                          <  2:535563/100‑1 (MQ=39)
         aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAt                                                                                          <  2:2028599/100‑1 (MQ=39)
         aCGCCATTGAACGCGCGAAGCACGACGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAt                                                                                          <  1:2763371/100‑1 (MQ=39)
          cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtt                                                                                         <  2:2812234/100‑1 (MQ=39)
          cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtt                                                                                         <  2:1263242/100‑1 (MQ=39)
          cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtt                                                                                         <  2:1926132/100‑1 (MQ=255)
           gctaTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAttt                                                                                        >  2:408116/4‑100 (MQ=255)
            cctatGAACTCGAGAAGCACGCCGCTGGATAGTTAAAAGAACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtttt                                                                                       <  1:1437003/96‑1 (MQ=25)
            ccATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                                                       >  1:127672/1‑100 (MQ=255)
            ccATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                                                       >  1:2725383/1‑100 (MQ=255)
              aTTGAACGCCCGAAGCACGCCAGCAGAGAGGGATAATTACAGGGTTAAACTCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAttt                                                                                        <  1:2666390/97‑1 (MQ=21)
                    cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                               >  2:648433/1‑100 (MQ=255)
                     gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                                                              >  2:1214084/1‑100 (MQ=255)
                     gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTGAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                                                              <  2:2415146/100‑1 (MQ=255)
                       gcgAAGCACGCCACCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACa                                                                             >  2:1916837/1‑99 (MQ=255)
                        cgcaGCACGCCGCCAGACAGGTAAAAGTACAGAATTAAACCCGTCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                                                              <  1:1602100/94‑1 (MQ=255)
                             cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                      >  1:1567430/1‑100 (MQ=255)
                             cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAg                                                                       >  2:649073/1‑99 (MQ=255)
                              aCGCCGCCGGACAGGTAAAAGTACAGGGTTAACCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAg                                                                       <  1:1722104/98‑1 (MQ=255)
                              aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACTCGGCAAAAAAACGCTACAAAAATGCAAGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                      >  1:552603/1‑99 (MQ=255)
                              aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACTCGGCAAAAAAACGCTACAAAAATGACCGATCTTCGATCGGGCATTTTGACTTTTACAGCTTAgcg                                                                     >  2:1762350/1‑100 (MQ=255)
                              aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                      <  1:1325799/99‑1 (MQ=255)
                               cgccgcCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcgc                                                                    >  1:883906/1‑100 (MQ=255)
                                gccgccGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGcc                                                                   <  2:116537/100‑1 (MQ=255)
                                 ccgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCt                                                                  >  2:2362770/1‑100 (MQ=255)
                                 ccgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCt                                                                  <  2:2269088/100‑1 (MQ=255)
                                 ccgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCt                                                                  >  2:2508198/1‑100 (MQ=255)
                                     cgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTc                                                                <  1:770555/98‑1 (MQ=255)
                                      ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                             <  2:695334/100‑1 (MQ=255)
                                      ggACAGGTAAAAGTACAGGGTAAAAAAAAAAAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTa                                                              <  1:2112186/99‑1 (MQ=255)
                                       gACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                            <  2:654670/100‑1 (MQ=255)
                                          aGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCt                                                         <  1:2192374/100‑1 (MQ=255)
                                           ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                        >  2:1703281/1‑100 (MQ=255)
                                             tAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCa                                                      >  1:2782539/1‑100 (MQ=255)
                                             tAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCa                                                      >  2:1185179/1‑100 (MQ=255)
                                             tAAAAGTACAGGGTTAAACCCGGCAAAAAAAAGCTACAAAAATGCATGATTATCGATCGGGCATTTTGACATTTACAGCTTATCGCCTTCTACAGCtt                                                        >  1:1822335/1‑98 (MQ=255)
                                                 aGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgcg                                                  >  2:312366/1‑100 (MQ=255)
                                                         gTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                          >  2:124283/1‑100 (MQ=255)
                                                         gTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                          >  2:2589928/1‑100 (MQ=255)
                                                             aaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTaa                                      <  1:2123615/100‑1 (MQ=255)
                                                             aaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTaa                                      >  2:880842/1‑100 (MQ=255)
                                                               cccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTa                                       <  2:2472175/97‑1 (MQ=255)
                                                               cccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATg                                    >  2:718114/1‑100 (MQ=255)
                                                                ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGc                                   >  2:2279403/1‑100 (MQ=255)
                                                                ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGc                                   <  1:1507771/100‑1 (MQ=255)
                                                                   gCAAAACAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGt                                >  1:2907307/1‑100 (MQ=255)
                                                                   gCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGt                                <  2:2478756/100‑1 (MQ=255)
                                                                    cAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTc                               >  2:1259162/1‑100 (MQ=255)
                                                                     aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCg                              <  2:392837/100‑1 (MQ=255)
                                                                      aaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGt                             >  2:194549/1‑100 (MQ=255)
                                                                      aaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGt                             <  2:883909/100‑1 (MQ=255)
                                                                        aaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTaa                           <  2:2510537/100‑1 (MQ=255)
                                                                         tgtaaaagtcaaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAAt                          >  2:1442216/11‑100 (MQ=25)
                                                                         aaaCGTTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAAt                          <  2:1289407/100‑1 (MQ=255)
                                                                             gCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGcc                      >  1:544371/1‑100 (MQ=255)
                                                                             gCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGcc                      <  1:2792142/100‑1 (MQ=255)
                                                                              cTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGccc                     <  2:843238/100‑1 (MQ=255)
                                                                               tACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                    >  2:822959/1‑100 (MQ=255)
                                                                               tACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                    <  2:76510/100‑1 (MQ=255)
                                                                               agtcaaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                    >  1:1327298/5‑100 (MQ=38)
                                                                                gtcaaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTTATCGCCCGc                   >  1:1605754/4‑100 (MQ=37)
                                                                                 tcaaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCt                  >  2:2759154/3‑100 (MQ=38)
                                                                                  caaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCtt                 >  2:2008595/2‑100 (MQ=38)
                                                                                  caaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTTAAGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCtt                 >  1:1684228/2‑100 (MQ=255)
                                                                                   aaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                >  1:2101698/1‑100 (MQ=38)
                                                                                    aaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCTCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTcc               >  2:1051649/1‑100 (MQ=38)
                                                                                    aaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTcc               <  2:1120340/100‑1 (MQ=255)
                                                                                     aaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa              >  1:107550/1‑100 (MQ=255)
                                                                                      aTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa              >  1:509621/1‑99 (MQ=255)
                                                                                      aTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg             >  2:2511247/1‑100 (MQ=39)
                                                                                      aTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTGCAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg             <  2:1736348/100‑1 (MQ=255)
                                                                                      aTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg             <  2:2782557/100‑1 (MQ=255)
                                                                                      aTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg             >  1:1748680/1‑100 (MQ=255)
                                                                                       tGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc            <  2:703986/100‑1 (MQ=255)
                                                                                       tGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc            <  2:520728/100‑1 (MQ=255)
                                                                                       tGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc            <  2:1734629/100‑1 (MQ=255)
                                                                                        gCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgcg           <  1:1627644/100‑1 (MQ=255)
                                                                                         cccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgcgc          >  1:484027/1‑100 (MQ=255)
                                                                                          ccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa         <  2:931537/100‑1 (MQ=255)
                                                                                            gATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATc       >  2:2051739/1‑100 (MQ=255)
                                                                                              tCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTg     <  2:2005447/100‑1 (MQ=255)
                                                                                               ccTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTGTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTg     <  1:1927936/99‑1 (MQ=255)
                                                                                               ccTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGc    <  1:1185174/100‑1 (MQ=255)
                                                                                                cTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc   >  2:318729/1‑100 (MQ=255)
                                                                                                cTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc   >  1:1190456/1‑100 (MQ=255)
                                                                                                 tCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCg  <  2:2725401/100‑1 (MQ=255)
                                                                                                 tCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCg  >  2:2634251/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCG  >  NC_000913/1931998‑1932194

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: