Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,942 G→C 2.2% H239D (CAT→GAT)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9420G→C 2.2% 6.0 ‑1.2 ‑0.3334H239D (CAT→GAT) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/2):  ref base (177/149):  total (183/151)
Fisher's exact test strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.17e-01

Only 100 of 364 total aligned reads displayed.

ACGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCC  >  NC_014034/221842‑222042
                                                                                                    |                                                                                                    
 gggTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGc                                                                                                     <  1:5183426/100‑1 (MQ=255)
 cGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGc                                                                                                     <  2:6392342/101‑1 (MQ=255)
 cGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGc                                                                                                     <  1:6189470/101‑1 (MQ=255)
  ggTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCg                                                                                                    >  1:2380313/1‑101 (MQ=255)
  ggTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATAGACGGTGGGCGTCGGATGCg                                                                                                    >  2:2408180/1‑101 (MQ=255)
       aGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGc                                                                                               >  2:6797928/1‑101 (MQ=255)
       aGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGc                                                                                               >  2:2931387/1‑101 (MQ=255)
         ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtc                                                                                             <  1:1151373/101‑1 (MQ=255)
           ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                           >  1:6062666/1‑101 (MQ=255)
             tttCGCCGGTGCCCGACACCGTCGCCNNNNNTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                         <  1:4494665/101‑1 (MQ=255)
             tttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                         <  2:1682678/101‑1 (MQ=255)
              ttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtcg                                                                                        <  2:202743/101‑1 (MQ=255)
               tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGGGGGCGTCGGACGCGGGGGCtcggccgg                                                                                       >  2:734305/1‑101 (MQ=255)
                  ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                    <  1:623620/101‑1 (MQ=255)
                          cgcccccGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatctgtcaggt                                                                            <  2:2631359/96‑1 (MQ=255)
                           gACGCCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                           <  1:2900653/101‑1 (MQ=255)
                            acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                          <  1:4369024/101‑1 (MQ=255)
                            acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggcc                                                                          >  1:5370426/1‑101 (MQ=255)
                              cccGTCGCGGAGGCTTGTGTGTTGGTGGTCAGTGTCGCCGGATCGCGCTCATCGACAGTGGGCGCCGGATGCGTGGGCTCGGTCggatcggtcgggtccgt                                                                        <  2:3561496/100‑1 (MQ=255)
                                 gTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggg                                                                     <  1:4266374/101‑1 (MQ=255)
                                 gTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggg                                                                     >  2:1105305/1‑101 (MQ=255)
                                   cGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggat                                                                   <  2:1418312/101‑1 (MQ=255)
                                    gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                  <  1:4961294/101‑1 (MQ=255)
                                     ccGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                                 <  2:2508399/101‑1 (MQ=255)
                                       gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                               >  1:628975/1‑101 (MQ=255)
                                       gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                               >  1:4726904/1‑101 (MQ=255)
                                        aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                              >  1:1049958/1‑101 (MQ=255)
                                           cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                           >  2:4251007/1‑101 (MQ=255)
                                            ttCGGTGTTTGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggttcggtcgggtccgggggttgggtcgggc                                                          >  1:6505431/1‑100 (MQ=255)
                                            ttCGGTGTTGGTGGTCAGCGTCGCCGGATCCCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggcccgggggatcggtcgggt                                                          >  1:4639262/1‑101 (MQ=255)
                                             tCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggagcggccgggtc                                                         >  1:4540889/1‑101 (MQ=255)
                                              ttccgatcttggtCAGCGTAGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccttgggatcggtcgggtcc                                                        <  1:4234736/92‑1 (MQ=38)
                                              cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                        >  1:3309786/1‑101 (MQ=255)
                                              cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                        <  1:4022613/101‑1 (MQ=255)
                                              cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                        >  2:562997/1‑101 (MQ=255)
                                              cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                        >  1:5913848/1‑101 (MQ=255)
                                              cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgggggatcggtcgggtcc                                                        >  2:3240311/1‑101 (MQ=255)
                                                gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                      <  1:2991642/101‑1 (MQ=255)
                                                 tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                     >  2:5546961/1‑101 (MQ=255)
                                                 tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                     >  1:220678/1‑101 (MQ=255)
                                                 tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                     >  1:2816279/1‑101 (MQ=255)
                                                 tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcgggcgggtccgtg                                                     >  1:2634251/1‑101 (MQ=255)
                                                   ttGGTGGTCAGCGTCGCCGGGTCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtgggatcggtcgggtccgtggg                                                   <  1:4917878/101‑1 (MQ=255)
                                                   ttGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtgggatcggtcgggtccgtggg                                                   <  1:88899/101‑1 (MQ=255)
                                                     ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGAtcggtcgggtccgtgggatcggtcgggtccgtgggac                                                 >  2:1872934/1‑100 (MQ=255)
                                                     ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATcggtcgggtccgtgggatcggtcgggtccgtgggat                                                 <  1:3894951/101‑1 (MQ=255)
                                                     ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATcggtcgggtccgtgggatcggtcgggtccgtgggat                                                 <  1:191040/101‑1 (MQ=255)
                                                      gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                <  1:4819107/101‑1 (MQ=255)
                                                       tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtggggtcggtcggggccgggggggcg                                               >  1:5437386/1‑97 (MQ=38)
                                                       tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                               >  2:3163479/1‑101 (MQ=255)
                                                       tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                               <  1:5019542/101‑1 (MQ=255)
                                                          tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                            <  2:16096/101‑1 (MQ=255)
                                                           gaGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                           <  2:3072634/100‑1 (MQ=255)
                                                           cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggcccgtgggatcggtcg                                           >  2:4308827/1‑101 (MQ=255)
                                                            aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                          <  2:4486206/101‑1 (MQ=255)
                                                            aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                          <  2:2188533/101‑1 (MQ=255)
                                                               gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                       <  1:1238061/101‑1 (MQ=255)
                                                               gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                       <  2:6250572/101‑1 (MQ=255)
                                                                tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                      <  1:2480739/101‑1 (MQ=255)
                                                                gcGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                      <  2:4809991/100‑1 (MQ=255)
                                                                 cGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGACCGGTCGGGTCCGTGGGATCGGCCGGGGCCGTGGGATCGGTGGGATCCg                                     >  2:6666204/1‑101 (MQ=37)
                                                                  gCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt                                    <  2:6062354/101‑1 (MQ=255)
                                                                    cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                  >  1:1786651/1‑101 (MQ=255)
                                                                    cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                  >  1:4187882/1‑101 (MQ=255)
                                                                    cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                  >  1:3601882/1‑101 (MQ=255)
                                                                    cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGGCCGGGGGATCGGGCGGGTCCCTGGGATCGGGCGGATCCGTgg                                  >  2:5928657/1‑101 (MQ=255)
                                                                      gATCGCGCTCATCGAGGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                >  1:6874613/1‑101 (MQ=255)
                                                                      gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGa                                >  1:7058460/1‑101 (MQ=255)
                                                                       aTCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                               <  1:3121754/101‑1 (MQ=255)
                                                                       aTCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGATCCGTggggt                               >  1:1071136/1‑99 (MQ=255)
                                                                        tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                              <  1:1635346/101‑1 (MQ=255)
                                                                        tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                              <  2:2361629/101‑1 (MQ=255)
                                                                        tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                              <  1:5545208/101‑1 (MQ=255)
                                                                           cgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                           >  1:2932394/1‑101 (MQ=255)
                                                                            cttccgatctggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                          <  1:1302581/91‑1 (MQ=255)
                                                                               gaggcgtcgTGGTCGGCGGATCGGTCGGTTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcgga                       <  1:2059995/93‑2 (MQ=11)
                                                                               cATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGACGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCggg                       >  2:3452329/1‑101 (MQ=255)
                                                                                 tCGACGGTGGGGGTCGGATGCGTGGGCTCGGTCGGGTCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                     <  2:4701319/101‑1 (MQ=255)
                                                                                    gtgggctcggTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGGCGGGTCCGTGGGATCCCGCGGGTCCgc                  >  1:3633312/10‑100 (MQ=14)
                                                                                    gcGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                  <  2:5370154/100‑1 (MQ=255)
                                                                                    aCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                  >  1:674970/1‑101 (MQ=255)
                                                                                      ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcagatcggaaa                >  2:1302533/1‑91 (MQ=255)
                                                                                      ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                >  1:6654342/1‑101 (MQ=255)
                                                                                       gTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg               <  2:6515892/101‑1 (MQ=255)
                                                                                        tGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa              >  2:3711931/1‑101 (MQ=255)
                                                                                        tGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGa              >  2:5170611/1‑101 (MQ=255)
                                                                                        ggggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa              <  1:5123482/100‑1 (MQ=255)
                                                                                         gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt             >  1:3432788/1‑101 (MQ=255)
                                                                                         ggccGTTGGCTGTGTCGGGTCCGTGGGCTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt             <  2:2811412/98‑1 (MQ=11)
                                                                                           gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg           <  2:1239692/101‑1 (MQ=255)
                                                                                           gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg           <  1:1591205/101‑1 (MQ=255)
                                                                                             gttggATCCGTCGGCTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt         <  2:5466559/98‑1 (MQ=11)
                                                                                             gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGTATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt         >  1:2839059/1‑101 (MQ=255)
                                                                                             gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggg         >  1:2721378/1‑100 (MQ=255)
                                                                                             gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGAgcggg         >  2:2737270/1‑100 (MQ=17)
                                                                                              tCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtccctg        >  2:5752656/1‑97 (MQ=11)
                                                                                              gcGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcagtc        <  1:3180154/100‑1 (MQ=255)
                                                                                                 gATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa     >  2:6849753/1‑101 (MQ=255)
                                                                                                 gATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa     >  2:3894691/1‑101 (MQ=255)
                                                                                                    gCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc  <  1:1353545/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
ACGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCC  >  NC_014034/221842‑222042

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG