Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 478,984 T→C 5.4% Q162R (CAG→CGG) ‡ cca ← multifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947478,9840TC 5.4% 218.5 / 2.5 74Q162R (CAG→CGG) ‡ccamultifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base C (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CGAGCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCC  >  NC_002947/478838‑479121
                                                                                                                                                    |                                                                                                                                         
cGAGCTCAGGCAGCAGCTCCTGCAGCCCCCCACACGCGCTCAGCACATTGAAGAATACCTGCCGTTGCACTTCCAACAGCGCGCGCTCGATTACCTTCCAGCTGCGCTCGGCAGTCAGTCCATGCAGTTCGCCTCAGGCGATGATCTgg                                                                                                                                         >  4:316109/1‑149 (MQ=255)
 gAGCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGTACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGgc                                                                                                                                        <  7:184480/149‑1 (MQ=255)
 gAGCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGgc                                                                                                                                        <  1:208443/149‑1 (MQ=255)
 gAGCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGgc                                                                                                                                        <  3:316109/149‑1 (MQ=255)
  aGCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGgc                                                                                                                                        >  1:96262/1‑149 (MQ=255)
   gCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGgcgc                                                                                                                                      <  5:24959/149‑1 (MQ=255)
   gCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGgcg                                                                                                                                       >  6:24959/1‑148 (MQ=255)
       aGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCa                                                                                                                                  <  4:136940/149‑1 (MQ=255)
       aGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCa                                                                                                                                  <  4:142699/149‑1 (MQ=255)
         gcagcaGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAgc                                                                                                                                >  4:200772/1‑149 (MQ=255)
           agcagcTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATGCGCCAGATCAgcgc                                                                                                                              <  4:158134/149‑1 (MQ=255)
           agcagcTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATGCGCCAGATCAgcgc                                                                                                                              <  3:216964/149‑1 (MQ=255)
            gcagcTCCTGCCGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTGCATCAGAGCGCGCTGGATTTCCTTCCAGCTGCGCTCGGGAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAgcgc                                                                                                                              >  3:247521/1‑148 (MQ=255)
               gcTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAg                                                                                                                          <  2:85431/149‑1 (MQ=255)
                 tCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGGACTTCCATCAGCGCGCGCTCGCTTTCCTTCCACCTGCGCTCGGGAGTGAGTGCCTGCAGTTCGCCTGCGGCGCTGATCTGGCGCCTCAGCGCCAAGGt                                                                                                                        >  1:286536/1‑149 (MQ=255)
                  ccTGCAGGGCCCCCCCCGCGCTCAGCACCTAGAAGAATACCGGCGGTTGCACGTCCATCCGCGCGCGCGCGCTTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGtt                                                                                                                       <  8:144655/149‑1 (MQ=255)
                       aGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGTGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATc                                                                                                                  >  1:254070/1‑149 (MQ=255)
                          gCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCgg                                                                                                                >  4:166849/1‑148 (MQ=255)
                           ccccACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGc                                                                                                               >  8:200564/1‑148 (MQ=255)
                            cccACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCaa                                                                                                             <  8:275702/149‑1 (MQ=255)
                                 acGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCg                                                                                                         >  2:23916/1‑148 (MQ=255)
                                      cTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAAccc                                                                                                   <  5:233776/149‑1 (MQ=255)
                                           cACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCCGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGTCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAgt                                                                                               <  8:244842/148‑1 (MQ=255)
                                           cACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAgtg                                                                                              >  8:231627/1‑149 (MQ=255)
                                           cACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCGCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCGGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAgtg                                                                                              <  2:286536/149‑1 (MQ=255)
                                              cTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAgtggtg                                                                                           >  7:110423/1‑149 (MQ=255)
                                               ttGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGc                                                                                          >  6:290748/1‑149 (MQ=255)
                                                           cTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGgcag                                                                               >  5:131984/1‑149 (MQ=255)
                                                            tGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGgcag                                                                               >  2:336972/1‑148 (MQ=255)
                                                              cGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGgcagca                                                                             <  2:39455/148‑1 (MQ=255)
                                                               ggTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCaaa                                                                           <  7:200564/149‑1 (MQ=255)
                                                                     aCTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTg                                                                      >  2:246805/1‑148 (MQ=255)
                                                                     aCTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGc                                                                     <  3:200772/149‑1 (MQ=255)
                                                                     aCTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGAGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCCGCGCCAAGGTTTCATCGGCAACCCGGAACCCGAGTGGTGCCTAACGGGCAGCAAAGCGtt                                                                      >  1:62898/1‑147 (MQ=255)
                                                                         ccATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAAc                                                                  <  6:131984/148‑1 (MQ=255)
                                                                          cATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTACGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAAcgc                                                                >  1:301746/1‑149 (MQ=255)
                                                                           aTCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAAcgcg                                                               >  3:139724/1‑149 (MQ=255)
                                                                            tCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCGAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTg                                                                      >  1:69015/1‑142 (MQ=255)
                                                                            tCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTg                                                                      <  2:69015/141‑1 (MQ=255)
                                                                            tCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACgcgc                                                              <  1:246805/149‑1 (MQ=255)
                                                                              aGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCg                                                                                         >  2:382925/1‑120 (MQ=255)
                                                                               gcgcgcgcTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCgcag                                                            >  6:133766/1‑148 (MQ=255)
                                                                                 gcgcgcTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCAc                                                         >  5:89439/1‑149 (MQ=255)
                                                                                               cttccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCtt                                           >  5:340271/1‑149 (MQ=255)
                                                                                                ttccAGCTGCGCTCGGCCCTGAGTGCCTGCAGTTCGCCGGAGGCGCGTCTCTGGCGCATCAGCGCCTAGGTTTCATCGGCAACCCGGCACCCCAGTGGTGCGTAACGGGCAGCAAAGCGGGCAACGCGCAGCACCCGCAAAGGATCTTc                                          <  4:247521/149‑1 (MQ=255)
                                                                                                ttccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGt                                                                       <  3:220202/120‑1 (MQ=255)
                                                                                                ttccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGt                                                                       >  4:220202/1‑120 (MQ=255)
                                                                                                ttccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTc                                          >  5:4745/1‑149 (MQ=255)
                                                                                                  ccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTc                                          <  1:62636/147‑1 (MQ=255)
                                                                                                  ccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTc                                          <  2:96262/147‑1 (MQ=255)
                                                                                                          gcTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGGGGTGCGTAACGGGCAGCAAAGCGTGCACCGCGCAGCACCCGCCAAGGATCTTCGGCGCATGcc                                >  2:233189/1‑149 (MQ=255)
                                                                                                           cTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCg                               <  5:290748/149‑1 (MQ=255)
                                                                                                           cTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCg                               <  6:4745/149‑1 (MQ=255)
                                                                                                            tCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCgg                              >  1:382118/1‑149 (MQ=255)
                                                                                                              ggCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCggg                             <  2:20012/148‑1 (MQ=255)
                                                                                                                       gtgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTACCGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCATAGGATCTTCGGCGAATGCCGGGGAAACATg                     <  1:233189/148‑1 (MQ=255)
                                                                                                                        tgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCCGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGc                   <  1:336972/149‑1 (MQ=255)
                                                                                                                          ccTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGATACATGGCg                  >  4:341659/1‑148 (MQ=255)
                                                                                                                          ccTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGt                 <  2:254070/149‑1 (MQ=255)
                                                                                                                          ccTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGt                 <  3:166849/149‑1 (MQ=255)
                                                                                                                           cTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGa                           <  7:318166/138‑1 (MQ=255)
                                                                                                                           cTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGa                           >  8:318166/1‑138 (MQ=255)
                                                                                                                              cAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGc              <  6:89439/148‑1 (MQ=255)
                                                                                                                                 ttCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAAt                                   <  6:118796/124‑1 (MQ=255)
                                                                                                                                 ttCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGc          >  7:207823/1‑149 (MQ=255)
                                                                                                                                     ccTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGtt       >  2:370425/1‑148 (MQ=255)
                                                                                                                                      cTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGAt                                              >  2:146982/1‑108 (MQ=255)
                                                                                                                                      cTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGAt                                              <  1:146982/108‑1 (MQ=255)
                                                                                                                                      cTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGa     <  4:139724/149‑1 (MQ=255)
                                                                                                                                        ttggCGCTGATGCGCCAGATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATc   >  3:21119/3‑149 (MQ=255)
                                                                                                                                         aGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATcc  <  5:316478/149‑1 (MQ=255)
                                                                                                                                          ggCGCTGATGCGCCAGATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATc   >  3:196850/1‑147 (MQ=255)
                                                                                                                                             gcTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGaa                                                                                                      >  2:177104/1‑45 (MQ=255)
                                                                                                                                             gcTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGaa                                                                                                      <  1:177104/45‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
CGAGCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCC  >  NC_002947/478838‑479121

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: