Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,155,466 G→A 38.3% V11M (GTC→ATC)  RCAP_rcc02000 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,155,4660G→A38.3% 285.5 0.0 ‑0.2227V11M (GTC→ATC) RCAP_rcc02000hypothetical protein
Reads supporting (aligned to +/- strand):  new base (49/38):  ref base (80/60):  total (129/98)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.68e-01

Only 100 of 240 total aligned reads displayed.

ACCTTCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTT  >  NC_014034/2155366‑2155566
                                                                                                    |                                                                                                    
cccTTCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGACAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCtcg                                                                                                      <  2:3342103/100‑1 (MQ=255)
  cTTCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCtcgtc                                                                                                    >  1:526955/1‑101 (MQ=255)
    tCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGa                                                                                                  >  2:1283786/1‑101 (MQ=255)
        aTGCCCTGGGGCCCCGCCCCCGTCCCCGCCCGCGGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGcc                                                                                              <  2:3089802/101‑1 (MQ=255)
          gCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCAt                                                                                            >  2:603609/1‑101 (MQ=255)
          ccccTCGGGCGCCGTCCCCGCCCCCGCCCGCGGGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCAt                                                                                            <  2:370273/100‑1 (MQ=255)
            ccTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCt                                                                                          >  2:3211237/1‑101 (MQ=255)
                 ggCGCCGCCCCCGTCCCCGCCCGCCGGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTggggc                                                                                     <  2:39119/101‑1 (MQ=255)
                        tcaccgtcaccgCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCg                                                                              >  2:2820582/1‑101 (MQ=255)
                         caccgtcaccgCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCgg                                                                             >  2:1729015/1‑100 (MQ=255)
                         caccgtcaccgCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGt                                                                             >  1:192955/1‑101 (MQ=255)
                          ccctcccccccccgcgGGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTc                                                                            <  1:2956785/91‑1 (MQ=255)
                             gtcaccgCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcacc                                                                         >  2:2997949/1‑101 (MQ=255)
                             gcccccgcccgcCGGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcacc                                                                         <  2:514284/97‑1 (MQ=255)
                             ccccccccccgcgGGGCGAAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcacc                                                                         <  2:827554/94‑1 (MQ=255)
                              cccccccccgccGGGCGAAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcacca                                                                        <  2:2707371/95‑1 (MQ=255)
                                 ccgcccgcCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGa                                                                     >  1:672984/1‑101 (MQ=255)
                                    cccgccGGGCGCAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATgg                                                                  <  2:4153661/101‑1 (MQ=255)
                                     ccgccgGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggc                                                                 <  2:4001591/101‑1 (MQ=255)
                                      cgccgGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggcg                                                                <  1:2505706/101‑1 (MQ=255)
                                          gggCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggcggcgg                                                            >  2:3789401/1‑101 (MQ=255)
                                              ggcAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCt                                                        >  1:1625843/1‑101 (MQ=255)
                                                 atctcgcgAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGc                                                     <  2:3634651/97‑1 (MQ=255)
                                                  aGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGcc                                                    <  1:3220321/101‑1 (MQ=255)
                                                  aGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGcc                                                    >  1:416624/1‑101 (MQ=255)
                                                  aGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGcc                                                    <  1:2423696/101‑1 (MQ=255)
                                                   gCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCg                                                   <  2:678858/101‑1 (MQ=255)
                                                     cgcgAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGAGGCCGTCACCACGATGGCGGGGGCGCTCGcagat                                                 >  1:3634904/1‑97 (MQ=255)
                                                        gAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcgg                                              >  1:813513/1‑101 (MQ=255)
                                                        gAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcgg                                              >  1:348982/1‑101 (MQ=255)
                                                        gAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcgg                                              >  2:3490116/1‑101 (MQ=255)
                                                         aaaGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcggc                                             <  2:2606885/101‑1 (MQ=255)
                                                         aaaGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcggc                                             <  1:1167286/101‑1 (MQ=255)
                                                          aaGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCt                                            >  2:255528/1‑101 (MQ=255)
                                                           aGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTg                                           >  1:2083659/1‑101 (MQ=255)
                                                             ggCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGAt                                         <  1:423630/101‑1 (MQ=255)
                                                              gCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGAtg                                        >  2:3569839/1‑101 (MQ=255)
                                                                aaaaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGAtgtg                                      <  2:1306654/95‑1 (MQ=255)
                                                                aGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGAtgtg                                      >  2:2180712/1‑101 (MQ=255)
                                                                aGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGAtgtg                                      <  1:2648126/101‑1 (MQ=255)
                                                                 ggATGATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTgg                                     >  2:3582362/1‑101 (MQ=255)
                                                                 ggATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTgg                                     <  2:2392573/101‑1 (MQ=255)
                                                                 ggATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTgg                                     <  2:2665430/101‑1 (MQ=255)
                                                                 aaaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTgg                                     >  1:3722990/6‑101 (MQ=255)
                                                                  gATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGc                                    <  1:3140995/101‑1 (MQ=255)
                                                                  aaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGc                                    <  1:2137272/97‑1 (MQ=255)
                                                                   aTCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCa                                   >  1:967259/1‑101 (MQ=255)
                                                                    aaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCAt                                  >  2:3817785/3‑101 (MQ=255)
                                                                     cATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGGTGTGGCAtt                                 >  1:145185/1‑101 (MQ=255)
                                                                     aaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCAtt                                 >  1:2891517/2‑101 (MQ=255)
                                                                     aaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCAtt                                 <  1:1829486/100‑1 (MQ=255)
                                                                      aTGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTc                                <  1:3321350/101‑1 (MQ=255)
                                                                      aTGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTc                                <  2:1711026/101‑1 (MQ=255)
                                                                      aTGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTc                                <  1:4025019/101‑1 (MQ=255)
                                                                      aTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTc                                <  1:1420867/101‑1 (MQ=255)
                                                                       tGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCg                               <  2:201665/101‑1 (MQ=255)
                                                                       tGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCg                               <  2:2307370/101‑1 (MQ=255)
                                                                       tGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCg                               <  2:1977537/101‑1 (MQ=255)
                                                                        gaccGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCCTTCgc                              >  2:296935/4‑100 (MQ=255)
                                                                         aGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGAc                             >  2:6146/1‑101 (MQ=255)
                                                                          gCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGAcc                            <  1:836526/101‑1 (MQ=255)
                                                                          ccGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCNTGATGTGGCATTCGAcc                            >  2:746756/2‑101 (MQ=255)
                                                                           cGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCg                           <  2:409275/101‑1 (MQ=255)
                                                                            gataccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGa                          >  2:2633700/10‑101 (MQ=255)
                                                                            gATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGa                          >  2:3529522/1‑101 (MQ=255)
                                                                             ataccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAg                         <  2:1428082/93‑1 (MQ=255)
                                                                             ataccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAg                         <  1:545536/93‑1 (MQ=255)
                                                                             ataccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAg                         >  1:4056956/9‑101 (MQ=255)
                                                                               accagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAggg                       >  2:874151/7‑100 (MQ=255)
                                                                               accagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGt                       >  2:2935578/7‑101 (MQ=255)
                                                                               accagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGt                       <  1:3754808/95‑1 (MQ=255)
                                                                                cAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTg                      >  1:1855685/1‑101 (MQ=255)
                                                                                cAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTg                      <  2:3437129/101‑1 (MQ=255)
                                                                                 cagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCCCGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGc                     >  2:1852480/5‑101 (MQ=255)
                                                                                  agaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCg                    <  1:1997935/98‑1 (MQ=255)
                                                                                  aGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCg                    <  1:2884948/101‑1 (MQ=255)
                                                                                    aggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCggg                  >  2:2254729/2‑101 (MQ=255)
                                                                                    aggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCggg                  <  2:1693469/100‑1 (MQ=255)
                                                                                    aggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCggg                  >  2:109298/2‑101 (MQ=255)
                                                                                       ccTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgc               >  1:10244/1‑101 (MQ=255)
                                                                                        cTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcg              >  2:1542037/1‑101 (MQ=255)
                                                                                         tCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGgcgcgc             >  2:46546/1‑101 (MQ=255)
                                                                                          cGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggc            <  2:1124918/101‑3 (MQ=255)
                                                                                          cGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggc            >  1:1729013/1‑99 (MQ=255)
                                                                                          cGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCa            <  2:423015/101‑1 (MQ=255)
                                                                                           gtcgagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcg           <  2:3255501/98‑4 (MQ=255)
                                                                                           gtcgagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcg           <  1:1200421/98‑4 (MQ=255)
                                                                                           gtcgagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcg           <  1:2684276/98‑4 (MQ=255)
                                                                                           gTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCaa           >  1:4199577/1‑101 (MQ=255)
                                                                                            tcgagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcgg          >  1:753460/3‑97 (MQ=255)
                                                                                            tAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAg          >  1:3949223/1‑101 (MQ=255)
                                                                                              gagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcggcg        <  2:941081/101‑7 (MQ=255)
                                                                                              gagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcggcg        >  1:1690267/1‑95 (MQ=255)
                                                                                              gagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcggcg        >  1:2874483/1‑95 (MQ=255)
                                                                                              gagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGCGCGGgcgcggcggcg        >  1:2220779/1‑95 (MQ=255)
                                                                                              gagCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAgcg        <  1:875645/101‑1 (MQ=255)
                                                                                               agCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAgcgc       >  1:3590648/1‑101 (MQ=255)
                                                                                                 cTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGcgccg     <  1:1555068/101‑1 (MQ=255)
                                                                                                  tcgtcgATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGcgccgc    >  1:3746507/1‑101 (MQ=255)
                                                                                                    gtcgATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCtt  >  1:3967934/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
ACCTTCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTT  >  NC_014034/2155366‑2155566

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG