Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,095 T→A 16.7% 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,0950TA16.7% 321.7 / 53.7 132intergenic (+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 (62/48);  new base A (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 = 6.46e-01

Only 100 of 132 total aligned reads displayed.

GTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCC  >  NC_000913/1931998‑1932193
                                                                                                 |                                                                                                   
gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCg                                                                                                  >  1:711719/1‑100 (MQ=255)
  gttgttGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGAt                                                                                                <  1:610862/100‑1 (MQ=255)
    tgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCg                                                                                              >  1:504148/1‑100 (MQ=255)
    tgtggACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTAAAAAAAAAAAAAAAAACGCTACAAAAATGCCCGATCCTCGATc                                                                                               <  2:931282/99‑1 (MQ=255)
       tgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAAACAAGAAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGc                                                                                           <  2:2998312/100‑1 (MQ=255)
        gccGCAATTGAACGCGCGAAGCACGCCGCCGGACAGTTAAAAGTACAGGGTCAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGc                                                                                           <  1:63595/97‑1 (MQ=255)
        gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                          <  1:1356667/100‑1 (MQ=255)
        gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                                                          <  1:2440163/100‑1 (MQ=255)
         aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAt                                                                                         <  1:2106981/100‑1 (MQ=39)
         aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGACCGGGCAt                                                                                         >  1:1025846/1‑100 (MQ=255)
          cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtt                                                                                        >  1:2763760/1‑100 (MQ=255)
          cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtt                                                                                        >  1:3729754/1‑100 (MQ=255)
          cGCCATTGAACGAGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCAGATCggaaga                                                                                         >  2:1268324/1‑96 (MQ=255)
           gCCATTGAACGCGCGAAGCACGCGGCCGGACAGGAAAAAGTACAGGGTTAAACCCGGAAAAAAAAAAAAAAAAAAATGCCCGATCGAGGATCGGGCAtt                                                                                        <  2:1821630/99‑1 (MQ=21)
           gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAt                                                                                         <  2:1312000/98‑1 (MQ=255)
            ccATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATAGGGCAtt                                                                                        >  2:238008/1‑98 (MQ=255)
               ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGAc                                                                                   >  2:2136526/1‑100 (MQ=255)
                tGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACt                                                                                  >  1:1716299/1‑100 (MQ=255)
                 gAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtttttgtag                                                                                 <  1:2818930/100‑6 (MQ=37)
                  aaCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACttt                                                                                >  1:2808839/1‑100 (MQ=255)
                   aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtttt                                                                               <  2:412938/100‑1 (MQ=255)
                   aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtttt                                                                               >  1:1791699/1‑100 (MQ=255)
                    cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                              >  1:1092469/1‑100 (MQ=255)
                    cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                              >  2:2671955/1‑100 (MQ=255)
                    cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGACCGATCCTCGATCGGGCATTTTGACTTTTa                                                                              >  2:3090917/1‑100 (MQ=255)
                    cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCAACAAAAATGAACGATCCTCAAACGGGCAAATTGACTTTTa                                                                              >  1:1268553/1‑100 (MQ=255)
                      cgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACa                                                                            >  2:1686472/1‑100 (MQ=255)
                            gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAg                                                                      >  1:969272/1‑100 (MQ=255)
                            gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAg                                                                      >  1:1822223/1‑100 (MQ=255)
                             cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                     >  2:176268/1‑100 (MQ=255)
                              aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcg                                                                    <  1:2609344/100‑1 (MQ=255)
                                   gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTc                                                               >  1:2138053/1‑100 (MQ=255)
                                   gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTc                                                               <  1:770509/100‑1 (MQ=255)
                                    ccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCt                                                              >  2:877514/1‑100 (MQ=255)
                                    ccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCt                                                              <  1:3508681/100‑1 (MQ=255)
                                      ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                            <  2:1437915/100‑1 (MQ=255)
                                      ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                            <  1:2186225/100‑1 (MQ=255)
                                           ggTAAAAGTACAGGGTTAAACCCGGCAAAACAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACCGCt                                                        >  2:1194535/1‑99 (MQ=255)
                                           ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                       >  1:1553190/1‑100 (MQ=255)
                                           ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                       >  1:3355377/1‑100 (MQ=255)
                                            gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                      >  1:3444797/1‑100 (MQ=255)
                                               aaaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcg                                                   >  1:1383968/1‑100 (MQ=255)
                                                aaGTACAGGGTTAAACCCGGCAAAAAAACGCCACAAAAATGCCCGATCCTCGATCGGGCATTTCCACTTTTACAGCTTAGCGCCTTCTACAGCTTCCcgc                                                  >  1:2450091/1‑100 (MQ=255)
                                                 aGTACAGGGTTAAACCCGGCAAAAACACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgcg                                                 <  1:2473669/100‑1 (MQ=255)
                                                 aGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgcg                                                 <  1:1950937/100‑1 (MQ=255)
                                                  gTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTAACTTTTACAGCTTACCGCCTTCTACAGCTTCCCgcgc                                                >  2:3237767/1‑100 (MQ=255)
                                                    aCAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCa                                              >  2:619015/1‑100 (MQ=255)
                                                     cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                             >  1:1544578/1‑100 (MQ=255)
                                                       gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGc                                            <  1:1722579/99‑1 (MQ=255)
                                                       gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                           <  2:1943581/100‑1 (MQ=255)
                                                       gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                           >  1:1056819/1‑100 (MQ=255)
                                                       gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                           <  2:3159493/100‑1 (MQ=255)
                                                        ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                          <  2:210948/100‑1 (MQ=255)
                                                        ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                          >  2:14256/1‑100 (MQ=255)
                                                        ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                          <  2:1006724/100‑1 (MQ=255)
                                                         gTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                         <  1:2862951/100‑1 (MQ=255)
                                                          ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAg                                        <  1:3200360/100‑1 (MQ=255)
                                                             aaCCCGGCAAAAAAACGCTACAAAAATGCCTGAGCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCACGCTTAGTaa                                     >  2:2503084/1‑100 (MQ=255)
                                                              aCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAAt                                    >  1:3028809/1‑100 (MQ=255)
                                                               nccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATg                                   >  2:3117218/2‑100 (MQ=255)
                                                               cccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATg                                   >  1:1910454/1‑100 (MQ=255)
                                                                ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGc                                  >  1:912196/1‑100 (MQ=255)
                                                                   gCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGt                               >  1:1028282/1‑100 (MQ=255)
                                                                     aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCg                             <  1:2795402/100‑1 (MQ=255)
                                                                      aaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGt                            <  1:360912/100‑1 (MQ=255)
                                                                       aaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTa                           <  1:256800/100‑1 (MQ=255)
                                                                        aaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTaa                          <  1:2866643/100‑1 (MQ=255)
                                                                          aaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                        <  2:3325140/100‑1 (MQ=255)
                                                                          aaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                        >  1:2069270/1‑100 (MQ=255)
                                                                            cGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACTGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGc                      >  2:2582080/1‑100 (MQ=255)
                                                                            cGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGc                      >  2:1258742/1‑100 (MQ=255)
                                                                              aagtcaaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGcc                     >  2:2730022/6‑99 (MQ=37)
                                                                               tACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                   <  1:830669/100‑1 (MQ=255)
                                                                               tACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                   >  2:3766280/1‑100 (MQ=255)
                                                                                gtcaaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                  >  2:3321583/4‑100 (MQ=38)
                                                                                aCAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                  >  1:1552301/1‑100 (MQ=255)
                                                                                aCAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                  >  1:290663/1‑100 (MQ=255)
                                                                                 tcaaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCt                 >  1:2075618/3‑100 (MQ=38)
                                                                                 tcaaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCGTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCt                 >  1:171457/3‑100 (MQ=38)
                                                                                 cAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCCCCCGc                  >  2:3551194/1‑99 (MQ=255)
                                                                                  caaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCtt                >  1:2130442/2‑100 (MQ=38)
                                                                                  caaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCGTAGTAATGCGGTCGTAATCGccc                    >  2:2409670/2‑96 (MQ=38)
                                                                                   aaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc               >  2:348453/1‑100 (MQ=38)
                                                                                   aaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc               <  2:234358/100‑1 (MQ=255)
                                                                                   aaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc               >  2:2667023/1‑100 (MQ=255)
                                                                                      aTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa             >  2:3500736/1‑99 (MQ=39)
                                                                                      aTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg            >  1:833512/1‑100 (MQ=39)
                                                                                      aTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg            >  2:3115871/1‑100 (MQ=39)
                                                                                      aTGCCCGATCGAGGAGCGGGCATTTTGACTTTTACAGCTTAGCTCGTTCCACAGCTTCAAGCGCCAGCGTAGTAATGCGGCCGTAATCGCCCGCGTCCa             >  2:3068943/1‑99 (MQ=255)
                                                                                       tGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc           <  1:572921/100‑1 (MQ=255)
                                                                                        gCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgcg          >  1:3463003/1‑100 (MQ=255)
                                                                                          ccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgcgc         >  2:1845536/1‑99 (MQ=255)
                                                                                          ccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa        <  1:3813731/100‑1 (MQ=255)
                                                                                          ccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa        <  1:1965115/100‑1 (MQ=255)
                                                                                           cGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCAt       <  1:36945/100‑1 (MQ=255)
                                                                                             aTCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt     <  1:2467221/100‑1 (MQ=255)
                                                                                              tCCTCGATCGAGCATCTTGACTTTTACAGCTTAGCGCGTTCTACAGCTTCACGCGCCAGCTGAGTAATGCGGTCGTAATCGTCCGCGTCCAGCGCATCTg    <  2:1722937/100‑1 (MQ=255)
                                                                                               ccTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGc   <  2:504143/100‑1 (MQ=255)
                                                                                               ccTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGc   >  1:3312147/1‑100 (MQ=255)
                                                                                                cTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc  >  1:1779756/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCC  >  NC_000913/1931998‑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: