Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 478,983 C→G 10.5% Q162H (CAG→CAC) ‡ 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,9830CG10.5% 229.8 / 13.6 76Q162H (CAG→CAC) ‡ccamultifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base G (6/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCCAGATCGTTC  >  NC_002947/478847‑479130
                                                                                                                                        |                                                                                                                                                   
gcagcaGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAgc                                                                                                                                         >  8:290853/1‑149 (MQ=255)
gcagcaGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAgc                                                                                                                                         >  5:197078/1‑149 (MQ=255)
               gggCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCg                                                                                                                          <  7:290853/149‑1 (MQ=255)
               gggCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCg                                                                                                                          <  6:197078/149‑1 (MQ=255)
                ggCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGAGTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTCAGGCGCGGATCTGGCGCATCAGCGCCAAGCTTTCATctg                                                                                                                         >  3:230330/1‑147 (MQ=255)
                      acacGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGttt                                                                                                                               <  4:339528/137‑1 (MQ=255)
                      acacGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGttt                                                                                                                               >  3:339528/1‑137 (MQ=255)
                      acacGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAAccc                                                                                                                   <  7:278914/149‑1 (MQ=255)
                       cacGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCg                                                                                                                  <  4:49623/149‑1 (MQ=255)
                              tCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAgc                                                                                                                                         >  8:1388/1‑119 (MQ=255)
                              tCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAgc                                                                                                                                         <  7:1388/119‑1 (MQ=255)
                                aGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCTCTCGCTCTATTTCCTTGCAGCTGCGCTCGGCAGTGAGTGCCTGTAGTTCGCCTGAGGCTCTGATCTGGCGCATCAGAGCCAAGGTTTCTTCGTTTACCTGTAACCCCAg                                                                                                         <  4:144657/149‑1 (MQ=255)
                                aGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCa                                                                                                          <  6:142745/148‑1 (MQ=255)
                                  cACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGt                                                                                                                                 >  6:81737/1‑123 (MQ=255)
                                  cACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGt                                                                                                                                 <  5:81737/123‑1 (MQ=255)
                                  cACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAgtg                                                                                                       >  1:185485/1‑149 (MQ=255)
                                        gaagaaTACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGt                                                                                                 >  2:200479/1‑148 (MQ=255)
                                         aagaaTACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGt                                                                                                 >  8:107390/1‑148 (MQ=255)
                                           gaaTACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAAc                                                                                              >  1:302868/1‑149 (MQ=255)
                                           gaaTACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATAGGCAACCCGGAACCCCAGTGGTGCGTAAc                                                                                              >  3:7897/1‑149 (MQ=255)
                                                ccTGGGGTTGCTCTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGgc                                                                                          <  4:230330/148‑1 (MQ=255)
                                                ccTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGgc                                                                                          <  7:107390/148‑1 (MQ=255)
                                                ccTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGgca                                                                                         <  1:200479/148‑1 (MQ=255)
                                                 cTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAgtg                                                                                                       <  2:203948/134‑1 (MQ=255)
                                                     ggTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCaaa                                                                                    <  7:12158/149‑1 (MQ=255)
                                                         gCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCaaa                                                                                    >  4:71726/1‑145 (MQ=255)
                                                         gCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCaaa                                                                                    <  3:71726/145‑1 (MQ=255)
                                                           aCTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTg                                                                               >  6:328879/1‑148 (MQ=255)
                                                                 aTCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAAcgcg                                                                        >  6:71680/1‑149 (MQ=255)
                                                                  tCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAAcgcg                                                                        >  4:158572/1‑148 (MQ=255)
                                                                          cgcTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCg                                                               <  2:185485/149‑1 (MQ=255)
                                                                                 tttCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGa                                                        >  2:324489/1‑149 (MQ=255)
                                                                                     cttccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCtt                                                    >  4:73001/1‑149 (MQ=255)
                                                                                     cttccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGAACtt                                                    >  2:119041/1‑149 (MQ=255)
                                                                                      ttccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTc                                                   >  8:129304/1‑149 (MQ=255)
                                                                                      ttccAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTc                                                   >  5:262231/1‑149 (MQ=255)
                                                                                          aGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGc                                                >  4:22736/1‑148 (MQ=255)
                                                                                          aGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGc                                                >  5:307096/1‑148 (MQ=255)
                                                                                             tGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGgcag                                                                                        <  3:347082/105‑1 (MQ=255)
                                                                                             tGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGgcag                                                                                        >  4:347082/1‑105 (MQ=255)
                                                                                             tGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAAt                                            >  7:200733/1‑149 (MQ=255)
                                                                                              gcgcTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATg                                           >  6:24919/1‑149 (MQ=255)
                                                                                                gcTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGcc                                         <  3:158572/149‑1 (MQ=255)
                                                                                                  tCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCgg                                       >  8:76503/1‑149 (MQ=255)
                                                                                                   cGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCggg                                      >  7:249091/1‑149 (MQ=255)
                                                                                                     gCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGa                                    <  8:249091/149‑1 (MQ=255)
                                                                                                       agtgagtgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGaaa                                  <  7:129304/149‑1 (MQ=255)
                                                                                                         tgagtgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACa                                >  5:84803/1‑149 (MQ=255)
                                                                                                           agtgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATg                              <  2:302868/149‑1 (MQ=255)
                                                                                                            gtgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCgcag                                                                     >  3:151934/1‑109 (MQ=255)
                                                                                                            gtgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCgcag                                                                     <  4:151934/109‑1 (MQ=255)
                                                                                                            gtgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCg                                                  <  3:24907/128‑1 (MQ=255)
                                                                                                            gtgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCg                                                  >  4:24907/1‑128 (MQ=255)
                                                                                                            gtgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGAACTTCGGCGAATGCCGGGGAAACATgg                             <  1:119041/149‑1 (MQ=255)
                                                                                                             tgCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGc                            <  1:324489/149‑1 (MQ=255)
                                                                                                                   cAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAg                        >  3:262037/1‑147 (MQ=255)
                                                                                                                   cAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCa                      <  5:71680/149‑1 (MQ=255)
                                                                                                                         gCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGt                 <  6:307096/148‑1 (MQ=255)
                                                                                                                         gCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGt                 <  8:200733/148‑1 (MQ=255)
                                                                                                                           cTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATg                              <  3:176010/133‑1 (MQ=255)
                                                                                                                           cTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATg                              >  4:176010/1‑133 (MQ=255)
                                                                                                                           cTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGa              <  6:262231/149‑1 (MQ=255)
                                                                                                                            cttggCGCTGATGCGCCAGATCAGCGCCTAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCgca                                                                      <  5:263238/89‑1 (MQ=37)
                                                                                                                            cttggCGCTGATGCGCCAGATCAGCGCCTAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCgca                                                                      >  6:263238/4‑92 (MQ=37)
                                                                                                                            cttggCGCTGATGCGCCAGATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGAt             >  1:86135/4‑149 (MQ=255)
                                                                                                                            cttggCGCTGATGCGCCAGATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGAt             >  1:66362/4‑149 (MQ=255)
                                                                                                                             ttggCGCTGATGCGCCAGATCAGCGCCTAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATc            >  4:190807/3‑149 (MQ=255)
                                                                                                                             ttggCGCTGATGCGCCAGATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGt                          <  1:227745/133‑1 (MQ=255)
                                                                                                                             ttggCGCTGATGCGCCAGATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGt                          >  2:227745/3‑135 (MQ=255)
                                                                                                                             ttggCGCTGATGCGCCAGATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATc            >  6:70892/3‑149 (MQ=255)
                                                                                                                              aGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATcc           <  4:262037/149‑1 (MQ=255)
                                                                                                                                 cgcTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCCAGa        <  7:76503/149‑1 (MQ=255)
                                                                                                                                    tGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCgg                                                 <  7:29508/105‑1 (MQ=255)
                                                                                                                                    tGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCgg                                                 >  8:29508/1‑105 (MQ=255)
                                                                                                                                    tGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCCAGATCg     >  8:151632/1‑149 (MQ=255)
                                                                                                                                     gATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGTAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCCAGATCGt    >  2:100525/1‑149 (MQ=255)
                                                                                                                                        cTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCCAGATCGTTc  <  3:22736/148‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                   
GCAGCAGCTCCTGCAGGGCCCCACACGCGCTCAGCACCTTGAAGAATACCTGCGGTTGCACTTCCATCAGCGCGCGCTCGATTTCCTTCCAGCTGCGCTCGGCAGTGAGTGCCTGCAGTTCGCCTGAGGCGCTGATCTGGCGCATCAGCGCCAAGGTTTCATCGGCAACCCGGAACCCCAGTGGTGCGTAACGGGCAGCAAAGCGTGCAACGCGCAGCACCCGCAAAGGATCTTCGGCGAATGCCGGGGAAACATGGCGTAGCAGGCGTTGATCCAGATCGTTC  >  NC_002947/478847‑479130

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