Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,328,527 +C 1.2% intergenic (‑48/+12) entB ← / ← argF isochorismatase/ornithine carbamoyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,328,5271.→C 1.2% 3.7 ‑2.2 0.0467intergenic (‑48/+12)entB/argFisochorismatase/ornithine carbamoyltransferase
Reads supporting (aligned to +/- strand):  new base (1/3):  ref base (306/155):  total (307/160)
Fisher's exact test strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 519 total aligned reads displayed.

CGCGCCTGCATCTCGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..A..C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAATGCGGCCCG  >  NC_014034/3328427‑3328623
                                                                                                      |                                                                                                    
cgcgCCTGCATCTCGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCa.a                                                                                                       >  1:296566/1‑101 (MQ=255)
 gcgcCTGCATCTCGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCa.a..a                                                                                                    <  2:6621095/101‑1 (MQ=255)
  cgcCTGCATCTCGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCa.a..ag                                                                                                   <  2:5187188/101‑2 (MQ=255)
  cgcCTGCATCTCGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCa.a..ag                                                                                                   <  1:1540407/101‑2 (MQ=255)
      tGCATCTCGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..A..g.cccc                                                                                            >  1:8887047/1‑101 (MQ=255)
      tGCATCTCGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..A..g.cccc                                                                                            >  1:6893109/1‑101 (MQ=255)
           ctcGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGcc                                                                                        >  2:1824589/1‑101 (MQ=255)
           ctcGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGcc                                                                                        >  1:3197497/1‑101 (MQ=255)
             cGCTTTGCAGGTGGATGATCAGATGCCAGGTCGTGCCCATCGCAGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCGGCCtt                                                                                      <  2:8310960/101‑1 (MQ=21)
              gttttGCATGTCGATGACCAGAGGGGCGGTCTGCCCCACCCCCGTGCTTTTCATCTTCCCAAAAATACCGGGGGGGTCGGGGGGCA.A..AC.C.CCCCAGCCTTc                                                                                     <  2:1953172/99‑1 (MQ=255)
               cTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCa                                                                                    >  2:8646752/1‑101 (MQ=255)
               cTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCa                                                                                    <  1:5815539/101‑1 (MQ=255)
               cTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCa                                                                                    >  1:5314342/1‑101 (MQ=255)
                tttGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTccc                                                                                   >  2:7275468/1‑99 (MQ=255)
                 ttGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACt                                                                                  >  2:8149242/1‑101 (MQ=255)
                 ttGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAAAATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACt                                                                                  >  1:543274/1‑101 (MQ=255)
                  tGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGCGTCTGGTGGCA.A..AG.C.CTCCAGCCTTCACtt                                                                                 >  2:3231121/1‑101 (MQ=37)
                   gCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTg                                                                                >  1:4234370/1‑101 (MQ=255)
                   gCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTg                                                                                >  1:1409321/1‑101 (MQ=255)
                     aTGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTATCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGcc                                                                              >  2:4423456/1‑101 (MQ=255)
                     aTGTCCATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGcc                                                                              >  2:4524052/1‑101 (MQ=255)
                         cGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCaaa                                                                          >  1:7711822/1‑101 (MQ=255)
                         cGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCaaa                                                                          >  2:6025808/1‑101 (MQ=255)
                         cGATGACCGGAAGCGCGGCCTTACCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCaaa                                                                          <  1:1670264/101‑1 (MQ=25)
                          gatgatCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAg                                                                         >  2:5055938/1‑101 (MQ=255)
                            tgaacggAAGCGCGGTCTTGCCCAGCGCCGGACCTTTCATCTTGCCCTAAATCTCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGcc                                                                       <  1:4031161/95‑1 (MQ=37)
                              atCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCa                                                                     >  2:55926/1‑101 (MQ=255)
                                cAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAgc                                                                   >  1:8150920/1‑101 (MQ=255)
                                cAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAgc                                                                   >  1:2647354/1‑101 (MQ=255)
                                cAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAgc                                                                   >  1:6099022/1‑101 (MQ=255)
                                cAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAgc                                                                   >  2:2805817/1‑101 (MQ=255)
                                  gAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAgcgc                                                                 >  1:1919443/1‑101 (MQ=255)
                                   aaGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCa                                                                >  2:1153559/1‑101 (MQ=255)
                                   aaGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCa                                                                <  2:4671011/101‑1 (MQ=255)
                                    ggcgcGTTCTCCCCCTCCGCCGTCCTTTTCATCTTACCATAAAGCTGGGGGGGGTGGGGGGACACC..C..C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCAt                                                               <  1:4402984/100‑1 (MQ=255)
                                     gcgcGGTCTTGCCCATCGCCTTGCCTTTCATCTTGCCATAAATACCCTGGGGGTCTGGGGGCA.A..AC.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATg                                                              <  2:6376401/101‑1 (MQ=37)
                                         ggTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCg                                                          >  1:9189335/1‑101 (MQ=255)
                                         ggTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCg                                                          >  2:5868384/1‑101 (MQ=255)
                                         ggTCTGCCCCCTCGCCGCCCTTTTCTTCTCCCCA.AAATTTCCGGGGGGGGGGGGGGCA.ACCC..C.CCCCGGCCTTCACTTGCCCAAAGCCCAGCGCATGATCg                                                          <  2:3023779/101‑1 (MQ=255)
                                          gTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGc                                                         >  1:753367/1‑101 (MQ=255)
                                            cTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCt                                                       >  1:5592235/1‑101 (MQ=255)
                                            cTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCt                                                       >  1:12286/1‑101 (MQ=255)
                                             ttccccATCGCCGTCCCTTTCATCTTGCCATAAATATCCGGGGGGTCTGGGGGCA.A..AC.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCtt                                                      <  1:6242340/98‑1 (MQ=25)
                                             ttGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCtt                                                      >  1:3625762/1‑101 (MQ=255)
                                             ttGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCtt                                                      >  2:2166202/1‑101 (MQ=255)
                                              tGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCttt                                                     >  2:6456998/1‑101 (MQ=255)
                                              tGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTACCCAAAGCCCAGCGCATGATCGCCttt                                                     >  2:8011931/1‑101 (MQ=255)
                                               gCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCtttt                                                    >  1:2821231/1‑101 (MQ=255)
                                                   aTCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTgcgc                                                <  2:7444299/101‑1 (MQ=255)
                                                    ccccccttccTTTCATCCTTCCAATAATTTCCGGGGGGGTGGGGGGCA.A..ACCC.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTgcgc                                                <  2:4332210/93‑1 (MQ=255)
                                                     cGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTgcgaga                                              >  2:2302803/1‑98 (MQ=255)
                                                     cGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCAt                                              <  2:4948388/101‑1 (MQ=255)
                                                     cGCCGTCCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCAt                                              <  2:2454001/101‑1 (MQ=255)
                                                     cGCCGGGCTTTTCATCTTGCAATAAATATCCGGGGGGTCGGGGGGCA.A..AC.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCAt                                              <  1:221660/101‑1 (MQ=255)
                                                      gCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATg                                             >  2:3341736/1‑101 (MQ=255)
                                                       ccttGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGc                                            <  1:8575552/98‑1 (MQ=255)
                                                       ccGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGc                                            >  2:6193056/1‑101 (MQ=255)
                                                        cGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCa                                           >  1:1051888/1‑101 (MQ=255)
                                                        cGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCGTTTGCGCATGCa                                           <  1:6531475/101‑1 (MQ=255)
                                                         gTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAg                                          >  1:3770938/1‑101 (MQ=255)
                                                         gTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAg                                          >  2:3575416/1‑101 (MQ=255)
                                                          tGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATAGCATTTTGCGCATGCAGc                                         >  2:4068127/1‑101 (MQ=37)
                                                           tccTTTCATCTTGCCATAAATATCCGGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGcc                                        <  2:7059729/100‑1 (MQ=255)
                                                           gCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGcc                                        >  2:906371/1‑101 (MQ=255)
                                                           gCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGcc                                        >  1:3321517/1‑101 (MQ=255)
                                                             cTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTCGCGCATGAAGCCgg                                      >  2:1413223/1‑101 (MQ=37)
                                                               ttCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGtt                                    <  1:7940081/101‑1 (MQ=255)
                                                                tCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTtc                                   <  2:2659180/101‑1 (MQ=255)
                                                                 cATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTtct                                  >  2:6278907/1‑101 (MQ=255)
                                                                 cATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTtct                                  <  2:6904497/101‑1 (MQ=255)
                                                                  aTCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTtctc                                 >  2:5390059/1‑101 (MQ=255)
                                                                  aTCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTtctc                                 >  1:4987358/1‑101 (MQ=255)
                                                                    cTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCg                               >  1:5142766/1‑101 (MQ=255)
                                                                     ttGCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGc                              <  1:4587860/101‑1 (MQ=255)
                                                                       gCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCt                            >  1:575312/1‑101 (MQ=255)
                                                                       gCCATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCt                            >  1:4834921/1‑101 (MQ=255)
                                                                        ccATAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCtc                           <  1:3911152/101‑1 (MQ=255)
                                                                           tAAATATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCtcgtc                        >  1:8771577/1‑101 (MQ=255)
                                                                            aaaTATCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCtcgtcg                       >  2:1584546/1‑101 (MQ=255)
                                                                              atatCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGaa                     <  2:8838020/101‑1 (MQ=255)
                                                                              atatCCTGGGGGGCTGGGGGCA.A..AG.C.CCGCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGGCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGaa                     >  1:8082813/1‑101 (MQ=255)
                                                                               tatCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAg                    >  1:6757245/1‑101 (MQ=255)
                                                                               tatCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGGCGatg                    >  2:4200777/1‑99 (MQ=255)
                                                                                ttCCTGGGGGACTGAGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGa                   <  2:1661889/100‑1 (MQ=255)
                                                                                atCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGa                   >  2:9175535/1‑101 (MQ=255)
                                                                                atCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGa                   >  2:6541293/1‑101 (MQ=255)
                                                                                atCCTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCCGTTCTCCGCCTCGTCGAAGa                   >  1:8815303/1‑101 (MQ=255)
                                                                                 tCCGGGGGGTCGGGGGGCA.A..AC.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGAt                  <  1:8417351/101‑1 (MQ=255)
                                                                                   cTGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCa                >  1:2507284/1‑101 (MQ=255)
                                                                                    tGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCAc               <  2:6823807/101‑1 (MQ=255)
                                                                                    tGGGGGTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTATCCGCCTCGTCGAAGATCAc               >  1:7212453/1‑101 (MQ=255)
                                                                                      ggggTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCg             >  2:7807892/1‑101 (MQ=255)
                                                                                         gTCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAAt          >  2:534878/1‑101 (MQ=255)
                                                                                         gTCTGGGGGCA.A..A..CACCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAAt          <  2:7582093/101‑1 (MQ=255)
                                                                                          tCTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAATg         >  2:5507815/1‑101 (MQ=255)
                                                                                           cTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAATGc        >  1:2398935/1‑101 (MQ=255)
                                                                                           cTGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAATGc        >  1:2101511/1‑101 (MQ=255)
                                                                                            tGGGGGCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAATGCg       >  2:8193974/1‑101 (MQ=255)
                                                                                             gggggCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAATGCgg      >  2:896930/1‑101 (MQ=255)
                                                                                               gggCA.A..AG.C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAATGCGGcc    <  2:8460690/101‑1 (MQ=255)
                                                                                                      |                                                                                                    
CGCGCCTGCATCTCGCTTTGCATGTCGATGATCAGAAGCGCGGTCTTGCCCATCGCCGTGCCTTTCATCTTGCCATAAATATCCTGGGGGTCTGGGGGCA.A..A..C.CCCCAGCCTTCACTTGCCCAAAGCCCAGCGCATGATCGCCTTTTGCGCATGCAGCCGGTTCTCCGCCTCGTCGAAGATCACCGAATGCGGCCCG  >  NC_014034/3328427‑3328623

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