Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,434,412 C→G 20.9% intergenic (‑191/‑767) RCAP_rcc01334 ← / → RCAP_rcc01335 hypothetical protein/Fis family GAF modulated sigma54 specific transcriptional regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,434,4120C→G20.9% 271.5 ‑0.8 0.0521intergenic (‑191/‑767)RCAP_rcc01334/RCAP_rcc01335hypothetical protein/Fis family GAF modulated sigma54 specific transcriptional regulator
Reads supporting (aligned to +/- strand):  new base (50/58):  ref base (171/236):  total (224/297)
Fisher's exact test strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 535 total aligned reads displayed.

AAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAA  >  NC_014034/1434312‑1434510
                                                                                                       |                                                                                                     
  aCAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGGCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCCCG.GCTGCGTCAACCtt                                                                                                     >  1:6568961/1‑101 (MQ=255)
   cAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..g                                                                                                  <  1:1976671/101‑1 (MQ=255)
    aGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..Gt                                                                                                 <  1:4530900/101‑1 (MQ=255)
        tgctTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCt                                                                                            <  2:4574614/101‑1 (MQ=255)
         gctTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTc                                                                                           >  2:7051810/1‑101 (MQ=255)
              gggCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTcccccc                                                                                      <  1:5766408/101‑1 (MQ=255)
               ggCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCcg                                                                                     <  2:3967990/101‑1 (MQ=255)
                 cGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCcccg                                                                                   >  2:1318303/1‑101 (MQ=255)
                   gAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCggg                                                                                 >  1:9083024/1‑101 (MQ=255)
                      cacCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGc                                                                              <  1:1370762/101‑1 (MQ=255)
                       acCAGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCggggggg                                                                             >  1:7213238/1‑99 (MQ=255)
                            aTCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGAt                                                                        >  2:1912120/1‑101 (MQ=255)
                             tCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATc                                                                       <  1:5893786/101‑1 (MQ=255)
                             cgggctatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTgg                                                                          <  1:1142919/92‑1 (MQ=2)
                             cgggctatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTgg                                                                          <  2:3512227/92‑1 (MQ=2)
                             cgggctatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTgg                                                                          <  2:2351921/92‑1 (MQ=2)
                              ccTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCg                                                                      <  1:4335561/101‑1 (MQ=255)
                              ccTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCg                                                                      <  1:3858817/101‑1 (MQ=255)
                              ccTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCg                                                                      <  2:4788773/101‑1 (MQ=255)
                               ggctatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGAt                                                                        <  2:6067518/94‑1 (MQ=2)
                               ggctatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGAt                                                                        >  2:4926267/8‑101 (MQ=2)
                               ggctatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGAt                                                                        <  1:5524151/94‑1 (MQ=2)
                               ctgctgGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGa                                                                     >  2:9182526/1‑101 (MQ=255)
                               ctgctgGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCCCGGGGGGCTGGATCGa                                                                     >  2:5422294/1‑101 (MQ=255)
                                gctatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGGTc                                                                       >  1:4488511/7‑101 (MQ=2)
                                 ctatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCg                                                                      <  1:178397/96‑1 (MQ=2)
                                  tatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGa                                                                     <  1:8486400/97‑1 (MQ=2)
                                  gatcttCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                                  <  1:8201281/96‑1 (MQ=255)
                                  ctgGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                                  <  1:1592449/101‑1 (MQ=255)
                                  ctgGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                                  <  1:4761698/101‑1 (MQ=255)
                                   tgGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAaca                                                                 >  1:3307249/1‑101 (MQ=255)
                                   tgGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAaca                                                                 <  1:2631443/101‑1 (MQ=255)
                                   atcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                    >  1:2309858/4‑101 (MQ=2)
                                    tcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                   <  2:184020/99‑1 (MQ=2)
                                    gggTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacac                                                                <  1:6925140/101‑1 (MQ=255)
                                    gggTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacac                                                                <  2:5156404/101‑1 (MQ=255)
                                     cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                                  >  2:6349877/2‑101 (MQ=2)
                                     cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGGGGATCGAAac                                                                  >  2:2982106/2‑101 (MQ=2)
                                     cgTCGCGCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                                  <  2:6537707/100‑1 (MQ=2)
                                       tCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGa                                                             >  2:6131633/1‑101 (MQ=255)
                                         gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGAtt                                                           <  2:8781499/101‑1 (MQ=255)
                                         gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGAtt                                                           <  2:5664748/101‑1 (MQ=255)
                                         gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGAtt                                                           >  2:79879/1‑101 (MQ=255)
                                          cTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgc                                                             >  2:4191838/1‑96 (MQ=2)
                                          cTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGCTTg                                                          >  2:7291494/1‑101 (MQ=255)
                                          cTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTg                                                          >  1:7614744/1‑101 (MQ=255)
                                          cTCCCCCAGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgc                                                             >  2:3633396/1‑96 (MQ=2)
                                           tCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgct                                                            >  1:5266945/1‑95 (MQ=2)
                                           tCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTgt                                                         >  2:3157882/1‑100 (MQ=255)
                                           tCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGc                                                         >  1:4913538/1‑101 (MQ=255)
                                           tCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGc                                                         >  2:3283867/1‑101 (MQ=255)
                                            ccccccGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgctg                                                           >  2:2698435/1‑94 (MQ=2)
                                            ccccccGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgctg                                                           <  2:5301059/101‑8 (MQ=2)
                                            ccccccGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgctg                                                           <  1:3490039/101‑8 (MQ=2)
                                             cccccGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgctga                                                          <  1:6024221/101‑9 (MQ=2)
                                             cccccGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCAt                                                       <  2:8590553/101‑1 (MQ=255)
                                              ccccGCGGGGGCGTGGATCGAAAC..TTGTCTGACGCGCTCATCA.GC.ACGTCAAGATTCAGT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACTCggcaaggac                                                       >  2:720983/1‑94 (MQ=37)
                                              ccccGCGGGGGCGTGGATCGAAAC..TTGTCTGACGCGCTCATCA.GC.ACGTCAAGATTCAGT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACTCggcaacgac                                                       <  1:8791791/101‑8 (MQ=37)
                                               cccGCGGGGGCGTTGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATgg                                                     <  2:7853303/101‑1 (MQ=255)
                                               cccGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGC..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgctgcca                                                        >  2:230445/1‑91 (MQ=0)
                                                   cgGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGcc                                                 <  1:3702824/101‑1 (MQ=255)
                                                   cgGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGcc                                                 >  2:5806103/1‑101 (MQ=255)
                                                    gggggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCg                                                >  2:1055954/1‑101 (MQ=255)
                                                    gggggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCg                                                >  1:5147426/1‑101 (MQ=255)
                                                     ggggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACGCGGATTGCATGGCGCCGa                                               >  1:6195534/1‑101 (MQ=255)
                                                      gggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAg                                              >  1:7487776/1‑101 (MQ=255)
                                                       ggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.TGCTCCCCCCGCGGGGGCGTCGATCGAAACACGGATTGCATGGCGCCGAgg                                             >  2:1599192/1‑101 (MQ=255)
                                                           tGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACa                                         <  1:3171976/101‑1 (MQ=255)
                                                           tGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACa                                         <  2:6752850/101‑1 (MQ=255)
                                                            ggATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAg                                        >  2:4482443/1‑101 (MQ=255)
                                                                cGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCtt                                    <  2:9048746/101‑1 (MQ=255)
                                                                cGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCtt                                    <  2:287973/101‑1 (MQ=255)
                                                                 gAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCttt                                   <  2:3895339/101‑1 (MQ=255)
                                                                 gAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCttt                                   >  2:1476317/1‑101 (MQ=255)
                                                                 gAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCttt                                   >  1:6301726/1‑101 (MQ=255)
                                                                 gAAAC..GGCTCCGCTCTTCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCttt                                   <  1:8696868/101‑1 (MQ=255)
                                                                 gAAAC..GGCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCttt                                   <  2:1110265/101‑1 (MQ=255)
                                                                  aaaC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTg                                  <  1:79879/101‑1 (MQ=255)
                                                                  aaaC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTg                                  <  2:6924661/101‑1 (MQ=255)
                                                                  aaaC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTg                                  >  1:6346989/1‑101 (MQ=255)
                                                                   aaC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGc                                 >  1:7878248/1‑101 (MQ=255)
                                                                    aC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCa                                <  1:7062775/101‑1 (MQ=255)
                                                                    aC..GACTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCATTGcg                                >  2:823431/1‑100 (MQ=255)
                                                                        gTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCtttagatc                              >  2:84865/1‑96 (MQ=255)
                                                                         tctcCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTc                             <  1:4099684/101‑1 (MQ=255)
                                                                         tctcCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTc                             <  2:6464249/101‑1 (MQ=255)
                                                                             cGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTc                         <  2:8890638/101‑1 (MQ=255)
                                                                               ctctGCTTCACG.GCTGCGTCAACCTT..GTCCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTcc                        <  1:8941057/101‑1 (MQ=255)
                                                                               ctctGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTccc                       <  1:7689059/101‑1 (MQ=255)
                                                                                      tCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCgg                >  2:7638701/1‑101 (MQ=255)
                                                                                         cG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCggggg             >  2:2142165/1‑101 (MQ=255)
                                                                                          g.gCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGc            >  1:2090193/1‑101 (MQ=40)
                                                                                          g.gCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGc            <  1:1055972/101‑1 (MQ=40)
                                                                                            gCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCg           >  2:173650/1‑101 (MQ=40)
                                                                                            gCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCg           >  2:220581/1‑101 (MQ=40)
                                                                                                cGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGCCAGGCTTGGCAGTCGCTCCCCCCGCGGGGGCGTGGa       >  1:2963334/1‑101 (MQ=255)
                                                                                                  tCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATc     <  2:7575545/101‑1 (MQ=33)
                                                                                                  tCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATc     <  1:8314662/101‑1 (MQ=33)
                                                                                                    aaCCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGCAGATcct   <  2:6669412/101‑3 (MQ=9)
                                                                                                     aCCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa  >  1:5663519/1‑101 (MQ=11)
                                                                                                       |                                                                                                     
AAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTT..GT.CGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAA  >  NC_014034/1434312‑1434510

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG