Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,495 C→G 13.2% T459T (ACC→TCG RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,4950C→G13.2% 93.0 ‑0.5 ‑0.5313T459T (ACC→ACGRCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (23/18):  ref base (138/132):  total (163/150)
Fisher's exact test strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.42e-01

Only 100 of 351 total aligned reads displayed.

GAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCG  >  NC_014034/1397395‑1397594
                                                                                                           |                                                                                                   
gAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca...cc                                                                                                     <  1:1912383/101‑1 (MQ=255)
  aaCTCTTCAATGGATATGGCAGAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca...ccac                                                                                                   >  1:3292526/1‑101 (MQ=255)
  aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca...ccac                                                                                                   >  2:2928224/1‑101 (MQ=255)
  aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCACCCTGca...ccac                                                                                                   >  2:4339955/1‑101 (MQ=255)
    ctctTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCg                                                                                                 <  1:4452120/101‑1 (MQ=255)
     tctTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCgg                                                                                                <  2:4566918/101‑1 (MQ=255)
      ctTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGt                                                                                               >  2:236964/1‑101 (MQ=255)
         cAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGc                                                                                            >  2:5773012/1‑101 (MQ=255)
          aaTGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCt                                                                                           >  2:4658635/1‑101 (MQ=255)
          aaTGGATATGACAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCt                                                                                           <  2:1172391/101‑1 (MQ=255)
             ggATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTccc                                                                                        >  2:2467511/1‑101 (MQ=255)
             ggATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTccc                                                                                        >  2:2020451/1‑101 (MQ=255)
             ggATATGGCAAAACCGTCCGCGACATGTCGCTCCACCCGC.GGGGGCGCGGATGGAAACACCG.TCATC.ACCACGCAA.CCGCCTTGCA...CCACCGGTCGCTccc                                                                                        >  1:3275496/1‑101 (MQ=255)
               atatGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTccccc                                                                                      <  2:823890/101‑1 (MQ=255)
                tatGGCAAACCCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTcccccc                                                                                     <  2:3935351/101‑1 (MQ=255)
                   ggCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCcgcg                                                                                  <  2:3793612/101‑1 (MQ=255)
                      aaaaCCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCgggg                                                                               >  1:2056491/1‑101 (MQ=255)
                      aaaaCCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCgggg                                                                               >  1:1969433/1‑101 (MQ=255)
                      aaaaCCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCgggg                                                                               >  2:1681765/1‑101 (MQ=255)
                       aaaCCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCggggg                                                                              <  2:2186609/101‑1 (MQ=255)
                       aaaCCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCggggg                                                                              >  2:5565104/1‑101 (MQ=255)
                        aaCCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGc                                                                             >  1:3512701/1‑101 (MQ=255)
                        aaCCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGc                                                                             <  1:1049347/101‑1 (MQ=255)
                         aCCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCg                                                                            >  1:2187944/1‑101 (MQ=255)
                            gTCCGCGACATGTCGCTCCCCCCGCGGGGGGCGTGGATCGAAACCCCG.TCATC.ACCACGCAA.CCGCCCTGCA...CAACCGGTCGCCCCCCCCCGggggggggg                                                                          >  1:1565255/1‑97 (MQ=25)
                            gTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTgg                                                                         <  1:2262689/101‑1 (MQ=255)
                            gTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTgg                                                                         <  1:5067364/101‑1 (MQ=255)
                              cgataccgcgTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCT...CGATC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                         >  1:5789188/10‑101 (MQ=0)
                              ccGCGGCATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                       <  2:5690091/101‑1 (MQ=255)
                               gataccgcgTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCT...CGATC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                        <  2:1406599/93‑1 (MQ=0)
                               cgcgACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                      <  2:1257351/101‑1 (MQ=255)
                                ataccgcgTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCT...CGATC.GTCGCTCCCCCCGCGGGGGCGGGGAt                                                                       >  1:2963205/8‑101 (MQ=0)
                                 taccgcgTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCT...CGATC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                      >  1:1777458/7‑101 (MQ=2)
                                  gACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                   <  2:5233409/101‑1 (MQ=255)
                                   aCATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                  >  1:4178730/1‑101 (MQ=255)
                                    cgcgTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCT...CGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                   >  1:2493458/4‑101 (MQ=2)
                                    cctGTCGCCGCCCCCGC.GTGGGGGAGGGACGAAACACCT.TCATC.GCCCCGCAA.CCGCCGTGCC...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                 <  2:922023/99‑1 (MQ=16)
                                    cATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                 <  2:4116071/101‑1 (MQ=255)
                                     aTGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaaca                                                                >  2:3144568/1‑101 (MQ=255)
                                     aTGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaaca                                                                >  2:5083007/1‑101 (MQ=255)
                                      cgTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCCGCCGGGCG...CGATC.GCCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                 >  2:5197803/2‑101 (MQ=2)
                                        tCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAcc                                                             >  1:97394/1‑101 (MQ=255)
                                        tCGCTCCCCCCGC.GGGGGCGTGGATCGAAAC.CCG.CCA.C.GCCTTGGAA.CAGCGCGCCATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGCGGGTCGAaaccac                                                              >  1:804522/1‑101 (MQ=0)
                                           ccccccccGC.GGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCT...CGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgt                                                            <  2:1899379/99‑6 (MQ=2)
                                           cTCCCCCCGC.GGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCT...CGATC.GTCGCTACCCCCGCGGGGGCGTGGATCGAcacggagt                                                            >  1:3616779/1‑96 (MQ=2)
                                            tCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGAt                                                         >  1:4420743/1‑101 (MQ=255)
                                             ccccccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATc                                                        <  1:5595589/101‑1 (MQ=255)
                                              cccccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                                                       >  1:3822330/1‑101 (MQ=255)
                                              cccccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                                                       >  1:5469949/1‑101 (MQ=255)
                                              cccccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                                                       <  1:5096106/101‑1 (MQ=255)
                                               ccccGC.GGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCT...CGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccag                                                        <  1:1663365/101‑10 (MQ=2)
                                               ccccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCAt                                                      <  2:3292326/101‑1 (MQ=255)
                                               ccccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCAt                                                      <  2:141552/101‑1 (MQ=255)
                                                cccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                     <  1:667877/101‑1 (MQ=255)
                                                cccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                     <  1:4831626/101‑1 (MQ=255)
                                                cccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                     <  2:1769706/101‑1 (MQ=255)
                                                cccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                     >  2:302400/1‑101 (MQ=255)
                                                 ccGC.GGGGGCGTGGATCGAAACACCG.TCATT.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATgg                                                    <  1:5250688/101‑1 (MQ=255)
                                                 ccGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATgg                                                    >  1:268016/1‑101 (MQ=255)
                                                   gc.gGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGcc                                                  <  1:1256653/101‑1 (MQ=255)
                                                   gc.gGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGcc                                                  <  1:776627/101‑1 (MQ=255)
                                                    c.gGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGccc                                                 <  2:3362940/101‑1 (MQ=255)
                                                    c.gGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGccc                                                 >  1:1134436/1‑101 (MQ=255)
                                                      gggggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATAATGGCCCg                                                >  1:892680/1‑101 (MQ=255)
                                                       ggggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGc                                               <  2:3999705/101‑1 (MQ=255)
                                                         ggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTTGATCGAAACAACCCGATCATGGCCCGCCt                                             >  1:3852407/1‑101 (MQ=255)
                                                          gCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTc                                            <  1:4104137/101‑1 (MQ=255)
                                                           cGTGGATCCAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCa                                           <  1:3395240/101‑1 (MQ=255)
                                                              ggATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgcg                                        <  1:1385918/101‑1 (MQ=255)
                                                               gATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgcgc                                       <  2:3433266/101‑1 (MQ=255)
                                                               gATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgcgc                                       >  1:3710900/1‑101 (MQ=255)
                                                               gATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgcgc                                       >  2:25384/1‑101 (MQ=255)
                                                                aTCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgcc                                      >  1:1714358/1‑101 (MQ=255)
                                                                 tCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgccg                                     <  1:1206052/101‑1 (MQ=255)
                                                                  cNNNNCACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCNAAACAACCCGATCATGGCCCGCCTCAGcgccgc                                    >  2:3212456/1‑101 (MQ=255)
                                                                   gAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCcgcg                                   <  1:2379563/101‑1 (MQ=255)
                                                                    aaaCACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCgcgc                                  >  2:74260/1‑101 (MQ=255)
                                                                       cacCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGa                               >  2:3625373/1‑101 (MQ=255)
                                                                       cacCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGa                               <  2:263783/101‑1 (MQ=255)
                                                                       cacCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGa                               >  1:2859937/1‑101 (MQ=255)
                                                                        cccG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGAc                              <  1:3801799/100‑1 (MQ=255)
                                                                          cG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGt                            >  2:5025190/1‑101 (MQ=255)
                                                                           g.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACgcc                           >  2:2835101/1‑99 (MQ=255)
                                                                               atc.aCCACGCAA.CCGCCCTGCA...CCCCCGGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCt                        <  1:1481508/101‑1 (MQ=255)
                                                                               atc.aCCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCt                        >  1:2607650/1‑101 (MQ=255)
                                                                                tc.aCCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTc                       <  2:2056374/101‑1 (MQ=255)
                                                                                cc.accacGCAA.CCCCCCGGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTc                       <  2:5481547/100‑1 (MQ=255)
                                                                                      acGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcccccc                  >  2:2313764/1‑101 (MQ=255)
                                                                                      acGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcccccc                  <  1:5364708/101‑1 (MQ=255)
                                                                                        gCAA.CCGCCCTGCC...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgc                <  2:3729998/101‑1 (MQ=255)
                                                                                         cAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgcg               <  1:4773674/101‑1 (MQ=255)
                                                                                          aa.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgg              >  2:5380486/1‑101 (MQ=255)
                                                                                          aa.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCAGCgg              <  2:4897496/101‑1 (MQ=255)
                                                                                          aa.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCCCCCCCCGCgg              >  2:4718922/1‑101 (MQ=255)
                                                                                              cGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCggggg           >  1:2345018/1‑101 (MQ=255)
                                                                                               gCCCGGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGc          <  1:3957730/101‑1 (MQ=255)
                                                                                                cccTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCg         >  2:5330201/1‑101 (MQ=255)
                                                                                                  cTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg       <  2:127250/101‑1 (MQ=255)
                                                                                                  cTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg       <  2:4245535/101‑1 (MQ=255)
                                                                                                  cTGCA...CAACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGGTCATGGCCCGCCTCAGTGCACCGCTGACGTCGCCCCCCCCCCGGGGGCggg       >  2:157671/1‑99 (MQ=25)
                                                                                                           |                                                                                                   
GAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGC.GGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCA...CCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCG  >  NC_014034/1397395‑1397594

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