Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,081,127 C→A 7.1% S463I (AGC→ATC) ‡ PP_2691 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,081,1270CA 7.1% 169.5 / 4.0 56S463I (AGC→ATC) ‡PP_2691hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/28);  new base A (2/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

CCGGGCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAGTCGT  >  NC_002947/3080979‑3081274
                                                                                                                                                    |                                                                                                                                                    
ccGGGCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAgcgcgc                                                                                                                                                      >  1:59394/1‑149 (MQ=255)
           ggTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCActgc                                                                                                                                           <  2:59394/149‑1 (MQ=255)
             ttGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCActgctg                                                                                                                                         <  5:274293/149‑1 (MQ=255)
                 cACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGACTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGata                                                                                                                                     <  3:233305/149‑1 (MQ=255)
                  aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGTGCCATAGCGCGCTGCTGGa                                                                                                                                       >  5:339991/1‑146 (MQ=255)
                  aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGTGCCATAGCGCGCTGCTGGa                                                                                                                                       <  6:339991/146‑1 (MQ=255)
                  aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCActgc                                                                                                                                           >  1:147751/1‑142 (MQ=255)
                  aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCActgc                                                                                                                                           >  6:343849/1‑142 (MQ=255)
                  aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCActgc                                                                                                                                           <  5:343849/142‑1 (MQ=255)
                  aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCActgc                                                                                                                                           <  2:147751/142‑1 (MQ=255)
                  aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGatat                                                                                                                                    >  5:39406/1‑149 (MQ=255)
                     taTTTCACCGTAGTTCCATTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCCATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGc                                                                                                                                 <  4:301163/148‑1 (MQ=255)
                       tttCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCgg                                                                                                                               <  7:274395/149‑1 (MQ=255)
                              gTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGAtgc                                                                                                                        >  4:180732/1‑147 (MQ=255)
                                   tCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGc                                                                                                                   >  4:201609/1‑149 (MQ=255)
                                            ccATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGTGCCATAGCGCGCTGCTGGagcacc                                                                                                                                  <  6:99511/125‑6 (MQ=38)
                                            ccATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGTGCCATAGCGCGCTGCTGGagcacc                                                                                                                                  >  5:99511/1‑120 (MQ=38)
                                             cATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑tt                                                                                                        >  4:250159/1‑149 (MQ=255)
                                             cATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑tt                                                                                                        >  5:233344/1‑149 (MQ=255)
                                               tGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑tttt                                                                                                      >  3:65568/1‑149 (MQ=255)
                                               tGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑tttt                                                                                                      <  3:136723/149‑1 (MQ=255)
                                                 ccTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTc                                                                                                    <  6:233344/149‑1 (MQ=255)
                                                       aggtgtcgTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCggctcc                                                                                              <  5:73631/143‑1 (MQ=255)
                                                       aGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCggctcc                                                                                              >  1:231758/1‑149 (MQ=255)
                                                             cgTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCggctcctggctc                                                                                        <  4:65568/149‑1 (MQ=255)
                                                              gTACCGCGCCCTTGAACATCGGTTGTGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCggctcctggctc                                                                                        <  5:56489/148‑1 (MQ=255)
                                                               tACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGATGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCggctcctggctc                                                                                        >  6:56489/1‑147 (MQ=255)
                                                                    cgcCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTgcg                                                                                 <  3:257786/149‑1 (MQ=255)
                                                                           gAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑ttttt                                                                                                     >  2:163529/1‑122 (MQ=255)
                                                                           gAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑ttttt                                                                                                     <  1:163529/122‑1 (MQ=255)
                                                                           gAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAAt                                                                            <  6:141643/147‑1 (MQ=255)
                                                                            aaCATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGcc                                                                         <  6:39406/149‑1 (MQ=255)
                                                                             aCATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCggctcctggctc                                                                                        <  4:36596/133‑1 (MQ=255)
                                                                             aCATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCggctcctggctc                                                                                        >  3:36596/1‑133 (MQ=255)
                                                                                       ggCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAg                                                                      <  8:281217/141‑1 (MQ=255)
                                                                                       ggCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAg                                                                      >  7:281217/1‑141 (MQ=255)
                                                                                           aTCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtcttc                                                          <  3:250159/149‑1 (MQ=255)
                                                                                            tCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtcttct                                                         >  6:159650/1‑149 (MQ=255)
                                                                                            tCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtcttc                                                          >  4:302492/1‑148 (MQ=255)
                                                                                            tCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtcttc                                                          <  3:302492/148‑1 (MQ=255)
                                                                                                         ctgcAGCCTCAGGGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGc                                            >  4:168556/4‑149 (MQ=255)
                                                                                                          gTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGgcg                                                                   >  1:69559/1‑126 (MQ=255)
                                                                                                          gTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGgcg                                                                   <  2:69559/126‑1 (MQ=255)
                                                                                                             aGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATGCCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTa                                        <  3:180732/149‑1 (MQ=255)
                                                                                                             aGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAg                                             >  1:116401/1‑144 (MQ=255)
                                                                                                             aGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAg                                             <  2:116401/144‑1 (MQ=255)
                                                                                                             aGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTa                                        >  5:164486/1‑149 (MQ=255)
                                                                                                               ccccTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTAc                                       >  1:238710/1‑148 (MQ=255)
                                                                                                                  cTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCag                                    >  4:25331/1‑148 (MQ=255)
                                                                                                                         tGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTg                            >  5:198715/1‑149 (MQ=255)
                                                                                                                           cAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTg                            <  3:25331/147‑1 (MQ=255)
                                                                                                                            aGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTc                         <  6:164486/149‑1 (MQ=255)
                                                                                                                                  aGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATTCCCATCTGCATTTC‑TTTTTCGGATCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGTCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATc                   <  4:251756/149‑1 (MQ=255)
                                                                                                                                  aGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATc                   >  5:261703/1‑149 (MQ=255)
                                                                                                                                     tGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCgg                 >  8:19555/1‑148 (MQ=255)
                                                                                                                                              agcGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGaa       <  6:198715/149‑1 (MQ=255)
                                                                                                                                              agcGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGaa       <  7:19555/149‑1 (MQ=255)
                                                                                                                                                cgcgcTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAGt     >  5:302886/1‑149 (MQ=255)
                                                                                                                                                 gcgcTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAGTc    <  3:168556/149‑1 (MQ=255)
                                                                                                                                                   gcTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAGTCGt  >  5:28948/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGGGCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTC‑TTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAGTCGT  >  NC_002947/3080979‑3081274

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑