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% 25.8 125intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (70/43):  ref base (1/9):  total (71/54)
Fisher's exact test strand distribution p-value = 1.85e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.60e-01
*NC_014034904,9342.→C100.0% 31.7 125intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (69/43):  ref base (1/11):  total (71/54)
Fisher's exact test strand distribution p-value = 4.42e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.11e-01

Only 100 of 132 total aligned reads displayed.

GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTGGCT..GC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACG  >  NC_014034/904834‑905025
                                                                                                     ||                                                                                              
gggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.CTTTGGCT..g                                                                                               >  1:1725835/1‑101 (MQ=255)
   cTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..Gcgc                                                                                            <  1:1519883/101‑3 (MQ=255)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGGCCCCCTC.TTTTGGCT..Gcgcc                                                                                           >  2:2663949/1‑98 (MQ=255)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..Gcgcc                                                                                           >  2:3799087/1‑98 (MQ=255)
        cGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..Gcgccaaaa                                                                                       <  1:543397/101‑8 (MQ=37)
         gTGCTGCGGTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCCCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..Gcgccaaaat                                                                                      <  1:3066743/101‑9 (MQ=25)
         gTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..Gcgccaaaat                                                                                      <  1:2864022/101‑9 (MQ=37)
         gTGCGGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..Gcgccaaaat                                                                                      <  1:3475124/101‑9 (MQ=25)
            ctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCa                                                                                   >  1:2996326/1‑101 (MQ=25)
             tgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCAc                                                                                  <  1:4189093/101‑1 (MQ=25)
             ggcgcTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCAc                                                                                  <  1:1592884/100‑1 (MQ=25)
              gcgcTGGACGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCAcc                                                                                 >  1:1330803/1‑101 (MQ=21)
                 cTGGCCGGTGCGGTGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCAccccc                                                                              >  1:3649430/1‑101 (MQ=21)
                 cTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CCAAATTCAccccc                                                                              >  2:438510/1‑101 (MQ=21)
                 cTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCAccccc                                                                              >  2:116562/1‑101 (MQ=21)
                  tGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCAcccccc                                                                             >  2:2118434/1‑101 (MQ=21)
                    gCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.gg                                                                          >  1:1830108/1‑101 (MQ=21)
                     ccGGTGCGGCGGCTCTTTTCCGCGATGCTCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTGGCTTC.GC..CAAAATTCACCCCCC.GGc                                                                         <  2:422682/101‑1 (MQ=18)
                     ccGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGc                                                                         <  2:2845610/101‑1 (MQ=21)
                        gTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTa                                                                      >  2:2064750/1‑101 (MQ=21)
                        gTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTa                                                                      <  2:1854791/101‑1 (MQ=21)
                           cggcggCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTAttt                                                                   >  1:2554385/1‑101 (MQ=21)
                                gcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCAcc                                                              >  1:55546/1‑101 (MQ=21)
                                gcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCAcc                                                              >  2:1479203/1‑101 (MQ=21)
                                  tctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAg                                                            >  2:4253125/1‑101 (MQ=21)
                                   ctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAgg                                                           <  1:4432635/101‑1 (MQ=21)
                                   ctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAgg                                                           <  1:2501909/101‑1 (MQ=21)
                                     ttGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGaa                                                         >  1:1605012/1‑101 (MQ=21)
                                     ttGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGaa                                                         <  2:2882401/101‑1 (MQ=21)
                                      tGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGaaa                                                        >  2:4529002/1‑101 (MQ=21)
                                        ccGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGGCTGACCCCCTCTTTGGGCTTC.GC..CAAAATTCACCCCCC.GGCGATTTTGCACCAGGagaaa                                                      >  1:3921020/1‑101 (MQ=17)
                                        ccGCCATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGaaaaa                                                      <  1:3664159/101‑1 (MQ=18)
                                           cgATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGcc                                                   >  2:2801335/1‑101 (MQ=21)
                                           cgATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGcc                                                   >  2:3971655/1‑101 (MQ=21)
                                            gATGCCCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGccc                                                  <  1:3283003/101‑1 (MQ=21)
                                             aTGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCg                                                 >  1:2693186/1‑101 (MQ=21)
                                              tGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCgg                                                <  2:140051/101‑1 (MQ=21)
                                              tGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCgg                                                >  1:3447674/1‑101 (MQ=21)
                                              tGCACCGGCGCGACTGAGGGCGGGGCAGAGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCgg                                                >  1:2438607/1‑101 (MQ=21)
                                                 cccGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGt                                             <  2:374482/100‑1 (MQ=21)
                                                      cgcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTc                                        <  2:4456483/101‑1 (MQ=21)
                                                       gcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTcc                                       >  2:147045/1‑101 (MQ=21)
                                                       gcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTcc                                       >  2:1756770/1‑101 (MQ=21)
                                                        cgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.CGCGCATTTGCACCACGAAAAAGCCCGGCGTTTTTCCt                                      >  1:7330/1‑101 (MQ=18)
                                                        cgACTGAGGGCGGCGCAGGAGGGCGCTGCCTCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCCTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCt                                      <  1:440122/101‑1 (MQ=18)
                                                         gACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCACGGCGTTCTTCATg                                     >  2:433502/1‑101 (MQ=18)
                                                         gACTGAGGGCGGGGCA.GGGGGCGCGGCCCCCCCCCTTTGGCCTC.GC..CAAAAATTCCCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTgg                                    <  2:2221254/101‑1 (MQ=25)
                                                          aCTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGGAAAACGCCGGCGTTTTTCCTgg                                    >  2:740969/1‑101 (MQ=18)
                                                           cttaGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGt                                   >  1:4615515/4‑101 (MQ=18)
                                                            tGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGGGC..CAAAATTCACCCCCC.GGCGTATTTGCACCACGAAAAAGCCCGGCGTTTTTCCTGGTc                                  >  2:1522456/1‑101 (MQ=21)
                                                             gAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCg                                 >  1:1305927/1‑101 (MQ=21)
                                                               gggCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.CGCGCCTTTGCACCAGGAAAAAGCCCGGCGCTTTTCCTGGTCGaa                               >  1:831130/1‑101 (MQ=18)
                                                               gggCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..GCGGCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGaa                               <  2:1757080/101‑1 (MQ=21)
                                                                       cAGGGGGGCGCTTCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCg                       <  2:359727/101‑1 (MQ=21)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                      >  1:1769828/1‑101 (MQ=25)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                      >  2:2368886/1‑101 (MQ=25)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                      <  2:3233534/101‑1 (MQ=25)
                                                                         ggggggCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCggg                     >  1:3325459/1‑101 (MQ=25)
                                                                         ggggggCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCggg                     >  2:4305059/1‑101 (MQ=25)
                                                                          gggggCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgggg                    <  1:3349761/101‑1 (MQ=25)
                                                                           ggggCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCggggg                   >  2:2125901/1‑101 (MQ=25)
                                                                           ggggCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACCCCggtcg                   >  1:2552987/1‑98 (MQ=21)
                                                                             ggCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGa                 >  1:3763516/1‑101 (MQ=25)
                                                                             ggCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGa                 <  1:3837191/101‑1 (MQ=25)
                                                                             ggCGCTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGa                 <  1:3868298/101‑1 (MQ=25)
                                                                              gcccTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTGCCCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAg                <  1:4032029/98‑1 (MQ=21)
                                                                               cgcTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGt               >  2:110867/1‑101 (MQ=37)
                                                                               cgcTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGt               >  1:3426220/1‑101 (MQ=37)
                                                                                gcTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTc              >  1:140895/1‑101 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc             >  2:2596578/1‑101 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc             >  1:2123970/1‑101 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc             <  2:1726491/101‑1 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc             >  1:1080083/1‑101 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGCAGTcc             >  2:725131/1‑101 (MQ=25)
                                                                                 cTGCCCCCCT.ATTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCAACGGGAAAAAGAACGGCGTTNNNNNTGGTCGTAATACACCGGGGGGAggtc             >  1:2393677/1‑98 (MQ=12)
                                                                                  tGCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCg            >  1:39017/1‑101 (MQ=37)
                                                                                   gCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGc           >  2:1924402/1‑101 (MQ=37)
                                                                                   gCCCCCCTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGCCAAAAAACGCCGGGGGGGGTCCGc           >  1:4033166/1‑101 (MQ=25)
                                                                                    ccccccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt          >  2:1549671/1‑101 (MQ=37)
                                                                                    ccccccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt          >  1:2051615/1‑101 (MQ=37)
                                                                                    ccccccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt          <  2:807826/101‑1 (MQ=37)
                                                                                     cccccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt         <  1:694534/101‑1 (MQ=37)
                                                                                     cccccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt         <  1:4526498/101‑1 (MQ=37)
                                                                                     cccccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt         <  2:1330710/101‑1 (MQ=37)
                                                                                     cccccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt         <  1:3799372/101‑1 (MQ=37)
                                                                                      ccccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTg        >  1:3595983/1‑101 (MQ=37)
                                                                                      ccccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGACCGCTTg        <  2:444012/101‑1 (MQ=37)
                                                                                       cccTC.TTTTGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGc       <  1:2745804/101‑1 (MQ=37)
                                                                                         ccttttttGGCC..GCCGCAAAATTCCCCCCCCGGGCTTTTTTCCACAAGAAAAAACCCGGGCGTTTTCCCTGGCCAAAATCCGCCGGGGGGAGTCCGCTTGc       <  2:1604907/98‑1 (MQ=255)
                                                                                          ctgttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCgg     <  1:2123661/98‑1 (MQ=38)
                                                                                           tcttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGa    <  1:2118584/99‑1 (MQ=37)
                                                                                           tcttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGa    >  2:3864476/3‑101 (MQ=37)
                                                                                           tcttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGa    >  1:2828849/3‑101 (MQ=37)
                                                                                           tcttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGa    >  1:3262243/3‑101 (MQ=37)
                                                                                            cttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc   >  1:2277608/2‑101 (MQ=255)
                                                                                            cttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc   >  1:1413104/2‑101 (MQ=255)
                                                                                            cttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc   >  1:1260567/2‑101 (MQ=255)
                                                                                            cttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc   <  2:2688617/100‑1 (MQ=255)
                                                                                             ttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACg  >  2:3545772/1‑101 (MQ=255)
                                                                                             ttttGGCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACg  >  2:2012855/1‑101 (MQ=255)
                                                                                                     ||                                                                                              
GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTGGCT..GC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACG  >  NC_014034/904834‑905025

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG