Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,545 C→A 5.3% A476D (GCC→GAC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,5450C→A 5.3% 53.6 ‑2.4 ‑0.0568A476D (GCC→GAC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/19):  ref base (276/249):  total (292/276)
Fisher's exact test strand distribution p-value = 9.49e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 600 total aligned reads displayed.

CGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCG  >  NC_014034/1397445‑1397645
                                                                                                    |                                                                                                       
cgcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc                                                                                                         <  1:8509675/101‑1 (MQ=255)
 gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c                                                                                                       <  2:6160776/101‑1 (MQ=255)
    ggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGc                                                                                                  >  2:7642827/1‑101 (MQ=255)
     gggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCGCCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGcc                                                                                                 <  2:1299458/101‑1 (MQ=255)
       gCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTc                                                                                               >  1:3790639/1‑101 (MQ=255)
        cGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCa                                                                                              <  2:7218210/101‑1 (MQ=255)
         gTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAg                                                                                             >  1:1674260/1‑101 (MQ=255)
           ggATCGAAACACCGTCATCCCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgcg                                                                                           <  1:1357456/101‑1 (MQ=255)
           ggATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgcg                                                                                           <  1:8262569/101‑1 (MQ=255)
           ggATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgcg                                                                                           >  1:6584254/1‑101 (MQ=255)
                  aaCACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCt                                                                                    >  1:4682910/1‑101 (MQ=255)
                   acacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCCCAGCGCCgcgccg                                                                                   >  2:2778097/1‑99 (MQ=255)
                    cacCGGCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTgc                                                                                  >  1:8845799/1‑100 (MQ=255)
                     cccGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGAc                                                                                 <  2:163760/100‑1 (MQ=255)
                     acCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGAc                                                                                 <  2:1193589/101‑1 (MQ=255)
                      ccGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACg                                                                                >  1:3674345/1‑101 (MQ=255)
                        gTCATCCCCACGCAACCGCCGTGCACCGCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTc                                                                              <  2:8478439/101‑1 (MQ=255)
                        gTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTc                                                                              <  1:2532663/101‑1 (MQ=255)
                          catcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGc                                                                            >  1:4847879/1‑101 (MQ=255)
                             ccccacGCAACCGCCCTGCACCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcc                                                                         <  2:1479124/99‑1 (MQ=255)
                              accacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccc                                                                        <  1:7498079/101‑1 (MQ=255)
                                 ccGCAACCGCCCTGCCCCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                                     <  2:7075912/100‑1 (MQ=255)
                                 ccGCAACCGCCCTGCACCTCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                                     <  2:5065247/100‑1 (MQ=255)
                                 ccGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                                     <  2:1383848/100‑1 (MQ=255)
                                 acGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                                     <  1:7929186/101‑1 (MQ=255)
                                    cAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcccg                                                                  >  1:6009572/1‑101 (MQ=255)
                                    cAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..AGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgcg                                                                  >  2:5899784/1‑101 (MQ=255)
                                     aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgg                                                                 >  1:962619/1‑101 (MQ=255)
                                     aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgg                                                                 >  1:6357932/1‑101 (MQ=255)
                                       ccGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgggg                                                               <  1:5305757/101‑1 (MQ=255)
                                        cGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCggggg                                                              <  1:9134068/101‑1 (MQ=255)
                                          cccTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCg                                                            <  1:2718576/101‑1 (MQ=255)
                                          cccTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCg                                                            >  1:7099717/1‑101 (MQ=255)
                                             ggaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGa                                                        <  2:8689713/91‑1 (MQ=2)
                                                accaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      >  1:5023549/1‑101 (MQ=255)
                                                 ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                     <  2:5322311/101‑1 (MQ=255)
                                                   accGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                   <  2:1049737/101‑1 (MQ=255)
                                                   accGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                   <  1:4036974/101‑1 (MQ=255)
                                                     gcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                  <  1:2464223/99‑1 (MQ=2)
                                                     atgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcg                                                <  2:4353875/99‑2 (MQ=2)
                                                      ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCGGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                >  2:5230697/1‑101 (MQ=255)
                                                      ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                >  2:8806686/1‑101 (MQ=255)
                                                      ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                <  2:2593351/101‑1 (MQ=255)
                                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATGGAAACAt                                               <  1:8680354/101‑1 (MQ=2)
                                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                               <  1:6876012/101‑1 (MQ=2)
                                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                               >  2:6447517/1‑101 (MQ=2)
                                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                               <  2:5637915/101‑1 (MQ=2)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCCCGGGGGCGCGGATCGaaaaggggg                                            >  2:178394/1‑95 (MQ=2)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                             >  2:6927780/1‑100 (MQ=2)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                             <  1:46087/101‑2 (MQ=2)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                              <  1:5591827/101‑3 (MQ=2)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcc                                             >  1:3516399/1‑101 (MQ=255)
                                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                           <  2:6738245/101‑7 (MQ=2)
                                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                           <  1:7292519/101‑7 (MQ=2)
                                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcga                                            >  1:2012761/1‑99 (MQ=2)
                                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                             >  1:2993132/1‑98 (MQ=2)
                                                            tCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCCCGGGGGCGGGGATCGAAAcggggcgaa                                         >  1:1503410/1‑93 (MQ=2)
                                                            tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatg                                          >  2:511441/1‑97 (MQ=2)
                                                            tCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgcg                                          >  2:7758192/1‑101 (MQ=255)
                                                             ccccccGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaag                                        <  2:2017199/101‑10 (MQ=2)
                                                             ccccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaagc                                          >  1:3468573/1‑95 (MQ=0)
                                                             ccccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCgg                                         >  2:8063310/1‑101 (MQ=255)
                                                              cccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGCCGCTCCCCCCGCGGGGGCGCGGATCGAAACATcggtgac                                        >  2:5759486/1‑95 (MQ=2)
                                                               ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatc                                       <  2:7515841/101‑8 (MQ=2)
                                                               ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatc                                       <  1:7819688/101‑8 (MQ=2)
                                                               ccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAt                                       <  1:5440802/101‑1 (MQ=255)
                                                                cccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatca                                      <  2:737448/101‑9 (MQ=0)
                                                                 ccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCg                                     <  2:8226873/101‑1 (MQ=255)
                                                                     gggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGt                                 <  2:9040339/101‑1 (MQ=255)
                                                                     gggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGt                                 >  1:222536/1‑101 (MQ=255)
                                                                      ggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                                >  2:3869220/1‑101 (MQ=255)
                                                                      ggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                                >  2:323507/1‑101 (MQ=255)
                                                                        ggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAc                              <  2:4701243/101‑1 (MQ=255)
                                                                        ggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAc                              <  1:1229637/101‑1 (MQ=255)
                                                                        ggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAc                              >  1:1824953/1‑101 (MQ=255)
                                                                          cGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAccc                            <  2:7929976/101‑1 (MQ=255)
                                                                           gTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                           >  1:5909767/1‑101 (MQ=255)
                                                                           gTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                           >  2:3538027/1‑101 (MQ=255)
                                                                            tGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgc                          >  1:4424096/1‑101 (MQ=255)
                                                                            tGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgc                          >  2:5707747/1‑101 (MQ=255)
                                                                              gATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccc                        >  1:6654968/1‑101 (MQ=255)
                                                                               aTCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                       <  2:1921761/101‑1 (MQ=255)
                                                                               aTCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                       >  2:3789171/1‑101 (MQ=255)
                                                                               aTCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                       >  2:4893927/1‑101 (MQ=255)
                                                                                 cGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccgc                     >  1:5157100/1‑101 (MQ=255)
                                                                                  gAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACccgccccgcc                    <  2:3438078/101‑1 (MQ=255)
                                                                                   aaaCAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc                   >  1:4310151/1‑101 (MQ=255)
                                                                                   aaaCAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc                   >  1:1037424/1‑101 (MQ=255)
                                                                                      caccccGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc                <  1:8047280/98‑1 (MQ=255)
                                                                                      caacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGGGCGGCCCGCCCCCCCCAtg                >  2:4224299/1‑100 (MQ=37)
                                                                                       acccggagcATGGC.C..CGGCTCAGCTCCGCGCTGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCATCCGCGGATCGGCGTGGCCCGCCCCGCCCATCg               <  2:5523785/93‑1 (MQ=25)
                                                                                       aacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCg               <  2:6654708/101‑1 (MQ=255)
                                                                                         cccGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGt             >  1:2880477/1‑101 (MQ=255)
                                                                                           cGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCg           >  1:7677695/1‑101 (MQ=255)
                                                                                             aTCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCt         >  1:5676605/1‑101 (MQ=255)
                                                                                              tCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTc        >  2:3076961/1‑101 (MQ=255)
                                                                                               cATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCCTCGGTCCCTcc       >  2:2322565/1‑101 (MQ=255)
                                                                                                aTGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTccc      >  1:6237215/1‑101 (MQ=255)
                                                                                                  ggC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTccccc    >  2:7874866/1‑101 (MQ=255)
                                                                                                   gC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTcccccc   >  2:456388/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
CGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCG  >  NC_014034/1397445‑1397645

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