Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,301,060 +C 93.6% intergenic (+17/‑92) urtA → / → urtC urea ABC transporter permease UrtA/urea ABC transporter permease UrtC

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,301,0601.→C93.6% 21.8 ‑3.2 ‑0.6140intergenic (+17/‑92)urtA/urtCurea ABC transporter permease UrtA/urea ABC transporter permease UrtC
Reads supporting (aligned to +/- strand):  new base (77/54):  ref base (2/7):  total (79/61)
Fisher's exact test strand distribution p-value = 4.11e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.42e-01

Only 100 of 158 total aligned reads displayed.

TGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.GCG.CC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGTCCCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACC  >  NC_014034/1300960‑1301157
                                                                                                      |                                                                                                   
tGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.ccg                                                                                                      <  1:3256790/101‑1 (MQ=255)
tGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.ccg                                                                                                      >  1:352852/1‑101 (MQ=255)
 ggATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCAGATGCCCCCCG.ccg.c                                                                                                    >  1:3936996/1‑101 (MQ=255)
  gatgatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCG.cc                                                                                                   <  2:625746/101‑1 (MQ=255)
  gatgatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCG.cc                                                                                                   <  2:3408840/101‑1 (MQ=255)
   atgatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCG.ccc                                                                                                  >  2:2731523/1‑100 (MQ=255)
   atgatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCG.ccc                                                                                                  <  2:3038671/101‑2 (MQ=255)
    tgatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCG.ccct                                                                                                 >  1:2336911/1‑99 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCG.CC.ctt                                                                                               <  1:2915595/101‑1 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCG.CC.ctt                                                                                               <  1:2153269/101‑1 (MQ=255)
      atCCTCTTCATCATCATCTTCATCCAGTTCCGCAGGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCG.CC.cttt                                                                                              <  1:4072449/101‑2 (MQ=255)
       tCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCG.CC.ctttc                                                                                             <  1:4026433/101‑3 (MQ=255)
       tCCGCTTCATCATCATCTTCATCTAGTTCCGCCCTCGCGACATCATCGAGCAGCGCGTGCACGCCGCGGGAGATTGTGCCGATCCCCCCCG.CCG.CC.ctttc                                                                                             <  2:653416/101‑3 (MQ=255)
        ccTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.a                                                                                            <  1:1671837/101‑1 (MQ=255)
        ccTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.a                                                                                            <  2:1354306/101‑1 (MQ=255)
        ccTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.a                                                                                            >  1:3875356/1‑101 (MQ=255)
        ccTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.a                                                                                            >  2:1128843/1‑101 (MQ=255)
        ccTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.a                                                                                            <  1:2958431/101‑1 (MQ=255)
            ttCATCATCATCTTCATNCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCtt                                                                                        <  2:1497426/101‑1 (MQ=255)
            ttCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCtt                                                                                        >  1:3690225/1‑101 (MQ=255)
             tcatcatcaTCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTg                                                                                       <  1:1690899/101‑1 (MQ=255)
               atcatcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGcc                                                                                     >  1:3776642/1‑101 (MQ=255)
                tcatcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGccc                                                                                    <  2:4483787/101‑1 (MQ=255)
                tcatcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGccc                                                                                    >  2:4131907/1‑101 (MQ=255)
                 catcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCa                                                                                   <  1:3780397/101‑1 (MQ=255)
                 catcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCa                                                                                   >  2:1222446/1‑101 (MQ=255)
                 catcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCa                                                                                   >  2:2604283/1‑101 (MQ=255)
                  atcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaa                                                                                  <  2:2401184/101‑1 (MQ=255)
                  atcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaa                                                                                  <  2:2207318/101‑1 (MQ=255)
                  atcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaa                                                                                  >  2:1249306/1‑101 (MQ=255)
                  atcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaa                                                                                  <  1:693949/101‑1 (MQ=255)
                  atcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaa                                                                                  <  2:1848156/101‑1 (MQ=255)
                    catcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaaaa                                                                                <  1:3200543/101‑1 (MQ=255)
                    catcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaaaa                                                                                >  1:1575326/1‑101 (MQ=255)
                    catcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaaaa                                                                                <  1:1704748/101‑1 (MQ=255)
                    catcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaaaa                                                                                >  1:4567592/1‑101 (MQ=255)
                    catcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCaaaa                                                                                <  1:4430259/101‑1 (MQ=255)
                     atcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAat                                                                               >  2:547907/1‑101 (MQ=255)
                     atcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAat                                                                               >  2:1946805/1‑101 (MQ=255)
                      tcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAata                                                                              >  1:1288192/1‑101 (MQ=255)
                       cTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAatat                                                                             >  1:3458168/1‑101 (MQ=255)
                        ttCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATc                                                                            <  2:2553710/101‑1 (MQ=255)
                         tCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATcc                                                                           <  2:3407447/101‑1 (MQ=255)
                         tCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATcc                                                                           <  1:3374497/101‑1 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATccc                                                                          <  1:2360148/101‑1 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATccc                                                                          >  2:3987989/1‑101 (MQ=255)
                           aTCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCcg                                                                         <  2:1009039/101‑1 (MQ=255)
                            tCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCcgc                                                                        >  1:25391/1‑101 (MQ=255)
                            tCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCcgc                                                                        <  2:1575226/101‑1 (MQ=255)
                              cAGGTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCgg                                                                      >  2:4022015/1‑101 (MQ=38)
                               aGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCggg                                                                     >  1:4230311/1‑101 (MQ=255)
                               aGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCggg                                                                     >  1:2947000/1‑101 (MQ=255)
                                 ttCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCggggg                                                                   >  1:267734/1‑101 (MQ=255)
                                  tCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCgggggc                                                                  >  1:4289844/1‑100 (MQ=255)
                                  tCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGt                                                                  >  1:1416966/1‑101 (MQ=255)
                                  tCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGt                                                                  >  1:611430/1‑101 (MQ=255)
                                  tCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATATTGCCCAAAATATCCCGCGGGGGt                                                                  <  1:4724337/101‑1 (MQ=38)
                                   ccgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGTc                                                                 <  1:896662/101‑1 (MQ=255)
                                   ccgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGTc                                                                 >  2:18366/1‑101 (MQ=255)
                                    cgcccgcGCGGCATCATCGCGCGGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGTcc                                                                <  1:4132153/101‑1 (MQ=38)
                                      cccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGTccgg                                                              <  1:3367403/101‑3 (MQ=255)
                                          cgcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGTCCGggggc                                                          <  2:2862514/101‑2 (MQ=37)
                                          cgcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGTCCGggggc                                                          <  1:2731677/101‑2 (MQ=37)
                                           gcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGgcg                                                        <  1:2794018/101‑1 (MQ=25)
                                               catcatCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGaa                                                    >  2:2767039/1‑101 (MQ=25)
                                               catcatCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGaa                                                    <  2:4676412/101‑1 (MQ=25)
                                               catcatCGCGCGGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGaa                                                    <  2:352818/101‑1 (MQ=25)
                                                atcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAg                                                   >  2:1368045/1‑101 (MQ=25)
                                                 tcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGc                                                  >  1:195771/1‑101 (MQ=25)
                                                 tcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGc                                                  >  2:221191/1‑101 (MQ=25)
                                                 tcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGc                                                  <  1:3243186/101‑1 (MQ=25)
                                                 tcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGc                                                  >  2:2440784/1‑101 (MQ=25)
                                                  catcGCGCTGCGCGGCCNCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGcc                                                 <  2:1148728/101‑1 (MQ=21)
                                                  catcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGcc                                                 <  1:3014417/101‑1 (MQ=25)
                                                  catcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGcc                                                 >  2:4754076/1‑101 (MQ=25)
                                                   atcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGccc                                                >  2:1883942/1‑101 (MQ=21)
                                                     cgcgCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGccccc                                              <  2:4239684/101‑1 (MQ=21)
                                                     cgcgCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGccccc                                              >  1:2050060/1‑101 (MQ=21)
                                                     cgcgCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGccccc                                              >  1:181476/1‑101 (MQ=21)
                                                         cggcccgcccgcgccgcCGGAGATGGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGgcg                                          <  2:3689995/93‑1 (MQ=17)
                                                         cTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGgcg                                          >  2:2891035/1‑101 (MQ=21)
                                                               ggCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCa                                    >  2:3206918/1‑101 (MQ=18)
                                                               gcccccgccgcGGGAGATTGAGCCGATCCCCCCCGCCCC.CC.TTTC.ATCTTGCCCAAAATATCCCGGGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCa                                    <  2:4545704/96‑1 (MQ=17)
                                                                        ccgggAGATTGAGCCGATCCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATACCCCGGGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGga                           <  1:1128929/98‑1 (MQ=18)
                                                                           ggAGATTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCGCCGGCGCCGCCAGATCCCGGAGAc                        >  2:2073043/1‑101 (MQ=18)
                                                                            gagaTTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGAcc                       >  2:4168766/1‑101 (MQ=18)
                                                                              gaTTGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGAcccc                     >  2:1638474/1‑101 (MQ=18)
                                                                                ttGAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGa                   >  1:2490759/1‑101 (MQ=18)
                                                                                 tgccccccccccccccC.CCG.CCTTTTC.CTCTTGCCAAAATTACCCGGCGGGGG.CCGGGGGGCGCGACGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGAt                  <  2:3730674/91‑1 (MQ=255)
                                                                                  gacccGATCCCCCCCG.CCGCCC.TTTC.ATCTTCCCAAAAATACCCCG.GGGGGGCCGGGGGGCGCGAACCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATg                 <  1:1249396/98‑1 (MQ=255)
                                                                                  gAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCTCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATg                 >  1:508158/1‑101 (MQ=18)
                                                                                  gAGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATACCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATg                 >  2:2428678/1‑101 (MQ=18)
                                                                                   aGCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGa                >  2:914175/1‑101 (MQ=18)
                                                                                    gCCGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGac               >  2:3541710/1‑101 (MQ=18)
                                                                                      cgaccccccccG.CCGCCC.TTTC.ATCTTCCCAAAAATACCCCGGGGGGGT.CGGGGGGCGCGAACCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGacac             <  1:1222536/96‑1 (MQ=255)
                                                                                      cGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGacac             >  2:4483/1‑101 (MQ=18)
                                                                                      cGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGacac             >  2:3758731/1‑101 (MQ=18)
                                                                                      cGATGCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGacac             <  1:1476167/101‑1 (MQ=18)
                                                                                          gCCCCCCG.CCGCCC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGcc         >  2:3139193/1‑101 (MQ=18)
                                                                                                  cccg.cC.CTTCAATCTTCCCAAAAATCCCCCG.GGGGGGCCGGGGGGCGCGACCCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGAcc  <  2:4567339/99‑1 (MQ=255)
                                                                                                      |                                                                                                   
TGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCG.GCG.CC.TTTC.ATCTTGCCCAAAATATCCCGCGGGGGTCCCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACC  >  NC_014034/1300960‑1301157

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG