Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,094 C→G 17.3% intergenic (+35/+21) 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,0940CG17.3% 335.3 / 61.5 127intergenic (+35/+21)purT/edaphosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase
Reads supporting (aligned to +/- strand):  ref base C (59/46);  new base G (13/9);  total (72/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 133 total aligned reads displayed.

AGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCC  >  NC_000913/1931995‑1932193
                                                                                                   |                                                                                                    
   gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑g                                                                                                  <  1:2482844/100‑1 (MQ=255)
   gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑g                                                                                                  >  1:711719/1‑100 (MQ=255)
     gttgttGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GAt                                                                                                <  1:610862/100‑1 (MQ=255)
       tgtggACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTAAAAAAAAAAAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATc                                                                                               <  2:931282/99‑1 (MQ=255)
           gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAG‑GATCGGGCa                                                                                          <  1:1356667/100‑1 (MQ=255)
           gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCa                                                                                          >  1:540585/1‑100 (MQ=255)
           gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCa                                                                                          <  1:2440163/100‑1 (MQ=255)
            aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAG‑GATCGGGCAt                                                                                         <  1:2106981/100‑1 (MQ=39)
            aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GACCGGGCAt                                                                                         >  1:1025846/1‑100 (MQ=255)
             cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCAGATCggaaga                                                                                         >  2:1104053/1‑96 (MQ=255)
             cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCAtt                                                                                        >  1:3729754/1‑100 (MQ=255)
             cGCCATTGAACGAGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCAGATCggaaga                                                                                         >  2:1268324/1‑96 (MQ=255)
              gcgaTTGAACGTGGGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAGAAAACGCTACAAAAATGCCCGATCGAG‑GATCGGGCAtt                                                                                        <  2:2351835/96‑1 (MQ=37)
              gCCATTGAACGCGCGAAGCACGCGGCCGGACAGGAAAAAGTACAGGGTTAAACCCGGAAAAAAAAAAAAAAAAAAATGCCCGATCGAG‑GATCGGGCAtt                                                                                        <  2:1821630/99‑1 (MQ=21)
              gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAG‑GATCGGGCAt                                                                                         <  2:1312000/98‑1 (MQ=255)
               ccATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATAGGGCAtt                                                                                        >  2:238008/1‑98 (MQ=255)
                cATTGAACGCGGGAGGCACGCCGCCAGAAAGGGAAAAGTACAGGGGTAAACCCGGCAAAAAAACGCTCCAAAAATGCCCGATCGAG‑GATCGGGCAtttt                                                                                      <  2:1507620/99‑1 (MQ=25)
                  ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGAc                                                                                   <  2:3099232/100‑1 (MQ=255)
                   tGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACt                                                                                  >  1:1716299/1‑100 (MQ=255)
                     aaCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACttt                                                                                >  1:2808839/1‑100 (MQ=255)
                      aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACtttt                                                                               >  1:1791699/1‑100 (MQ=255)
                      aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACtttt                                                                               <  2:412938/100‑1 (MQ=255)
                       cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTa                                                                              >  1:1092469/1‑100 (MQ=255)
                       cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTa                                                                              >  2:2671955/1‑100 (MQ=255)
                       cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGACCGATCCTC‑GATCGGGCATTTTGACTTTTa                                                                              >  2:3090917/1‑100 (MQ=255)
                       cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCAACAAAAATGAACGATCCTC‑AAACGGGCAAATTGACTTTTa                                                                              >  1:1268553/1‑100 (MQ=255)
                         cgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACa                                                                            >  2:1686472/1‑100 (MQ=255)
                           cgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGc                                                                          >  1:902969/1‑100 (MQ=255)
                               gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAg                                                                      >  1:1822223/1‑100 (MQ=255)
                                cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAgc                                                                     >  2:176268/1‑100 (MQ=255)
                                 aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAgcg                                                                    <  1:2609344/100‑1 (MQ=255)
                                  cgccgcCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAgcgc                                                                   >  1:3504618/1‑100 (MQ=255)
                                      gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTc                                                               >  1:2138053/1‑100 (MQ=255)
                                      gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTc                                                               <  1:770509/100‑1 (MQ=255)
                                       ccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCt                                                              >  2:877514/1‑100 (MQ=255)
                                       ccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCt                                                              <  1:3508681/100‑1 (MQ=255)
                                         ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                            <  2:1437915/100‑1 (MQ=255)
                                         ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                            <  1:2186225/100‑1 (MQ=255)
                                          gACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                           <  2:1770422/100‑1 (MQ=255)
                                              ggTAAAAGTACAGGGTTAAACCCGGCAAAACAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACCGCt                                                        >  2:1194535/1‑99 (MQ=255)
                                              ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                       >  1:3354073/1‑100 (MQ=255)
                                              ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                       >  1:3355377/1‑100 (MQ=255)
                                              ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                       >  1:1553190/1‑100 (MQ=255)
                                               gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                      >  1:3444797/1‑100 (MQ=255)
                                                  aaaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcg                                                   >  1:1383968/1‑100 (MQ=255)
                                                   aaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgc                                                  >  2:2626650/1‑100 (MQ=255)
                                                   aaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgc                                                  >  2:1276496/1‑100 (MQ=255)
                                                   aaGTACAGGGTTAAACCCGGCAAAAAAACGCCACAAAAATGCCCGATCCTC‑GATCGGGCATTTCCACTTTTACAGCTTAGCGCCTTCTACAGCTTCCcgc                                                  >  1:2450091/1‑100 (MQ=255)
                                                    aGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgcg                                                 <  1:1950937/100‑1 (MQ=255)
                                                     gTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTAACTTTTACAGCTTACCGCCTTCTACAGCTTCCCgcgc                                                >  2:3237767/1‑100 (MQ=255)
                                                      tACAGGGCTAAACCCGGAAAAAAAACGCTACAAAAGTGACTGAAACTC‑GAGCCGGCATTTTGACTTTGACAGCTTAGCGCCTTCTCAAGCTGCAtgcg                                                 <  2:2093997/98‑1 (MQ=255)
                                                       aCAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCa                                              >  2:619015/1‑100 (MQ=255)
                                                        cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                             >  1:1544578/1‑100 (MQ=255)
                                                          gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGc                                            <  1:1722579/99‑1 (MQ=255)
                                                          gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                           <  2:3159493/100‑1 (MQ=255)
                                                          gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                           >  1:1056819/1‑100 (MQ=255)
                                                          gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                           <  2:1943581/100‑1 (MQ=255)
                                                           ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                          <  2:1006724/100‑1 (MQ=255)
                                                           ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                          <  2:210948/100‑1 (MQ=255)
                                                            gTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                         <  1:2862951/100‑1 (MQ=255)
                                                             ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAg                                        <  1:3200360/100‑1 (MQ=255)
                                                              tAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGt                                       <  1:2932200/100‑1 (MQ=255)
                                                                aaCCCGGCAAAAAAACGCTACAAAAATGCCTGAGCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCACGCTTAGTaa                                     >  2:2503084/1‑100 (MQ=255)
                                                                 aCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAAt                                    >  1:3028809/1‑100 (MQ=255)
                                                                  nccGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATg                                   >  2:3117218/2‑100 (MQ=255)
                                                                  cccGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATg                                   >  1:1910454/1‑100 (MQ=255)
                                                                   ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGc                                  >  1:912196/1‑100 (MQ=255)
                                                                   ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGc                                  <  2:273879/100‑1 (MQ=255)
                                                                      gCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGt                               >  1:1028282/1‑100 (MQ=255)
                                                                        aaaaaaaCGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCg                             <  1:2795402/100‑1 (MQ=255)
                                                                         aaaaaaCGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGt                            <  1:360912/100‑1 (MQ=255)
                                                                          aaaaaCGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTa                           <  1:256800/100‑1 (MQ=255)
                                                                           aaaaCGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTaa                          <  1:2866643/100‑1 (MQ=255)
                                                                             aaCGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                        >  1:2069270/1‑100 (MQ=255)
                                                                             aaCGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                        <  2:3325140/100‑1 (MQ=255)
                                                                               cGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACTGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGc                      >  2:2582080/1‑100 (MQ=255)
                                                                                 aagtcaaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGcc                     >  2:2730022/6‑99 (MQ=37)
                                                                                  tACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                   >  2:3766280/1‑100 (MQ=255)
                                                                                   gtcaaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                  >  2:3321583/4‑100 (MQ=38)
                                                                                   aCAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                  >  1:1552301/1‑100 (MQ=255)
                                                                                    tcaaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCt                 >  1:2075618/3‑100 (MQ=38)
                                                                                    tcaaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCGTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCt                 >  1:171457/3‑100 (MQ=38)
                                                                                    cAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCCCCCGc                  >  2:3551194/1‑99 (MQ=255)
                                                                                     caaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCtt                >  1:2130442/2‑100 (MQ=38)
                                                                                      aaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc               >  2:348453/1‑100 (MQ=38)
                                                                                      aaaaTGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc               >  1:2791739/1‑100 (MQ=255)
                                                                                      aaaaTGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc               <  2:234358/100‑1 (MQ=255)
                                                                                        aaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTcc              >  1:748623/1‑99 (MQ=39)
                                                                                         aTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa             >  2:3500736/1‑99 (MQ=39)
                                                                                         aTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg            >  2:3115871/1‑100 (MQ=39)
                                                                                         aTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg            >  1:833512/1‑100 (MQ=39)
                                                                                         aTGCCCGATCGAG‑GAGCGGGCATTTTGACTTTTACAGCTTAGCTCGTTCCACAGCTTCAAGCGCCAGCGTAGTAATGCGGCCGTAATCGCCCGCGTCCa             >  2:3068943/1‑99 (MQ=255)
                                                                                          tGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc           <  1:572921/100‑1 (MQ=255)
                                                                                             ccGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa        <  1:3813731/100‑1 (MQ=255)
                                                                                             ccGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa        <  1:1965115/100‑1 (MQ=255)
                                                                                                aTCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt     <  1:2467221/100‑1 (MQ=255)
                                                                                                  ccTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGc   <  2:504143/100‑1 (MQ=255)
                                                                                                   cTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc  >  1:2855841/1‑100 (MQ=255)
                                                                                                   cTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc  <  1:3700054/100‑1 (MQ=255)
                                                                                                   cTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc  >  1:1779756/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
AGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCC  >  NC_000913/1931995‑1932193

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: ‑