Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 478,986 G→C 6.0% R161R (CGC→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,9860GC 6.0% 218.7 / 4.1 67R161R (CGC→CGGccamultifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base C (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

CGAGCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCCAGATCGTTCTTG  >  NC_002947/478838‑479133
                                                                                                                                                      |                                                                                                                                                   
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                                                      <  2:96262/147‑1 (MQ=255)
                                                                                                  ccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTc                                                      <  1:62636/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                             <  3:166849/149‑1 (MQ=255)
                                                                                                                          ccTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGt                             <  2:254070/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)
                                                                                                                                                     ggCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCCAGATCGTTCTTg  >  6:123574/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGAGCTCAGGCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCCAGATCGTTCTTG  >  NC_002947/478838‑479133

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: