Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,890 T→A 2.6% intergenic (+194/+1704) RCAP_rcc01307 → / ← RCAP_rcc01309 hypothetical protein/IS66 family transposase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,8900T→A 2.6% 6.4 ‑0.4 ‑0.1218intergenic (+194/+1704)RCAP_rcc01307/RCAP_rcc01309hypothetical protein/IS66 family transposase
Reads supporting (aligned to +/- strand):  new base (2/4):  ref base (98/114):  total (100/118)
Fisher's exact test strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.08e-01

Only 100 of 231 total aligned reads displayed.

ATCGAAACGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGC  >  NC_014034/1397790‑1397990
                                                                                                     |                                                                                                    
 tCGAAACGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTg                                                                                                     <  2:1474037/101‑1 (MQ=255)
 tCGAAACGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTg                                                                                                     <  1:4675206/101‑1 (MQ=255)
 tCGAAACGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTg                                                                                                     <  1:3880000/101‑1 (MQ=255)
   gAAACGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGtc                                                                                                   >  1:2350042/1‑101 (MQ=255)
   gAAACGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGtc                                                                                                   >  1:3525987/1‑101 (MQ=255)
   gAAACGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGtc                                                                                                   <  2:55919/101‑1 (MQ=255)
    aaaCGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGtct                                                                                                  >  1:1065768/1‑101 (MQ=255)
      aCGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCa                                                                                                <  1:4439105/101‑1 (MQ=255)
        ggTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCAcc                                                                                              <  1:1875640/101‑1 (MQ=255)
        ggTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCAcc                                                                                              >  2:1609060/1‑101 (MQ=255)
          tgtCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCaa                                                                                            >  2:4169281/1‑101 (MQ=255)
          tgtCCAGCGTCACATGCAGCGGCGCCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCGCGACCGCTGTCTCACCaa                                                                                            >  2:1973961/1‑101 (MQ=255)
             ccAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAg                                                                                         >  1:1760657/1‑101 (MQ=255)
             ccAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAg                                                                                         <  2:166677/101‑1 (MQ=255)
               aGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTc                                                                                       <  1:4420514/101‑1 (MQ=255)
                 cGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGc                                                                                     >  1:4976946/1‑101 (MQ=255)
                  gTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCt                                                                                    <  1:1829035/101‑1 (MQ=255)
                   tCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTc                                                                                   >  2:2690142/1‑101 (MQ=255)
                    cacaTGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTcc                                                                                  >  1:2818433/1‑101 (MQ=255)
                    cacaTGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCCCGGGGGCGTGGATCGAAAC.CCGGCACCTAAGCCCCGCCCCCTGTCTCACCAACAGCCGCTcc                                                                                  >  1:3984013/1‑101 (MQ=38)
                     acaTGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTccc                                                                                 >  1:3262795/1‑101 (MQ=255)
                     acaTGCAGCGGCGGCGCATCGCCATGGCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTccc                                                                                 <  1:3333408/101‑1 (MQ=255)
                       aTGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTccccc                                                                               <  2:4386786/101‑1 (MQ=255)
                        tGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTcccccc                                                                              <  1:2011444/101‑1 (MQ=255)
                         gCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTTTCTCACCAACAGTCGCTCCCCCcg                                                                             <  2:3964829/101‑1 (MQ=255)
                         gCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCcg                                                                             <  1:2500196/101‑1 (MQ=255)
                          cAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCCCCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCcgc                                                                            <  1:19229/101‑1 (MQ=255)
                           aGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCcgcg                                                                           >  1:377590/1‑101 (MQ=255)
                           aGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCcgcg                                                                           >  2:4948651/1‑101 (MQ=255)
                              ggcgTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCgggg                                                                        <  1:4876987/101‑1 (MQ=255)
                               gcgTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCggggg                                                                       >  2:4819879/1‑101 (MQ=255)
                                cgccGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGc                                                                      <  2:3462976/98‑1 (MQ=255)
                                cgTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGc                                                                      >  2:103111/1‑101 (MQ=255)
                                 gTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCg                                                                     >  2:2188028/1‑101 (MQ=255)
                                 gTCGCATCGCCATGTCGCCCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCg                                                                     <  1:3449729/101‑1 (MQ=255)
                                  tCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGt                                                                    <  2:2550524/101‑1 (MQ=255)
                                   cGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCggg                                                                   >  2:2884200/1‑99 (MQ=255)
                                   cGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTg                                                                   >  2:1895833/1‑101 (MQ=255)
                                    gCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTgg                                                                  >  2:3714932/1‑101 (MQ=255)
                                     ccatcgagcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGCCAAGGCCGCAGCAGCGGCTGCGGTCGCTCCCCCCGCGGGGGCGTgg                                                                  <  1:4194597/92‑1 (MQ=0)
                                     ccatcgagcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGCCAAGGCCGCAGCAGCGGCTGCGGTCGCTCCCCCCGCGGGGGCGTgg                                                                  >  2:4663277/10‑101 (MQ=0)
                                     cATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                 <  1:69878/101‑1 (MQ=255)
                                      aTCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGGggcg                                                                >  2:3551976/1‑99 (MQ=255)
                                       tCGCCATGTCGCCCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  2:3670259/101‑1 (MQ=255)
                                        cGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              <  2:1708149/101‑1 (MQ=255)
                                        cGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              >  1:1571498/1‑101 (MQ=255)
                                           cATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                           >  2:87090/1‑101 (MQ=255)
                                            aTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCCGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          >  1:547523/1‑101 (MQ=255)
                                            aTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          <  2:2032406/101‑1 (MQ=255)
                                             tGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                         <  2:3061794/101‑1 (MQ=255)
                                              gtccccccccccGCGGGGGCGTGGATCGAACC.CCGGCACCTCAGCCCCGCCCGCTGTCTCCCCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                        <  1:963978/95‑1 (MQ=38)
                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActct                                                       <  1:3936034/101‑1 (MQ=255)
                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTAt                                                     >  2:3899245/1‑101 (MQ=255)
                                                  cTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATg                                                    <  1:1022328/101‑1 (MQ=255)
                                                   tCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgg                                                   >  1:3273607/1‑101 (MQ=255)
                                                   tCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCAGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACCCTATgg                                                   >  1:4275481/1‑101 (MQ=255)
                                                   cccccccGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgg                                                   <  1:3109770/100‑1 (MQ=255)
                                                    ccccccGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATggg                                                  >  2:3357801/1‑101 (MQ=255)
                                                    ccccccGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATggg                                                  <  1:937956/101‑1 (MQ=255)
                                                    ccccccGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATggg                                                  <  2:4939160/101‑1 (MQ=255)
                                                     cccccGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATggga                                                 >  2:2345322/1‑101 (MQ=255)
                                                      ccccGCGGGGGCGTGGATCGAGAC.CCGGCACCTCACCCCCGCCCGCTGTCTCGCCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggag                                                <  1:2199132/101‑1 (MQ=255)
                                                      ccccGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggag                                                <  1:1043226/101‑1 (MQ=255)
                                                          gcgGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggcta                                            <  1:4236849/101‑1 (MQ=255)
                                                           cgGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctag                                           >  1:4756740/1‑101 (MQ=255)
                                                           cgGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctag                                           >  1:3005354/1‑101 (MQ=255)
                                                                 cGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggaggc                                     >  1:583368/1‑101 (MQ=255)
                                                                  gTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggaggct                                    <  2:3342581/101‑1 (MQ=255)
                                                                  gTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggaggct                                    <  1:3385419/101‑1 (MQ=255)
                                                                  gTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTAGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATggggggctagggaggct                                    >  2:3410166/1‑101 (MQ=255)
                                                                   tGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggaggcta                                   <  2:913728/101‑1 (MQ=255)
                                                                   gtgATCGTAAC.CGGGCATCTCAGCCCCGCCCGCTGTCTCACCAACAGCCGCTCCCCCCGCGGGTGCGTGGATCGAAATTCTATgggtggttagggaggcga                                   >  2:804752/3‑99 (MQ=37)
                                                                    ggATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTAc                                  <  2:5029748/101‑1 (MQ=255)
                                                                     gATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACg                                 <  1:3001536/101‑1 (MQ=255)
                                                                      aTCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGc                                <  2:2354848/101‑1 (MQ=255)
                                                                      aTCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGc                                <  2:2070125/101‑1 (MQ=255)
                                                                      aTCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGc                                <  2:4501707/101‑1 (MQ=255)
                                                                       tCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCa                               <  1:4059972/101‑1 (MQ=255)
                                                                        cGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGTGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCAc                              >  1:297448/1‑101 (MQ=255)
                                                                         gAACC.CCGGCCCCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCAcc                             <  1:2140899/101‑1 (MQ=255)
                                                                          aaaC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCt                            >  2:3927258/1‑101 (MQ=255)
                                                                          aaaC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCt                            >  2:4066837/1‑101 (MQ=255)
                                                                           aaC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCtt                           <  2:1571368/101‑1 (MQ=255)
                                                                            aC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTc                          <  1:3602514/101‑1 (MQ=255)
                                                                               ccGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCaa                        >  2:4836633/1‑101 (MQ=255)
                                                                               ccGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCaa                        >  2:134563/1‑101 (MQ=255)
                                                                                cGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCaaa                       <  1:2876786/101‑1 (MQ=255)
                                                                                cGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCaaa                       <  1:4948981/101‑1 (MQ=255)
                                                                                 ggCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCaaaa                      <  2:4741332/101‑1 (MQ=255)
                                                                                   cACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGt                    <  1:2385853/101‑1 (MQ=255)
                                                                                   cACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGt                    <  2:954841/101‑1 (MQ=255)
                                                                                      cTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGc                 >  1:1293423/1‑101 (MQ=255)
                                                                                       tCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCt                >  2:4283044/1‑101 (MQ=255)
                                                                                         aGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGAACCTTCAAAAGTCGATcc              >  2:538843/1‑101 (MQ=255)
                                                                                           ccccGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTcccc            <  2:2203671/101‑1 (MQ=255)
                                                                                            cccgcccgCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTccccc           <  1:1688239/101‑1 (MQ=255)
                                                                                             ccgcccgcTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTcccccc          <  2:4423481/101‑1 (MQ=255)
                                                                                              cgcccgcTGTCTCACCAACAGTCGCTCCCCCCGTGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCcg         <  1:4309877/101‑1 (MQ=255)
                                                                                              cgcccgcTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCcg         <  1:4283361/101‑1 (MQ=255)
                                                                                                cccgcTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCcgcg       <  2:1088259/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
ATCGAAACGGTGTCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCACCTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGC  >  NC_014034/1397790‑1397990

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