Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,000,990 +C 1.9% intergenic (+34/+89) bhbA → / ← bfr methylmalonyl‑CoA mutase/bacterioferritin

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,000,9901.→C 1.9% 2.2 ‑3.3 ‑0.2159intergenic (+34/+89)bhbA/bfrmethylmalonyl‑CoA mutase/bacterioferritin
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (51/105):  total (54/105)
Fisher's exact test strand distribution p-value = 3.77e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.92e-01

Only 100 of 214 total aligned reads displayed.

GCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTATTCGGCCGC  >  NC_014034/1000890‑1001090
                                                                                                        |                                                                                                       
  cgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGC.TG.cc                                                                                                       >  1:971708/1‑101 (MQ=255)
   gggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.ccc                                                                                                      <  1:3385289/101‑1 (MQ=255)
     gcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.ccccc                                                                                                    <  2:4605458/101‑1 (MQ=255)
       cccAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCg                                                                                                  <  2:1801932/100‑1 (MQ=255)
       aCCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGC.TG.ccccccc                                                                                                  >  1:2183791/1‑100 (MQ=255)
        ccAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.cccccccc                                                                                                 >  1:3238059/1‑99 (MQ=255)
         cAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGG.TG.CCCCCCG.cc                                                                                               >  2:4868302/1‑101 (MQ=255)
         cAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.cc                                                                                               <  1:715081/101‑1 (MQ=255)
          aaCATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGAAGATGACACGGGCGAGGGGGGCGC.CG.CCCCCCG.CCt                                                                                              >  2:640856/1‑101 (MQ=255)
           aCATTCCGGGGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGGGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCtt                                                                                             <  1:3653813/101‑1 (MQ=255)
            cttccgggcggcggcGCGCGAGATCCTGGATCTGATCCGC.GCGGCGGGGGGCTG.ACAGGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCttt                                                                                            <  2:339076/95‑1 (MQ=255)
            cATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCC.CCctt                                                                                            >  2:2547440/1‑101 (MQ=255)
            cATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCACG.CCttt                                                                                            <  1:1648358/101‑1 (MQ=255)
             ctcccGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTg                                                                                           <  2:4941600/98‑1 (MQ=255)
             atcccGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTg                                                                                           <  1:4849812/98‑1 (MQ=255)
                ccGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCCTTGG..cc                                                                                      >  2:1506277/1‑101 (MQ=255)
                ccGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCACCCG.CCCTTgggc                                                                                        >  2:484662/1‑99 (MQ=255)
                  ggcggcggcGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGG.GG.CCCCCCC.CCCCTCG..cccc                                                                                    >  1:1052216/1‑101 (MQ=255)
                  ggcggcggcGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..cccc                                                                                    <  1:2476166/101‑1 (MQ=255)
                  ggcggcggcGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..cccc                                                                                    >  2:1451992/1‑101 (MQ=255)
                  ggcggcggcGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..cccc                                                                                    <  2:660499/101‑1 (MQ=255)
                  ggcggcggcGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCCTTGG..gccc                                                                                    >  2:2381594/1‑101 (MQ=255)
                  ggcagcggcGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCCTTGG..cccc                                                                                    >  1:5015431/1‑101 (MQ=255)
                   gcggcggcgCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGC.TG.CCCCCCG.CCCTTGG..ccccc                                                                                   >  2:380273/1‑101 (MQ=255)
                   gagggggcgCGCGAGATCCTGGATCTGACCCGC.GCGGCGCGCGACTG.ACACGGGGCAGCTGACCCGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..ccccc                                                                                   <  2:4996736/96‑1 (MQ=255)
                     ggcggcgCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCg                                                                                 >  1:3129391/1‑101 (MQ=255)
                     ggcggcgCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCg                                                                                 >  2:580240/1‑101 (MQ=255)
                     ggcggcgCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCg                                                                                 <  1:4008688/101‑1 (MQ=255)
                          ggcgcgAGATCCTGGATCTGATCCGC.GGGGCGCGGGGCTG.ACACGGGGCAGCTCACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGata                                                                            <  2:4370000/100‑1 (MQ=255)
                           gggcgAGGTCGTGGATCTGAGCCGC.GGGGCG.GGGGCTGAACACGGGGCAGCTGAAACGGGGCTGGGGGGGGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGatat                                                                           <  2:2233339/99‑1 (MQ=255)
                            cgcgAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATAtt                                                                          <  1:4695250/101‑1 (MQ=255)
                             gggagaTCCTGGATCTGACCCGC.GCGGGGGGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATAttt                                                                         <  1:1329789/99‑1 (MQ=255)
                             gcgAGATCGTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATAttt                                                                         <  2:129450/101‑1 (MQ=255)
                             gcgAGAGCTTGGAT.TGGTCCCCGGCGGCGTGGGACTG.ACACGGGGAAGCGGCCACGGGGCGGGGGGGCAC.TT.CCCCCCG.CCTTTGG..CCCCCCGGGATAttt                                                                         <  2:2822863/101‑1 (MQ=255)
                              cgAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGTAGGGGGGCAC.CG.TCCCCCC.CCCTTTG..TGCCCCGGGACAtttt                                                                        >  2:1983874/1‑101 (MQ=255)
                              cgAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGCGGG.CCCCCCC.CCCTTGT..CCCCCCgcggccttt                                                                         >  2:15469/1‑93 (MQ=255)
                              cgAGAGCCTGGATCTGACCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATAtttt                                                                        <  1:4675976/101‑1 (MQ=255)
                                agaTCCTGGATCTGACCCGC.GGGGCGGGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTgg                                                                      <  2:1718161/101‑1 (MQ=255)
                                  aTCCTGGATCGGATCCGC.GCGGCGGGGGGCTG.ACACGGGGCAGCTGACACGGGGCGGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGc                                                                    <  1:1887239/101‑1 (MQ=255)
                                     cTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAg                                                                 >  2:2474230/1‑101 (MQ=255)
                                     cTGGATCTGACCCGC.GCGGCGGGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAg                                                                 <  1:2722195/101‑1 (MQ=255)
                                      tGGATCTGATCCGC.GCGGCGGGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGa                                                                <  1:669956/101‑1 (MQ=255)
                                      tGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGa                                                                <  2:1688848/101‑1 (MQ=255)
                                       ggATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGC.GG.CCCCCCGCCCCTCGGC.CCCCCCGGGATTTTTTgtgaagt                                                                 >  2:3643206/1‑95 (MQ=255)
                                        gATCTGATCCGC.GGGGCGGGGGGCTG.CCCGGGGGCGGCGGACAGGGGGCAGGGGGGGGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGAtt                                                              <  1:3659196/101‑1 (MQ=255)
                                        gATCGGACCCGC.GCGGGGCGGGGCTG.ACACGGGGCGGCTGACACGGGGCAGGGGGGCGC.GG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGAtt                                                              <  1:1487307/101‑1 (MQ=255)
                                         aTCTGACCCGC.GGGGGGGGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTa                                                             <  1:204995/101‑1 (MQ=38)
                                          tCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAc                                                            >  1:2154747/1‑101 (MQ=255)
                                           cTGATCCGC.GCGGCGCGGTGCTG.ACTCGGGCCAGCTGACGCGGGGCGGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACt                                                           <  2:305118/101‑1 (MQ=38)
                                            tGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGC.GG.CCCCCCC.CCCTCTG..CCCCCCCGGGTGTTTTTGGGGAGATTAATg                                                          >  1:2612770/1‑101 (MQ=255)
                                            tGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACAAGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGCTATTTTGGGCAAGATTACTg                                                          >  2:5032853/1‑101 (MQ=255)
                                             gATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGCGGG.CCCCCCC.CCTTTGGC.CCCCCCGGGCTTTTTTGGGAAAGTtttta                                                           >  2:3459989/1‑98 (MQ=255)
                                             gATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGa                                                         >  1:311884/1‑101 (MQ=255)
                                             gATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGa                                                         >  1:658642/1‑101 (MQ=255)
                                             gATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCCTTGG..CCCCCCCGGGTATTTTTGGCAAGATTACTGa                                                         >  2:1949846/1‑101 (MQ=39)
                                              gtCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGaa                                                        <  2:135467/100‑1 (MQ=255)
                                              gcccGC.GGGGGGGGGGGCTG.ACGCGGGGCAGCTGACGGGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGaa                                                        <  1:2750880/99‑1 (MQ=38)
                                              aTCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTTG..CCCCCCCGGATATTTTGGGCAAGATTACTGaa                                                        >  2:1579813/1‑101 (MQ=255)
                                              aTCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGaa                                                        >  2:3931638/1‑101 (MQ=255)
                                                ccGC.GGGGGGCGGGGCTG.ACGCGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGaaaa                                                      <  1:1506407/101‑1 (MQ=255)
                                                ccGC.GGGGGGCGGGGCTG.ACCCGGGGCGGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGaaaa                                                      <  1:1307723/101‑1 (MQ=39)
                                                ccGC.GCGGCGCGGGGCGG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGaaaa                                                      <  2:2203242/101‑1 (MQ=255)
                                                 cgc.gcgGCGCGGGGCTG.ACACGGGGCAGCGGACAGGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGaaaaa                                                     <  2:3057624/101‑1 (MQ=255)
                                                  gc.gcggcgcgGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAg                                                    >  1:2336556/1‑101 (MQ=255)
                                                      cggcgcggGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGgcg                                                 >  1:1584341/1‑101 (MQ=255)
                                                       ggggcggGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGgcgc                                                <  2:3729988/98‑1 (MQ=255)
                                                          gcggGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTcc                                             <  2:2696104/101‑1 (MQ=255)
                                                           cggGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGC.GG.GCCCCCC.CCGCTGG..GCCCCCGGGGTATTTTGGGGGAGGTTTCTGaaaaaagggggccg                                            >  1:2082455/1‑92 (MQ=255)
                                                            ggggCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGGGG.GGCCCCCCCC.CCTTTGG..CCCCCCGGG.GATTTTTGGGGAGAATAATGAAAAAAGAGCGCcgc                                           >  1:1926887/1‑101 (MQ=255)
                                                             gggCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCcgcg                                          <  2:3502074/101‑1 (MQ=255)
                                                               gctg.acacggggcagctgacacggggcaGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAg                                        >  2:711804/1‑101 (MQ=255)
                                                                 tg.acacggggcagctgacacggggcagGGGGGGGC.TGCCCCCCCC.CCTTTGG..CCCCCCGGG.TATTTTTGGGAAGATTACTGAAAAAGGCGCTCCGCGCggg                                      >  2:1128740/1‑101 (MQ=255)
                                                                    acacggggcagctgacacggggcagGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGgc                                    >  2:4689507/1‑101 (MQ=255)
                                                                        ggggcagctgacagggggcagGGGGGCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCt                                <  2:2653090/101‑1 (MQ=255)
                                                                        ggggcagctgacacggggcagGGGGGGGG.GG.CCCCCCC.CCCCTCGCCCCCCCCGGGAT.TTTTGGGC.AGAATACATAAAAAAGAGAGCCGCGCgggggggcc                                >  1:2583962/1‑96 (MQ=255)
                                                                         gggcagctgacacggggcagGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCtt                               <  1:1565626/101‑1 (MQ=255)
                                                                         gggcagctgacacggggcagGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCtt                               >  1:3149139/1‑101 (MQ=255)
                                                                          ggcagctgacaggggggggGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCttt                              <  2:1052106/101‑1 (MQ=255)
                                                                          ggcagctgacagggggcggGGGGGCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCttt                              <  1:3909843/101‑1 (MQ=255)
                                                                          gccagtcgccacggggcggGGGGGCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCttt                              <  2:2076169/92‑1 (MQ=255)
                                                                           gcggcggacacggggcagGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTg                             <  1:382910/95‑1 (MQ=255)
                                                                                tgacacggggcagGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCtt                        <  1:3326790/101‑1 (MQ=255)
                                                                                 gacacggggcagGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCttt                       >  1:4788306/1‑101 (MQ=255)
                                                                                  acagggggcagGGGGGCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtg                      <  2:3237808/101‑1 (MQ=255)
                                                                                  acacggggcagGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtg                      >  1:2836355/1‑101 (MQ=255)
                                                                                    ccggggcagGGGGGCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgtg                    <  1:4898802/100‑1 (MQ=255)
                                                                                    acggggcagGGGGGCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgtg                    <  1:4447749/101‑1 (MQ=255)
                                                                                     cggggcagGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGGGGGGGGCCTTTGACCTTtgtgg                   >  2:3784750/1‑100 (MQ=255)
                                                                                      ggggcagGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCg                  <  2:2901077/101‑1 (MQ=255)
                                                                                           agGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCt             >  1:4596027/1‑101 (MQ=255)
                                                                                            ggggggCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCtt            <  2:489041/101‑1 (MQ=255)
                                                                                            ggggggCCC.CC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCtt            <  1:3513783/101‑1 (MQ=255)
                                                                                             gggggCGC.CG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTa           <  2:4191052/101‑1 (MQ=255)
                                                                                              ggggCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTAt          <  2:2014472/101‑1 (MQ=255)
                                                                                               gggCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTAtt         <  2:1076790/101‑1 (MQ=255)
                                                                                                ggCGC.TC.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTATTc        <  1:3643478/101‑1 (MQ=255)
                                                                                                  cgc.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTATTCgg      <  2:2506946/101‑1 (MQ=255)
                                                                                                   gc.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTATTCGGc     >  2:4153911/1‑101 (MQ=255)
                                                                                                   gc.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTATTCGGc     >  1:2963823/1‑101 (MQ=255)
                                                                                                      tG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTATTCGGCcg   <  2:517190/101‑1 (MQ=255)
                                                                                                        |                                                                                                       
GCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGC.GCGGCGCGGGGCTG.ACACGGGGCAGCTGACACGGGGCAGGGGGGCGC.TG.CCCCCCG.CCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCTTATTCGGCCGC  >  NC_014034/1000890‑1001090

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 19 ≤  ATCG  < 31 ≤  ATCG  < 39 ≤  ATCG