Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,157,926 2 bp→AC 100% intergenic (‑298/‑483) RCAP_rcc02004 ← / → RCAP_rcc02005 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,157,9260C→A100.0% 57.1 255intergenic (‑298/‑484)RCAP_rcc02004/RCAP_rcc02005hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (106/126):  ref base (20/2):  total (127/128)
Fisher's exact test strand distribution p-value = 2.95e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.58e-01
*NC_0140342,157,9270G→C100.0% 49.5 255intergenic (‑299/‑483)RCAP_rcc02004/RCAP_rcc02005hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (106/126):  ref base (20/2):  total (127/128)
Fisher's exact test strand distribution p-value = 2.95e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.77e-01

Only 100 of 263 total aligned reads displayed.

GAGGCTCTCCCAGCATACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.ACGGCCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTCCGCCGTCGTC  >  NC_014034/2157828‑2158021
                                                                                                   ||                                                                                              
      ctcCCAGCATACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCg                                                                                         >  1:584131/1‑101 (MQ=255)
        cccAGCATACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGt                                                                                       >  1:1794879/1‑101 (MQ=255)
         ccAGCATACCAAACCGCTGGCGCCCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGtt                                                                                      <  2:388138/101‑1 (MQ=38)
         ccAGCATACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGtt                                                                                      <  1:2387446/101‑1 (MQ=255)
           aGCATACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCa                                                                                    >  1:3175/1‑101 (MQ=255)
            gCATACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAg                                                                                   >  2:1035207/1‑101 (MQ=255)
             cATACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGt                                                                                  <  2:3037953/101‑1 (MQ=255)
              aTACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTc                                                                                 >  1:1550735/1‑101 (MQ=255)
               tACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATCAACGG.CTTCGGTTCAGTcc                                                                                <  2:2696500/101‑1 (MQ=255)
               tACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTcc                                                                                >  2:4594511/1‑101 (MQ=255)
               tACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTcc                                                                                <  2:364154/101‑1 (MQ=255)
               tACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTcc                                                                                <  2:2081763/101‑1 (MQ=255)
               tACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTcc                                                                                >  2:1636430/1‑101 (MQ=255)
               tACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTcc                                                                                <  1:707096/101‑1 (MQ=255)
                  cAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTccgcc                                                                             <  1:3135474/101‑1 (MQ=255)
                   aaaCCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGGCCGATGAAAATC.AACGGCTTCGGTTCAGTtcggcg                                                                            >  2:3627235/1‑101 (MQ=37)
                   aaaCCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTccgccg                                                                            >  2:2368366/1‑101 (MQ=255)
                     aCCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCcgtc                                                                          <  2:3586597/101‑1 (MQ=255)
                     aCCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCcgtc                                                                          >  1:3384834/1‑101 (MQ=255)
                     aCCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCcgtc                                                                          <  2:918308/101‑1 (MQ=255)
                        gCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCgtcgtc                                                                       <  1:3211869/101‑1 (MQ=255)
                        gCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCAGCcgtcgta                                                                       >  1:3294048/1‑100 (MQ=37)
                          tGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCtt                                                                     >  2:3075053/1‑101 (MQ=255)
                             cGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGc                                                                  <  1:8483/101‑1 (MQ=255)
                               aCCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCgg                                                                <  2:3276619/101‑1 (MQ=255)
                                 cATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCgggg                                                              >  1:1356119/1‑99 (MQ=255)
                                  aTCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGa                                                             >  2:3210507/1‑101 (MQ=255)
                                  aTCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGa                                                             <  2:124944/101‑1 (MQ=255)
                                   tCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGAc                                                            <  2:3384628/101‑1 (MQ=255)
                                   tCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGAc                                                            <  1:1973135/101‑1 (MQ=255)
                                      ggTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTg                                                         <  2:3537068/101‑1 (MQ=255)
                                      ggTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTg                                                         >  2:1474464/1‑101 (MQ=255)
                                       gTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGa                                                        <  2:467666/101‑1 (MQ=255)
                                        ttCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAg                                                       >  1:4565517/1‑101 (MQ=255)
                                        ttCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAg                                                       >  1:2090461/1‑101 (MQ=255)
                                        ttCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAg                                                       >  2:4436676/1‑101 (MQ=255)
                                         tCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAgg                                                      >  2:1719394/1‑101 (MQ=255)
                                         tCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AAAGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGGGACGTGAgg                                                      >  2:1857036/1‑101 (MQ=38)
                                          cAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGc                                                     <  1:4031451/101‑1 (MQ=255)
                                          cAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGc                                                     >  2:469496/1‑101 (MQ=255)
                                           aGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCa                                                    <  1:2940995/101‑1 (MQ=255)
                                           aGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCa                                                    <  1:4596668/101‑1 (MQ=255)
                                           aGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGCTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCa                                                    >  1:3903597/1‑101 (MQ=38)
                                             tCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCAcc                                                  <  1:2913987/101‑1 (MQ=255)
                                              ccgccgTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTGGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCGCCt                                                 <  2:4055347/101‑1 (MQ=37)
                                              ccgccgTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCt                                                 >  2:2173284/1‑101 (MQ=255)
                                              ccgccgTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCt                                                 >  1:1424522/1‑101 (MQ=255)
                                              ccgccgTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCt                                                 <  2:3932315/101‑1 (MQ=255)
                                                gccgTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTaa                                               >  1:1197173/1‑101 (MQ=255)
                                                gccgTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCATCTTGGCGGTGACGTGAGGCACCTaa                                               <  2:1166733/101‑1 (MQ=38)
                                                 ccgTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTaaa                                              <  1:2655455/101‑1 (MQ=255)
                                                  cgtcgtCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAAc                                             >  1:4153160/1‑101 (MQ=255)
                                                  cgtcgtCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAAc                                             >  1:3407287/1‑101 (MQ=255)
                                                     cgtcTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGGGACGTGAGGCACCTAAACAag                                          >  2:1195090/1‑101 (MQ=38)
                                                       tcTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAagag                                        >  1:4373319/1‑101 (MQ=255)
                                                         ttGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGt                                      >  1:2262507/1‑101 (MQ=255)
                                                          tGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGtt                                     <  1:1007620/101‑1 (MQ=255)
                                                            gCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCt                                   <  1:3738166/101‑1 (MQ=255)
                                                            gCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCt                                   <  1:1652795/101‑1 (MQ=255)
                                                             cGGTGATGTGATGCACAGGAGAGACCGATGAAAATCAACGG.CTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTa                                  >  2:673269/1‑101 (MQ=255)
                                                             cGGGGATGTGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTa                                  >  2:2847999/1‑101 (MQ=38)
                                                              ggTGATGGGATGCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTAc                                 <  2:2882289/101‑1 (MQ=38)
                                                               gtgatgtgatGCACAGGAGAGACCGATGAAAATCAACGG.CTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACg                                >  1:4492943/1‑101 (MQ=255)
                                                                tgatgtgatgCACAGGAGAGACCGATGAAAATCAACGG.CTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGa                               >  1:3607047/1‑101 (MQ=255)
                                                                tgatgtgatgCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGa                               >  2:4462330/1‑101 (MQ=255)
                                                                tgatgtgatgCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGa                               <  2:1323923/101‑1 (MQ=255)
                                                                tgatgtgatgCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGa                               <  2:3401745/101‑1 (MQ=255)
                                                                 gatgtgatgCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGAt                              <  2:2020888/101‑1 (MQ=255)
                                                                 gatgtgatgCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGAt                              >  2:1466121/1‑101 (MQ=255)
                                                                  atgtgatgCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATg                             <  1:3844970/101‑1 (MQ=255)
                                                                   tgtgatgCACAGGAGAGACCGATGAAAATCAACGG.CTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGc                            >  2:2133296/1‑101 (MQ=255)
                                                                    gtgatgCACAGGAGAGACCGATGAAAATCAACGG.CTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGcc                           >  2:384594/1‑101 (MQ=255)
                                                                     tgatgCACAGGAGAGACCGATGAAAATCAACGG.CTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCg                          >  2:3450006/1‑101 (MQ=255)
                                                                     tgatgCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCg                          <  1:2343579/101‑1 (MQ=255)
                                                                      gatgCACAGGAGAGACCGATGAAAATT.AACGGCTTTGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGATGCTCCTAAACAAGAGGTTCTACGATGCCGa                         <  2:764042/101‑1 (MQ=37)
                                                                      gatgCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGa                         <  1:1024460/101‑1 (MQ=255)
                                                                       atgCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGaa                        <  1:3575801/101‑1 (MQ=255)
                                                                         gCACAGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAg                      <  1:431034/101‑1 (MQ=255)
                                                                         ctacaGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAg                      <  1:1131751/99‑1 (MQ=255)
                                                                           acaGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAggg                    <  1:3444980/101‑1 (MQ=255)
                                                                             aGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGt                  <  2:2433411/101‑1 (MQ=255)
                                                                             aGGAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGt                  >  1:244081/1‑101 (MQ=255)
                                                                              ggAGAGACCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGtt                 >  1:3501192/1‑101 (MQ=255)
                                                                              ggAGAGACCGATGAAAATC.AACGGATTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGtt                 <  2:3722029/101‑1 (MQ=38)
                                                                               gagagaCCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTc                <  1:4390738/101‑1 (MQ=255)
                                                                                agagaCCGATGAAAATCAACGG.CTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCa               >  2:4206505/1‑101 (MQ=255)
                                                                                agagaCCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCa               <  2:1550634/101‑1 (MQ=255)
                                                                                  agaCCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGt             >  1:1080195/1‑101 (MQ=255)
                                                                                  agaCCGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGt             <  1:3075241/101‑1 (MQ=255)
                                                                                     ccGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTccg          <  2:2990221/101‑1 (MQ=255)
                                                                                      cGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTccgc         >  1:792757/1‑101 (MQ=255)
                                                                                      cGATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTAAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTccgc         >  2:3158161/1‑101 (MQ=38)
                                                                                       gATGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTccgcc        >  1:1292643/1‑101 (MQ=255)
                                                                                        aTGAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTccgccg       >  1:4584009/1‑101 (MQ=255)
                                                                                          gAAAATC.AACGGTTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTCCGCcgtc     <  2:2625855/101‑1 (MQ=38)
                                                                                          gAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTCCGCcgtc     >  1:1573003/1‑101 (MQ=255)
                                                                                          gAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTCCGCcgtc     >  1:2682363/1‑101 (MQ=255)
                                                                                          gAAAATC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTCCGCcgtc     >  1:363093/1‑101 (MQ=255)
                                                                                            aaaTC.AACGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTCCGCcgtcgt   >  2:4754624/1‑101 (MQ=255)
                                                                                             aaTC.AAGGGCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTTAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTCCGCCgtcgtc  <  2:1337696/101‑1 (MQ=255)
                                                                                                   ||                                                                                              
GAGGCTCTCCCAGCATACCAAACCGCTGGCGACCATCAGGTTCAGTCCGCCGTCGTCTTGGCGGTGATGTGATGCACAGGAGAGACCGATGAAAATC.ACGGCCTTCGGTTCAGTCCGCCGTCGTCTTGGCGGTGACGTGAGGCACCTAAACAAGAGGTTCTACGATGCCGAAAGGGGTTCAGTCCGCCGTCGTC  >  NC_014034/2157828‑2158021

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG