Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 904,934 +GC 100% intergenic (+38/‑158) RCAP_rcc00821 → / → murB hypothetical protein/UDP‑N‑acetylmuramate dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034904,9341.→G100.0% 94.7 146intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (59/58):  ref base (3/26):  total (62/84)
Fisher's exact test strand distribution p-value = 5.23e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.66e-01
*NC_014034904,9342.→C100.0% 93.9 146intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (59/58):  ref base (3/26):  total (62/84)
Fisher's exact test strand distribution p-value = 5.23e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.14e-01

Only 100 of 159 total aligned reads displayed.

GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTT.TG.GCT..G..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAC  >  NC_014034/904834‑905024
                                                                                                        ||                                                                                              
  gCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cg                                                                                           <  1:4492416/101‑2 (MQ=255)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCGCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgcc                                                                                         <  1:4607222/101‑4 (MQ=255)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGGCTGCCCCCCTC.TTT.TG.GCT..G..cgcc                                                                                         >  2:4736587/1‑98 (MQ=255)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgcc                                                                                         <  2:4805324/101‑4 (MQ=255)
    tGAGCGTGCTGCGCGGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgcc                                                                                         <  2:3703741/101‑4 (MQ=255)
    ggggCGTGCTGCGCTGGCCGGTGCGG.GGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgcca                                                                                        <  1:4661159/98‑5 (MQ=37)
     gAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgcca                                                                                        <  2:2475222/101‑5 (MQ=255)
      aGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCC.CCTTT.TG.GCT..G..cgccaa                                                                                       >  1:4075995/1‑96 (MQ=38)
       gCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaa                                                                                      <  2:4133850/101‑7 (MQ=255)
        cGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaa                                                                                     <  1:1414180/101‑8 (MQ=37)
        cGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaa                                                                                     <  1:4369760/101‑8 (MQ=37)
         gggctggGCTGGCCGGGGCGGGGGCTCTTTGCC.GCGATGCCCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaat                                                                                    <  1:1012994/99‑9 (MQ=21)
         gTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgcccaaat                                                                                    >  2:4164367/1‑93 (MQ=37)
         gTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaat                                                                                    <  1:3861166/101‑9 (MQ=37)
         gTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaat                                                                                    <  1:3833559/101‑9 (MQ=37)
          tgctggGCTGGCCGGTGCGGCGGCTCTTTGCC.GCTATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaatt                                                                                   <  2:1124661/101‑10 (MQ=25)
          tgctgcGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGTGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaatt                                                                                   <  2:630250/101‑10 (MQ=25)
          tgctgcGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaatt                                                                                   <  1:2923359/101‑10 (MQ=25)
          tgctgcGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaatt                                                                                   <  1:3152759/101‑10 (MQ=25)
          tgctgcGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..G..cgccaaaatt                                                                                   <  1:329164/101‑10 (MQ=25)
           gctgcGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTc                                                                                  <  2:229295/101‑1 (MQ=25)
            ctgcGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCa                                                                                 <  1:1935753/101‑1 (MQ=25)
             tgcGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCAc                                                                                <  2:4602988/101‑1 (MQ=25)
             tgcGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCAc                                                                                <  2:4986527/101‑1 (MQ=25)
              gcgcTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CCAAATTCAcc                                                                               >  2:3286486/1‑101 (MQ=25)
              gcgcTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCAcc                                                                               <  2:1490501/101‑1 (MQ=25)
               cgcTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCAccc                                                                              <  2:4923352/101‑1 (MQ=25)
                tcgggCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCAcccc                                                                             <  2:4382248/98‑1 (MQ=21)
                gcTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCAcccc                                                                             <  2:2939298/101‑1 (MQ=25)
                  tGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCAcccccc                                                                           >  2:2111291/1‑101 (MQ=21)
                  tGGCCGGTGCGGCGGCTCTTGGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCAcccccc                                                                           <  1:991909/101‑1 (MQ=21)
                   ggCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCg                                                                          >  1:3255591/1‑101 (MQ=21)
                   ggCCGGTGCGGCCGCTCTTTGCC.GCGTTGCCCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCg                                                                          <  1:2727414/101‑1 (MQ=18)
                    gCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTGTC.TGC..G..C..CAAAATTCACCCCCCgg                                                                         <  2:3840833/101‑1 (MQ=21)
                      cGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGTGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCg                                                                       <  2:2694125/101‑1 (MQ=21)
                       ggTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCCGCgg                                                                      >  1:1016686/1‑100 (MQ=21)
                       cgtgggggcgcTCTTTGCCGGCGTGGCCCCGGCGCGCCGGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTT.GGCGGC..G..C..CAAAATTC.CCCCCCGGCGt                                                                      <  2:1697909/92‑1 (MQ=9)
                         tGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTAt                                                                    <  1:3750859/101‑1 (MQ=21)
                          gcggcgGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCggaccttt                                                                   >  1:2075033/1‑95 (MQ=18)
                          gcggcgGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTCtt                                                                   >  1:292082/1‑101 (MQ=21)
                           cggcggCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTAttt                                                                  >  1:936372/1‑101 (MQ=21)
                           cggcggCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATACACCCCCCGGCGTAttt                                                                  >  2:4539227/1‑101 (MQ=21)
                            ggcggcTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTg                                                                 >  1:1934338/1‑101 (MQ=21)
                            ggcggcTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTg                                                                 >  1:3129949/1‑101 (MQ=21)
                            ggcggcTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTg                                                                 >  1:3079952/1‑101 (MQ=21)
                             gcggcTCTTTGCC.TCGTTGCCCCGGCGCGAGTGAGGGCTGAGAATGGTGGCGCTTCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGc                                                                <  1:130820/101‑1 (MQ=12)
                               gtctctTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCAc                                                              <  2:3548813/99‑1 (MQ=18)
                                 ctctTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCa                                                            >  2:3424034/1‑101 (MQ=21)
                                   ctTTGCC.GCGATGCACCGGAGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAgg                                                          <  2:3527426/101‑1 (MQ=21)
                                     ttGCC.GCGATGTACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGaa                                                        <  1:484115/101‑1 (MQ=18)
                                     ttGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGaa                                                        <  2:3649886/101‑1 (MQ=21)
                                       gCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGaaaa                                                      <  1:4615926/101‑1 (MQ=21)
                                           gggATGCACCGGCGCGACTGAGGGCCGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGc                                                   <  1:1671118/99‑1 (MQ=18)
                                             gATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGccc                                                 >  1:2650606/1‑101 (MQ=21)
                                              gtggaCCGGGGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCg                                                <  2:3893015/97‑1 (MQ=18)
                                              aTGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCg                                                >  2:3868711/1‑101 (MQ=21)
                                                  aCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGt                                            >  2:192098/1‑101 (MQ=21)
                                                   ccGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGtt                                           <  1:3174275/101‑1 (MQ=21)
                                                    cGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGttt                                          >  2:212519/1‑101 (MQ=21)
                                                     ggCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCCGCGTATTTGCACCCGGAAAAAACCCGGCGtttt                                         >  1:4125215/1‑101 (MQ=18)
                                                      gcgcGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGttttt                                        <  1:1131153/101‑1 (MQ=21)
                                                       cgcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTc                                       >  2:690568/1‑101 (MQ=21)
                                                       cgcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTc                                       <  1:4011803/101‑1 (MQ=21)
                                                         ggACTGATGGCTGGGCCGGGGGGCCCTTCCCCCCTC.TTT.TG.GCT..GCTC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCt                                     <  1:214754/100‑1 (MQ=21)
                                                            cTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGt                                  >  2:4709159/1‑101 (MQ=21)
                                                              gAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGCTTTTCCTGGTCg                                >  1:4813507/1‑101 (MQ=18)
                                                               aGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGa                               <  2:334652/101‑1 (MQ=21)
                                                                gggCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGaa                              <  2:4001336/101‑1 (MQ=21)
                                                                  gCGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAAt                            <  1:1144777/101‑1 (MQ=21)
                                                                   cGGGGCAGGGTGGCGCTGCCCCCC.CCTTT.TG.GCT..G..CGGCAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGTCGTTTTTCCTGGTCGAGATa                           <  2:242754/101‑1 (MQ=21)
                                                                   cGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATa                           <  1:4395998/101‑1 (MQ=21)
                                                                   cGGGGCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..GCCC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCGGGGGGTTTTCCTGGTCGAAATa                           >  2:527156/1‑101 (MQ=18)
                                                                   cGGGGCAGGGGAGCTCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCCCCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATa                           <  2:2803203/101‑1 (MQ=21)
                                                                    ggggCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATAc                          >  2:1080786/1‑101 (MQ=21)
                                                                     gggCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCT..GATC..CGAAATTCACCCCCCGGCGTACTTGCACCAGGAAAAAGCCCGGCGTTTTTCATGGTCGAAATACg                         <  1:4321433/101‑1 (MQ=25)
                                                                       gCAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGcc                       <  1:2344656/101‑1 (MQ=21)
                                                                        cAGGGGGGCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCg                      >  1:744000/1‑101 (MQ=21)
                                                                           gggggCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTCCCTGGTCGAAATACGCCgggg                   <  1:3717463/101‑1 (MQ=21)
                                                                             gggCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgggggg                 <  2:4341886/101‑1 (MQ=25)
                                                                              ggCGCTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGa                <  2:292029/101‑1 (MQ=25)
                                                                              ggCGCTGCCCACCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGa                <  1:848384/101‑1 (MQ=25)
                                                                                cgcTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAgc              >  1:4896555/1‑100 (MQ=25)
                                                                                 gcTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTc             >  2:2706114/1‑101 (MQ=37)
                                                                                  cTGCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc            <  1:2168750/101‑1 (MQ=37)
                                                                                    gCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGc          >  2:1788048/1‑101 (MQ=37)
                                                                                    gCCCCCCTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTNCCTGGTCGAAATACGCCGGGGGGAGTCCGc          >  1:574146/1‑101 (MQ=37)
                                                                                     ccccccTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCAGCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt         <  1:2197539/101‑1 (MQ=37)
                                                                                      cccccTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt        <  1:5016027/101‑1 (MQ=37)
                                                                                      cccccT.TTTT.TG.GCTGCG..C..CAAAATTCCCCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt        <  1:3644587/101‑1 (MQ=25)
                                                                                       ccccTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTg       >  2:1226729/1‑101 (MQ=37)
                                                                                        cccTC.TTT.TG.GCTGCG..C..CAAAATTCCCCCCCCGGCGTATTTGCCCCAGGAAAAAGCCCGGCTTTTTTCCTGGGCGAAATACGCCGGGGGGAGTCCGCTTGc      <  1:1824105/101‑1 (MQ=25)
                                                                                        cccTC.TTT.TG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGc      <  1:1315122/101‑1 (MQ=37)
                                                                                          cctcttt.tG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCg     >  2:1759580/5‑101 (MQ=38)
                                                                                           ctcttt.tG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCgg    <  2:1559695/98‑1 (MQ=38)
                                                                                           ctcttt.tG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCgg    >  1:4397741/4‑101 (MQ=37)
                                                                                           ctcttt.tG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCgg    <  1:4519426/98‑1 (MQ=38)
                                                                                            ttttt.tG.GCT..GCCC..CAAAATTCCCCCCCCGGCGTATTTGCACCAGAAAAAAGCCCGGCTTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGa   <  1:3424298/100‑1 (MQ=25)
                                                                                            tcttt.tG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGa   >  2:1914444/3‑101 (MQ=37)
                                                                                             cttt.tG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc  >  1:4219990/2‑101 (MQ=255)
                                                                                             cttt.tG.GCTGCG..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc  >  1:2204341/2‑101 (MQ=255)
                                                                                                        ||                                                                                              
GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCC.GCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTT.TG.GCT..G..C..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAC  >  NC_014034/904834‑905024

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 39 ≤  ATCG