Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 279,645 +A 3.0% intergenic (‑47/‑33) RCAP_rcc00242 ← / → dapD hypothetical protein/2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑succinyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034279,6450.→A 3.0% 31.0 0.0 ‑0.3336intergenic (‑47/‑33)RCAP_rcc00242/dapDhypothetical protein/2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑succinyltransferase
Reads supporting (aligned to +/- strand):  new base (5/5):  ref base (172/154):  total (177/159)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.41e-01

Only 100 of 364 total aligned reads displayed.

TTTTGCATCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAACCCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCTCACGCC  >  NC_014034/279545‑279742
                                                                                                                                                                                                      
 tttGCATCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAAcc                                                                                                  >  2:3832635/1‑101 (MQ=255)
  ttGCATCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATCTTCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAAcca                                                                                                 >  1:3101837/1‑100 (MQ=255)
   tGCATCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAAccat                                                                                                >  2:699057/1‑99 (MQ=255)
    gCATCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAAccatg                                                                                               >  2:2598299/1‑98 (MQ=255)
    gCATCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAAccatg                                                                                               >  1:2583574/1‑98 (MQ=255)
     cATCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCGTATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAAccatga                                                                                              <  1:4880614/101‑5 (MQ=255)
      aTCCATGATGCTGTCACGGAACGGGTGGGAACGTCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGaa                                                                                            <  1:4342298/101‑1 (MQ=255)
       tCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGGTTTCCCCAACTATAGGGGGAAAA.CCATGAAt                                                                                           <  1:4554084/101‑1 (MQ=255)
       tCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAAt                                                                                           <  2:4047513/101‑1 (MQ=255)
        ccATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATg                                                                                          <  1:2424081/101‑1 (MQ=255)
         cATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGc                                                                                         >  2:4315532/1‑101 (MQ=255)
         cATGATACTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGc                                                                                         >  1:4870813/1‑101 (MQ=255)
          atgatgCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGcc                                                                                        <  2:3320846/101‑1 (MQ=255)
           tgatgCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCa                                                                                       >  2:1825683/1‑101 (MQ=255)
           tgatgCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCa                                                                                       >  2:921693/1‑101 (MQ=255)
            gatgCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCAc                                                                                      <  1:177497/101‑1 (MQ=255)
             atgCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACt                                                                                     >  2:4943060/1‑101 (MQ=255)
             atgCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACt                                                                                     <  1:4227923/101‑1 (MQ=255)
              tgCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACtt                                                                                    <  2:3036769/101‑1 (MQ=255)
              tgCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACtt                                                                                    <  1:2725583/101‑1 (MQ=255)
               gCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTg                                                                                   <  1:360604/101‑1 (MQ=255)
                cTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGa                                                                                  >  1:1300243/1‑101 (MQ=255)
                cTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGa                                                                                  >  1:1632109/1‑101 (MQ=255)
                 tGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGaa                                                                                 <  1:1673637/101‑1 (MQ=255)
                  gTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGaaa                                                                                <  1:601442/101‑1 (MQ=255)
                  gTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGaaa                                                                                >  1:4102438/1‑101 (MQ=255)
                  gTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGaaa                                                                                >  2:3542123/1‑101 (MQ=255)
                    cACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAaga                                                                              >  2:1458210/1‑101 (MQ=255)
                       ggAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGa                                                                           <  1:2018448/101‑1 (MQ=255)
                       ggAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGa                                                                           <  1:2479150/101‑1 (MQ=255)
                         aaCGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGAcc                                                                         >  2:1592565/1‑101 (MQ=255)
                          aCGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCg                                                                        <  1:1129197/101‑1 (MQ=255)
                                 ggAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGCAAGAGGAGCGTTCCGag                                                                 >  1:1453580/1‑100 (MQ=255)
                                 ggAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGAt                                                                 >  1:3454241/1‑101 (MQ=255)
                                   aaCGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGt                                                               <  1:1401741/101‑1 (MQ=255)
                                    aCGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTc                                                              <  2:2372575/101‑1 (MQ=255)
                                    aCGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTc                                                              <  1:3556924/101‑1 (MQ=255)
                                      ggCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGa                                                            <  2:2151827/101‑1 (MQ=255)
                                      ggCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGa                                                            <  2:2003442/101‑1 (MQ=255)
                                        cTTTCGTCTTTCATATCNTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAAc                                                          <  2:77934/101‑1 (MQ=255)
                                        cTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAAc                                                          <  2:2487454/101‑1 (MQ=255)
                                        cTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAAc                                                          <  2:2649203/101‑1 (MQ=255)
                                         tttCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACg                                                         >  1:2188388/1‑101 (MQ=255)
                                         tttCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACg                                                         >  1:4704455/1‑101 (MQ=255)
                                          ttCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGGACGc                                                        >  2:3568333/1‑101 (MQ=255)
                                          ttCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGc                                                        >  1:4449326/1‑101 (MQ=255)
                                          ttCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGc                                                        <  1:3225030/101‑1 (MQ=255)
                                           tCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGcc                                                       >  2:2903560/1‑101 (MQ=255)
                                           tCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGcc                                                       >  1:4636571/1‑101 (MQ=255)
                                           tCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGcc                                                       <  1:3832920/101‑1 (MQ=255)
                                            cGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCa                                                      <  2:1209918/101‑1 (MQ=255)
                                             gTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCAt                                                     <  2:596723/101‑1 (MQ=255)
                                              tCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATg                                                    <  1:4014697/101‑1 (MQ=255)
                                               cTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGTACCGTTCCGATGTCGAACGCCATGc                                                   >  2:5005760/1‑101 (MQ=255)
                                               cTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGc                                                   >  2:34417/1‑101 (MQ=255)
                                               cTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGc                                                   >  1:2622754/1‑101 (MQ=255)
                                                tttCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCt                                                  >  2:4173260/1‑101 (MQ=255)
                                                tttCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCt                                                  >  1:2289556/1‑101 (MQ=255)
                                                  tCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCg                                                >  2:3503659/1‑101 (MQ=255)
                                                  tCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCg                                                >  2:3812181/1‑101 (MQ=255)
                                                   cATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGa                                               >  1:1976006/1‑101 (MQ=255)
                                                    atatCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGaa                                              <  2:1421547/101‑1 (MQ=255)
                                                    atatCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATCCTCGaa                                              >  2:4693003/1‑101 (MQ=255)
                                                     tatCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAg                                             <  1:4777481/101‑1 (MQ=255)
                                                      atCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGc                                            <  2:3654146/101‑1 (MQ=255)
                                                      atCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGc                                            >  1:2904614/1‑101 (MQ=255)
                                                        ccTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGATGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTg                                          <  2:3405355/101‑1 (MQ=255)
                                                        ccTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTg                                          >  1:3421336/1‑101 (MQ=255)
                                                        ccTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTg                                          >  1:3908614/1‑101 (MQ=255)
                                                          tGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGcc                                        <  1:2801991/101‑1 (MQ=255)
                                                            cTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCAt                                      <  2:3137647/101‑1 (MQ=255)
                                                             tCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGTTCCCATc                                     >  1:3044844/1‑101 (MQ=255)
                                                             tCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATc                                     <  1:312854/101‑1 (MQ=255)
                                                             tCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATc                                     >  1:4140706/1‑101 (MQ=255)
                                                             tCATTGCGCGAGGAGTTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATc                                     <  1:719961/101‑1 (MQ=255)
                                                               aTTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGa                                   >  1:2237823/1‑101 (MQ=255)
                                                                 tGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAg                                 <  2:1664250/101‑1 (MQ=255)
                                                                  gcgcGAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAgc                                <  1:505848/101‑1 (MQ=255)
                                                                   cgcgAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAgcc                               <  1:3503929/101‑1 (MQ=255)
                                                                      gAGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAgccgcc                            >  2:5047949/1‑101 (MQ=255)
                                                                       aGGGGCTTTCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCt                           <  2:4722640/101‑1 (MQ=255)
                                                                              ttCCCCAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGt                    >  2:3904520/1‑101 (MQ=255)
                                                                                 cccAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTcg                 <  1:3887795/101‑1 (MQ=255)
                                                                                 cccAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTcg                 <  1:3333578/101‑1 (MQ=255)
                                                                                 cccAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTcg                 >  2:1892062/1‑101 (MQ=255)
                                                                                  ccAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTcgc                >  1:2859823/1‑101 (MQ=255)
                                                                                   cAACTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTcgcg               <  2:1075998/101‑1 (MQ=255)
                                                                                    aaCTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGa              >  2:3259669/1‑101 (MQ=255)
                                                                                    aaCTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGa              >  2:2727916/1‑101 (MQ=255)
                                                                                    aaCTATAGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGa              <  1:1012855/101‑1 (MQ=255)
                                                                                        ataGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACAcc          >  2:2763578/1‑101 (MQ=255)
                                                                                         taGGGGGAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACAccc         >  1:1879350/1‑101 (MQ=255)
                                                                                           gggggAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCTc       >  1:944086/1‑101 (MQ=255)
                                                                                            ggggAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCTCa      <  2:1300122/101‑1 (MQ=255)
                                                                                             gggAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCTCAc     <  2:4566895/101‑1 (MQ=255)
                                                                                             gggAAAA.CCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCTCAc     <  2:2237668/101‑1 (MQ=255)
                                                                                              gggaaaGCCATTAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCTCAc     <  2:2467116/98‑1 (MQ=255)
                                                                                               ggaaaACCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCTCACg    >  2:3965442/3‑101 (MQ=255)
                                                                                               ggaaaACCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCTCACg    <  2:2821618/99‑1 (MQ=255)
                                                                                                gaaaACCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCCCACGc   >  1:80126/2‑101 (MQ=255)
                                                                                                                                                                                                      
TTTTGCATCCATGATGCTGTCACGGAACGGGTGGGAACGGCTTTCGTCTTTCATATCCTGCTCATTGCGCGAGGGGCTTTCCCCAACTATAGGGGGAAAACCCATGAATGCCACTTGAAAGAGGACCGTTCCGATGTCGAACGCCATGCTCGAAGCTGCCATCGAAGCCGCCTGGGAGGTTCGCGACACCCTCACGCC  >  NC_014034/279545‑279742

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG