Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,949 C→A 10.3% E236D (GAG→GAT RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9490C→A10.3% 39.9 ‑3.2 ‑1.3205E236D (GAG→GATRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (15/6):  ref base (86/97):  total (101/104)
Fisher's exact test strand distribution p-value = 3.94e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.02e-02

Only 100 of 246 total aligned reads displayed.

AGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221849‑222049
                                                                                                    |                                                                                                    
aGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGc                                                                                                      <  2:4380411/101‑1 (MQ=255)
 gCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCt                                                                                                     >  1:4041276/1‑101 (MQ=255)
  ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtc                                                                                                    <  1:2800200/101‑1 (MQ=255)
  ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtc                                                                                                    <  2:3022950/101‑1 (MQ=255)
  ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGGGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGCCGGATGCGGGGGCtc                                                                                                    >  1:1634603/1‑101 (MQ=255)
   cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcg                                                                                                   <  1:508946/101‑1 (MQ=255)
   cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcg                                                                                                   <  1:2161191/101‑1 (MQ=255)
   cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGGGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGGGGGCGCCGGATGCGTGGGCtcg                                                                                                   >  2:594138/1‑101 (MQ=255)
     gTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggt                                                                                                 <  2:1021567/101‑1 (MQ=255)
        tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCAACGACGGTGGGCGTCGGATGCGTGGGCtcggtcgg                                                                                              <  1:4098587/101‑1 (MQ=255)
         cGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCgga                                                                                             <  1:4203163/101‑1 (MQ=255)
         cGCCGGTGCCCGACACCGTCGCCGAGGCTGCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCgga                                                                                             <  1:3169516/101‑1 (MQ=255)
          gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                            <  2:3095386/101‑1 (MQ=255)
          gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGCGTTGGGGGTCAGCGTCGCCGGATCCCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCtcggtcggtc                                                                                            >  1:3734544/1‑99 (MQ=255)
           ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                           <  2:2418473/101‑1 (MQ=255)
            cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                          >  1:1831024/1‑101 (MQ=255)
            cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                          <  1:3387036/101‑1 (MQ=255)
               ggCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                       <  2:4330155/100‑1 (MQ=255)
                 cccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgg                                                                                     >  1:2323019/1‑101 (MQ=255)
                 cccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgg                                                                                     <  2:1138242/101‑1 (MQ=255)
                    gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggcc                                                                                  >  1:2047417/1‑99 (MQ=255)
                     acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                 <  1:4118551/101‑1 (MQ=255)
                        ccGTCTCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtg                                                                              <  2:144848/101‑1 (MQ=255)
                        ccGTCGCGGAGGTGTCGGTGGTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtg                                                                              <  2:881120/101‑1 (MQ=255)
                         cGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgg                                                                             >  2:4114083/1‑101 (MQ=255)
                             gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggctc                                                                         >  2:2720975/1‑101 (MQ=255)
                                  ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcg                                                                    <  2:1576670/101‑1 (MQ=255)
                                   gCTTCGGAGTTGGGGGGCAGCGTCGCCGGATCGCGCTGATCGCCGGGGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatcggtcgg                                                                   >  2:1624758/1‑101 (MQ=255)
                                    cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                                  >  2:4039131/1‑101 (MQ=255)
                                     ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcggtgggatcggtcgggt                                                                 >  1:756821/1‑101 (MQ=255)
                                     ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                 >  2:3426362/1‑101 (MQ=255)
                                      tCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtc                                                                >  1:221183/1‑101 (MQ=255)
                                       cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               <  2:1083251/101‑1 (MQ=255)
                                       cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               <  2:1405995/101‑1 (MQ=255)
                                       cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGGGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               >  1:1344715/1‑101 (MQ=255)
                                        ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatcggtcgggcccg                                                              >  1:1323444/1‑101 (MQ=255)
                                         gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                             >  2:3614997/1‑101 (MQ=255)
                                         gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                             >  2:3328580/1‑101 (MQ=255)
                                          tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                            >  2:3148172/1‑101 (MQ=255)
                                             tggtggTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGAtcggtcgggtccgtgggatcggtcgggtccgtggga                                                         <  1:534991/101‑1 (MQ=255)
                                              ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATcggtcgggtccgtgggatcggtcgggtccgtgggat                                                        <  1:1735420/101‑1 (MQ=255)
                                                 ggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGtcgggtccgtgggatcggtcgggtccgtgggatcgg                                                     <  2:2331178/101‑1 (MQ=255)
                                                     aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGccgggtccgtgggctcggtcgggtccgcgggatcggcggg                                                 >  2:2652984/1‑97 (MQ=255)
                                                     aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtggg                                                 >  2:3927896/1‑101 (MQ=255)
                                                     aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                                 >  2:3802314/1‑101 (MQ=255)
                                                     aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcgggcggggccgtgggatcggtcgg                                                 >  2:800779/1‑101 (MQ=255)
                                                      gCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa                                                <  2:416416/101‑1 (MQ=255)
                                                      gCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa                                                >  1:1946164/1‑101 (MQ=255)
                                                       cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGAt                                               >  2:3217367/1‑101 (MQ=255)
                                                         tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTTGGATCGGTCGGATcc                                             <  1:3754101/101‑1 (MQ=255)
                                                          cGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg                                            >  1:4322777/1‑101 (MQ=255)
                                                           gCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt                                           <  1:342407/101‑1 (MQ=255)
                                                            ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                          >  2:3874188/1‑101 (MQ=255)
                                                            ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                          <  1:662528/101‑1 (MQ=255)
                                                             cGGATCACGCTCATCGACGGGGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                         <  2:2794219/101‑1 (MQ=255)
                                                                    ggcTCATGGGCGGTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                  <  1:1279094/100‑1 (MQ=255)
                                                                     gcTCATCGGCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                 <  1:3515409/101‑1 (MQ=255)
                                                                     gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                 >  2:1211352/1‑101 (MQ=255)
                                                                     gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                 >  2:1788020/1‑101 (MQ=255)
                                                                     gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                 >  1:387690/1‑101 (MQ=255)
                                                                      cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                >  2:1555913/1‑101 (MQ=255)
                                                                      cTCATCGACGGGGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                <  2:3528853/101‑1 (MQ=255)
                                                                       tCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                               >  2:2535860/1‑101 (MQ=255)
                                                                         accGACGGTTGCCCTCGGATGCGTGGCCTCGTCCGGTCCGGTCGGGCCCTTGGAATCGGCCGGG.CCGGTGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                            <  1:2653156/99‑1 (MQ=9)
                                                                          tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                            <  1:2576961/101‑1 (MQ=255)
                                                                           cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                           <  2:2270314/101‑1 (MQ=255)
                                                                             aCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                         <  2:2936354/101‑1 (MQ=255)
                                                                              tgggctctgTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGtc                        <  2:3171346/93‑2 (MQ=11)
                                                                              cGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGGtcgggcggcgcccgg                        >  2:784879/1‑94 (MQ=21)
                                                                               ggTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                       <  1:3243246/101‑1 (MQ=255)
                                                                                gTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                      >  1:3352891/1‑101 (MQ=255)
                                                                                  gggCGTCGGAGGGGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                    <  2:1830900/101‑1 (MQ=255)
                                                                                   tcggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtc                   <  1:4164851/98‑1 (MQ=1)
                                                                                   ggCGTCGGATGCGTGGGCTCGGGCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                   <  2:3399409/101‑1 (MQ=255)
                                                                                   ggCGCCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                   <  2:1323347/101‑1 (MQ=255)
                                                                                    gCGTGGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                  <  2:3272893/101‑1 (MQ=255)
                                                                                    gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                  >  1:237556/1‑101 (MQ=255)
                                                                                    gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                  <  1:2234615/101‑1 (MQ=255)
                                                                                    cggtggggtccGTGGGATCGGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                  <  1:3119913/91‑1 (MQ=11)
                                                                                      gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                <  2:2047280/101‑1 (MQ=255)
                                                                                       tgggttccGTTGGGTCCGTCGGATCGGTTGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               <  1:450795/94‑1 (MQ=11)
                                                                                       tgggATGGGTCGGGTCGCTCGGATCGGTGAGATCCGTGGGAGCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTCGGGTCGGTCGGGTCCGTGGGAtcggtc               <  2:888880/99‑1 (MQ=2)
                                                                                       tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               <  1:1116024/101‑1 (MQ=255)
                                                                                       tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               <  2:2328403/101‑1 (MQ=255)
                                                                                        tgtttccGTGGGATCCGTCGGTTCCGTGGGATCCGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGAtcggtcg              <  2:2393703/95‑1 (MQ=11)
                                                                                        ggggtccGTGGGGTCGGTGGGATCGGT.GGGTCCCTGGGGTTCGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccgtcgg             <  1:950923/95‑1 (MQ=2)
                                                                                         ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             <  1:2771770/101‑1 (MQ=255)
                                                                                         ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGCGCGAACGCACGGATCCGAGGGATCGGTCGGGTCCGTGGGAtcggtcgg             >  2:3681979/1‑101 (MQ=18)
                                                                                           gggggTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt           <  1:1029514/97‑1 (MQ=255)
                                                                                             gcggtggCTCGGTGGGGGCGGTGGGGGGGGTGGGGCCGGGAGGGGCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc         <  1:3050787/96‑1 (MQ=9)
                                                                                              gctcggCTCCGTAGGATCGGTCGGGTCCGTGGGTTCGCTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg        <  2:2822447/97‑1 (MQ=11)
                                                                                               gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt       <  2:4057471/101‑1 (MQ=255)
                                                                                               gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCGGTCGGATCGGTCGGGTCCGTGGGGtcgggggggaccgt       >  2:450755/1‑91 (MQ=11)
                                                                                               gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt       >  1:2108375/1‑101 (MQ=1)
                                                                                                tGGGCGCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      <  2:1807368/101‑1 (MQ=255)
                                                                                                 gggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg     <  1:4380854/101‑1 (MQ=255)
                                                                                                 gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg     <  2:1974226/101‑3 (MQ=11)
                                                                                                 gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg     >  1:3651038/1‑99 (MQ=11)
                                                                                                 gggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGAccgctcggggcggcgg     >  2:4067658/1‑93 (MQ=11)
                                                                                                  ggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg    <  2:3822963/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
AGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221849‑222049

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