Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,472 A→T 5.1% G118G (GGA→GCT RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,4720A→T 5.1% 26.0 0.0 ‑1.2239G118G (GGA→GGTRCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/6):  ref base (112/114):  total (118/121)
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.86e-02

Only 100 of 252 total aligned reads displayed.

CGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAGCCGCCCGCGCCA  >  NC_014034/1396372‑1396571
                                                                                                    |                                                                                                    
cGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGa                                                                                                      >  1:1425523/1‑101 (MQ=255)
cGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGGGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGGGGATCGAAAGGGa                                                                                                      >  1:3257563/1‑101 (MQ=255)
 gCTCCCCCCGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGAc                                                                                                     <  1:3397282/101‑1 (MQ=255)
 gCTCCCCCCGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGAc                                                                                                     <  2:3921422/101‑1 (MQ=255)
   tCCCCCCGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCt                                                                                                   >  1:4206776/1‑101 (MQ=255)
   tCCCCCCGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCt                                                                                                   <  1:476245/101‑1 (MQ=255)
    ccccccGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTg                                                                                                  <  1:1982812/101‑1 (MQ=255)
    ccccccGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTg                                                                                                  <  1:2037939/101‑1 (MQ=255)
     cccccGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGa                                                                                                 <  1:262265/101‑1 (MQ=255)
      ccccGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTgga                                                                                                >  2:325454/1‑99 (MQ=255)
       cccGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGaaa                                                                                               <  2:3740313/101‑1 (MQ=255)
       cccGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGaaa                                                                                               <  1:1819040/101‑1 (MQ=255)
        ccGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAAt                                                                                              >  1:908240/1‑101 (MQ=255)
         cgcgGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAAtt                                                                                             <  2:2597176/101‑1 (MQ=255)
          gcgGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTc                                                                                            <  2:4018675/101‑1 (MQ=255)
           cgGGGGCGTGGCTCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTcc                                                                                           <  2:2030573/101‑1 (MQ=255)
            gggggCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCcg                                                                                          <  1:3322555/101‑1 (MQ=255)
            gggggCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCcg                                                                                          >  1:4630325/1‑101 (MQ=255)
               ggCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCgg                                                                                       <  1:4520174/101‑1 (MQ=255)
                gCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCggg                                                                                      <  2:313958/101‑1 (MQ=255)
                  gTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGaa                                                                                    <  1:2370376/101‑1 (MQ=255)
                   tGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGaaa                                                                                   <  1:3132151/101‑1 (MQ=255)
                   tGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGaaa                                                                                   >  1:3627364/1‑101 (MQ=255)
                   tGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGaaa                                                                                   <  2:2013305/101‑1 (MQ=255)
                   tGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGaaa                                                                                   <  2:459561/101‑1 (MQ=255)
                    ggATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAAc                                                                                  <  1:1786399/101‑1 (MQ=255)
                    ggATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAAc                                                                                  >  1:3655181/1‑101 (MQ=255)
                      aTCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAAcgc                                                                                <  1:4142964/101‑1 (MQ=255)
                       tCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAAcgcg                                                                               <  1:88343/101‑1 (MQ=255)
                        cGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACgcgg                                                                              >  1:3831140/1‑101 (MQ=255)
                        cGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACgcgg                                                                              <  2:887170/101‑1 (MQ=255)
                        cGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACgcgg                                                                              >  1:122207/1‑101 (MQ=255)
                          aaaCCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACgcggcg                                                                            >  2:3383001/1‑101 (MQ=255)
                          aaaCCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACgcggcg                                                                            >  2:2148136/1‑101 (MQ=255)
                              cgcgCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGAGGGAAACGCGGCGGCGa                                                                        >  2:656668/1‑101 (MQ=255)
                               gcgcCCTGCGGGTGATGGCGACGGGCACCANGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAg                                                                       <  1:2287833/101‑1 (MQ=255)
                                cgcGCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGc                                                                      <  2:3045405/101‑1 (MQ=255)
                                    cTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.Gt                                                                 >  2:695276/1‑101 (MQ=255)
                                        gggTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCt                                                             >  1:1670954/1‑101 (MQ=255)
                                          gTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGATCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTcc                                                           <  2:1175158/101‑1 (MQ=255)
                                          gTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTcc                                                           >  1:61539/1‑101 (MQ=255)
                                          gTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGAAC.GTCGCTcc                                                           >  2:2856351/1‑101 (MQ=255)
                                           tGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTccc                                                          <  2:1017312/101‑1 (MQ=255)
                                           tGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGCGAAACGCGGCGGCGAGCAC.GTCGCTccc                                                          >  2:4423832/1‑101 (MQ=255)
                                              tGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTcccccc                                                       >  2:2212410/1‑101 (MQ=255)
                                              tGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTcccccc                                                       <  1:3574424/101‑1 (MQ=255)
                                                gCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCcgc                                                     >  2:3276942/1‑101 (MQ=255)
                                                 cGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCcgcg                                                    <  2:1086654/101‑1 (MQ=255)
                                                   gcGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCggg                                                  <  2:3484370/100‑1 (MQ=255)
                                                    cGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCgggg                                                 <  2:2786883/101‑1 (MQ=255)
                                                     gggCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCggggg                                                >  2:792162/1‑101 (MQ=255)
                                                     gggCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCggggg                                                >  2:361150/1‑101 (MQ=255)
                                                     gggCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCggggg                                                <  2:1400194/101‑1 (MQ=255)
                                                     gggCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCggggg                                                >  2:1933956/1‑101 (MQ=255)
                                                      tctaccacGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGc                                               <  1:672878/98‑1 (MQ=255)
                                                        caccacGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGt                                             >  1:3023148/1‑101 (MQ=255)
                                                         accgcGGCGGTCGCTCCCCCTGCGGGGGCGTGGATCGAAACCGGGCTGAAAATGCGCGGGCAACGCGGCGTCGGGCAG.GTCGCCCCCCCGGCGTGGGCGTg                                            >  2:293710/1‑101 (MQ=16)
                                                         accacGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTg                                            <  2:3830855/101‑1 (MQ=255)
                                                         accacGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTg                                            <  2:2348878/101‑1 (MQ=255)
                                                         accacGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTg                                            <  1:392249/101‑1 (MQ=255)
                                                          ccacGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTgg                                           >  1:3769262/1‑101 (MQ=255)
                                                          ccacGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTgg                                           <  1:3276261/101‑1 (MQ=255)
                                                            tgcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACA..CCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATc                                        <  2:1819199/95‑1 (MQ=2)
                                                            acGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                         <  2:2995310/101‑1 (MQ=255)
                                                            acGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                         >  1:3928672/1‑101 (MQ=255)
                                                               gcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                      >  2:4143666/1‑101 (MQ=255)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACA..CCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgct                                <  1:852934/101‑4 (MQ=2)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGc                                <  2:1355256/101‑1 (MQ=255)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACA..CCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgctg                               <  2:2527058/101‑5 (MQ=0)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACA..CCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgctg                               <  1:350411/101‑5 (MQ=0)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACA..CCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgctg                               <  1:4303982/101‑5 (MQ=0)
                                                                        ccccccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCttt                             <  1:1650489/101‑1 (MQ=255)
                                                                         cccccGGGGGGGGGTGGATCGAAACGGACCGGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTc                            <  2:470791/101‑1 (MQ=255)
                                                                         cccccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTc                            >  2:855656/1‑101 (MQ=255)
                                                                         cccccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTc                            >  2:4371832/1‑101 (MQ=255)
                                                                          ccccGCGGGGGCGTGGATCGAAACGGGCCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCg                           <  1:2938967/101‑1 (MQ=255)
                                                                          ccccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCg                           >  2:3441548/1‑101 (MQ=255)
                                                                            ccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTg                         >  1:4380170/1‑101 (MQ=255)
                                                                              gcgGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTggg                       <  1:4040400/101‑1 (MQ=255)
                                                                               cgGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGa                      >  1:3189712/1‑101 (MQ=255)
                                                                               cgGGGGCGTGGATCGAAACGGACATGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGa                      >  2:1147954/1‑101 (MQ=255)
                                                                               cgGGGGCGGGGATCGAAACGGACCTGAAATGCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGa                      <  2:1849670/101‑1 (MQ=255)
                                                                                gggcgCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGAc                     >  2:160251/1‑101 (MQ=255)
                                                                                 ggggCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCTGCGGGGAGCAC.GACGCTCCCCCCGCGGGGGCGTGGATCGTAACCGCTTTcaccggcct                    >  1:1882481/1‑93 (MQ=37)
                                                                                 ggggCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACt                    >  1:1920714/1‑101 (MQ=255)
                                                                                     cGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGcc                >  1:396221/1‑101 (MQ=255)
                                                                                      gTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCa               <  2:1663943/101‑1 (MQ=255)
                                                                                        ggATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAgc             >  1:1312275/1‑101 (MQ=255)
                                                                                        ggATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAgc             >  1:200545/1‑101 (MQ=255)
                                                                                         gATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAgcc            <  1:1051343/101‑1 (MQ=255)
                                                                                         gATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAgcc            <  1:282457/101‑1 (MQ=255)
                                                                                            tctaaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAgccgcc         <  2:3285960/98‑1 (MQ=255)
                                                                                            cGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAgccgcc         <  1:2399354/101‑1 (MQ=255)
                                                                                             gAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAGccgccc        >  2:420593/1‑101 (MQ=255)
                                                                                             gAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAGccgccc        <  1:2913636/101‑1 (MQ=255)
                                                                                              aaaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAGccgcccg       >  2:2257356/1‑101 (MQ=255)
                                                                                              aaaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAGccgcccc       >  1:1007394/1‑100 (MQ=255)
                                                                                                aCGGACCTGAAATTACGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAGCCGCCcgcg     >  2:189183/1‑101 (MQ=255)
                                                                                                   gACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAGCCGCCCGCGCCa  >  1:1556299/1‑101 (MQ=255)
                                                                                                   gACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAGCCGCCCGCGCCa  >  1:507578/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCGCCCTGCGGGTGATGGCGACGGGCACCACGGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCAC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTCGTGGGACTGGCCAGCCGCCCGCGCCA  >  NC_014034/1396372‑1396571

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG