Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,096 C→G 16.7% intergenic (+37/+19) 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,0960CG16.7% 351.0 / 59.4 132intergenic (+37/+19)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 (62/48);  new base G (13/9);  total (75/57)
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.62e-01

Only 100 of 135 total aligned reads displayed.

GTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGG  >  NC_000913/1931998‑1932195
                                                                                                  |                                                                                                    
gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑g                                                                                                    <  1:2482844/100‑1 (MQ=255)
gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑g                                                                                                    >  1:711719/1‑100 (MQ=255)
  gttgttGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GAt                                                                                                  <  1:610862/100‑1 (MQ=255)
    tgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCg                                                                                                >  1:504148/1‑100 (MQ=255)
    tgtggACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTAAAAAAAAAAAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATc                                                                                                 <  2:931282/99‑1 (MQ=255)
       tgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAAACAAGAAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGc                                                                                             <  2:2998312/100‑1 (MQ=255)
        gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAG‑GATCGGGCa                                                                                            <  1:1847630/100‑1 (MQ=255)
        gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAG‑GATCGGGCa                                                                                            <  1:1356667/100‑1 (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)
          cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCAtt                                                                                          >  1:2763760/1‑100 (MQ=255)
          cGCCATTGAACGAGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCAGATCggaaga                                                                                           >  2:1268324/1‑96 (MQ=255)
           gcgaTTGAACGTGGGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAGAAAACGCTACAAAAATGCCCGATCGAG‑GATCGGGCAtt                                                                                          <  2:2351835/96‑1 (MQ=37)
             cATTGAACGCGGGAGGCACGCCGCCAGAAAGGGAAAAGTACAGGGGTAAACCCGGCAAAAAAACGCTCCAAAAATGCCCGATCGAG‑GATCGGGCAtttt                                                                                        <  2:1507620/99‑1 (MQ=25)
              atagAACGCGCCAACCACGCGCCCGGCCAGGTAAAAGTACAGCGTTAATCCCGGCAAAAAAACGCTACAAACATGCCCGATCCTC‑GATCGGGCATTTTg                                                                                       <  2:339608/96‑1 (MQ=255)
              aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGa                                                                                      <  1:2422151/100‑1 (MQ=255)
               ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGAc                                                                                     >  2:2136526/1‑100 (MQ=255)
               ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGAc                                                                                     <  2:3099232/100‑1 (MQ=255)
                tGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACt                                                                                    >  1:1716299/1‑100 (MQ=255)
                 gAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAG‑GATCGGGCAtttttgtag                                                                                   <  1:2818930/100‑6 (MQ=37)
                 gAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGa                                                                                      <  2:3217232/97‑1 (MQ=255)
                 gAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACtt                                                                                   <  1:2109908/100‑1 (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                                                                                >  2:2671955/1‑100 (MQ=255)
                    cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCAACAAAAATGAACGATCCTC‑AAACGGGCAAATTGACTTTTa                                                                                >  1:1268553/1‑100 (MQ=255)
                     gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTAc                                                                               <  2:3561028/100‑1 (MQ=255)
                        cgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGc                                                                            >  1:902969/1‑100 (MQ=255)
                             cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAgc                                                                       <  1:1200842/100‑1 (MQ=255)
                             cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAgc                                                                       >  2:176268/1‑100 (MQ=255)
                              aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAgcg                                                                      <  1:2609344/100‑1 (MQ=255)
                              aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAgcg                                                                      <  1:2063067/100‑1 (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)
                                           ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                         >  1:3355377/1‑100 (MQ=255)
                                            gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                        >  1:3444797/1‑100 (MQ=255)
                                            gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                        >  1:1700967/1‑100 (MQ=255)
                                                aaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgc                                                    >  2:1276496/1‑100 (MQ=255)
                                                 aGTACAGGGTTAAACCCGGCAAAAACACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgcg                                                   <  1:2473669/100‑1 (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)
                                                     cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                               >  1:1544578/1‑100 (MQ=255)
                                                     cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                               >  2:572915/1‑100 (MQ=255)
                                                       gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                             <  2:1943581/100‑1 (MQ=255)
                                                       gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                             >  1:1056819/1‑100 (MQ=255)
                                                        ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                            <  2:1006724/100‑1 (MQ=255)
                                                        ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                            >  2:14256/1‑100 (MQ=255)
                                                          ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAg                                          <  1:3200360/100‑1 (MQ=255)
                                                           tAAAGCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGt                                         <  1:143501/100‑1 (MQ=255)
                                                           tAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGt                                         <  1:2932200/100‑1 (MQ=255)
                                                             aaCCCGGCAAAAAAACGCTACAAAAATGCCTGAGCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCACGCTTAGTaa                                       >  2:2503084/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)
                                                                      aaaaaaCGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGt                              <  1:360912/100‑1 (MQ=255)
                                                                        aaaaCGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTaa                            <  1:2866643/100‑1 (MQ=255)
                                                                          aaCGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                          >  1:2069270/1‑100 (MQ=255)
                                                                            cGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACTGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGc                        >  2:2582080/1‑100 (MQ=255)
                                                                            cGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGc                        >  2:1258742/1‑100 (MQ=255)
                                                                               tACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                     >  2:3766280/1‑100 (MQ=255)
                                                                                gtcaaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                    >  2:3321583/4‑100 (MQ=38)
                                                                                aCAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                    >  1:1552301/1‑100 (MQ=255)
                                                                                aCAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                    >  1:290663/1‑100 (MQ=255)
                                                                                 tcaaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCGTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCt                   >  1:171457/3‑100 (MQ=38)
                                                                                  caaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCtt                  >  1:2130442/2‑100 (MQ=38)
                                                                                  caaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCGTAGTAATGCGGTCGTAATCGccc                      >  2:2409670/2‑96 (MQ=38)
                                                                                   aaaaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                 >  2:348453/1‑100 (MQ=38)
                                                                                   aaaaTGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                 >  2:2667023/1‑100 (MQ=255)
                                                                                   aaaaTGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                 <  2:234358/100‑1 (MQ=255)
                                                                                   aaaaTGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                 >  1:2791739/1‑100 (MQ=255)
                                                                                     aaTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTcc                >  1:748623/1‑99 (MQ=39)
                                                                                      aTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa               >  2:3500736/1‑99 (MQ=39)
                                                                                      aTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg              >  1:833512/1‑100 (MQ=39)
                                                                                      aTGCCCGATCGAG‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg              >  2:3115871/1‑100 (MQ=39)
                                                                                      aTGCCCGATCGAG‑GAGCGGGCATTTTGACTTTTACAGCTTAGCTCGTTCCACAGCTTCAAGCGCCAGCGTAGTAATGCGGCCGTAATCGCCCGCGTCCa               >  2:3068943/1‑99 (MQ=255)
                                                                                       tGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc             <  1:572921/100‑1 (MQ=255)
                                                                                        gCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgcg            >  1:3463003/1‑100 (MQ=255)
                                                                                          ccGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgcgc           >  2:1845536/1‑99 (MQ=255)
                                                                                          ccGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa          <  1:3813731/100‑1 (MQ=255)
                                                                                          ccGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa          >  1:176836/1‑100 (MQ=255)
                                                                                           cGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCAt         <  1:36945/100‑1 (MQ=255)
                                                                                             aTCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt       <  1:2467221/100‑1 (MQ=255)
                                                                                              tCCTC‑GATCGAGCATCTTGACTTTTACAGCTTAGCGCGTTCTACAGCTTCACGCGCCAGCTGAGTAATGCGGTCGTAATCGTCCGCGTCCAGCGCATCTg      <  2:1722937/100‑1 (MQ=255)
                                                                                               ccTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGc     >  1:3312147/1‑100 (MQ=255)
                                                                                               ccTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGc     <  2:504143/100‑1 (MQ=255)
                                                                                                cTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc    >  1:2855841/1‑100 (MQ=255)
                                                                                                cTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc    >  1:1779756/1‑100 (MQ=255)
                                                                                                  c‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCgg  <  2:3580364/100‑1 (MQ=255)
                                                                                                  c‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCgg  <  1:586221/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
GTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTC‑GATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGG  >  NC_000913/1931998‑1932195

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