Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,525,992 C→A 6.2% P308Q (CCA→CAA)  mnmA → tRNA 2‑thiouridine(34) synthase MnmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,525,9920CA 6.2% 263.9 / 5.9 81P308Q (CCA→CAA) mnmAtRNA 2‑thiouridine(34) synthase MnmA
Reads supporting (aligned to +/- strand):  ref base C (37/39);  new base A (3/2);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

CCAGTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAGT  >  NC_002947/4525845‑4526139
                                                                                                                                                    |                                                                                                                                                   
ccAGTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCa                                                                                                                                                    >  2:362861/1‑149 (MQ=255)
  aGTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACg                                                                                                                                                  >  3:371008/1‑149 (MQ=255)
   gTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACgg                                                                                                                                                 <  5:133067/149‑1 (MQ=255)
    tGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACgg                                                                                                                                                 >  8:203235/1‑148 (MQ=255)
        gtgCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTCGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGtt                                                                                                                                            <  1:98619/147‑1 (MQ=255)
        gAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGtt                                                                                                                                            <  8:335243/149‑1 (MQ=255)
         aGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTg                                                                                                                                           <  7:203235/149‑1 (MQ=255)
          gCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGa                                                                                                                                          <  1:47407/149‑1 (MQ=255)
             gTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTg                                                                                                                                       >  7:302047/1‑149 (MQ=255)
                gTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGcc                                                                                                                                     >  5:146820/1‑148 (MQ=255)
                   cGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGTCAACTGCCGTGGGCTGCCGAgg                                                                                                                                 <  1:272458/149‑1 (MQ=255)
                   cGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGGCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAg                                                                                                                                  <  4:215244/148‑1 (MQ=255)
                     tgctgcACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTg                                                                                                                               >  2:310008/1‑149 (MQ=255)
                     tgcagcACAAGGACCTGGCACGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCCCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTg                                                                                                                               >  2:219641/1‑149 (MQ=255)
                         gcACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTTGGCAGCCCACGGCAGTTGACTGCCAAGGTgcgc                                                                                                                            >  5:338946/1‑148 (MQ=255)
                          cacaAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTgcgc                                                                                                                            <  1:362861/147‑1 (MQ=255)
                           acaAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTAc                                                                                                                         <  6:146820/149‑1 (MQ=255)
                            caAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTTGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTAcc                                                                                                                        <  6:338946/149‑1 (MQ=255)
                                      gCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGc                                                                                                                <  6:117750/147‑1 (MQ=255)
                                       cccGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGAc                                                                                                             <  8:223014/149‑1 (MQ=255)
                                           cAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGc                                                                                                         <  6:151431/149‑1 (MQ=255)
                                            aaCGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCg                                                                                                        <  1:310008/149‑1 (MQ=255)
                                                        gtCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTgg                                                                                             >  6:166293/1‑148 (MQ=255)
                                                            gCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACg                                                                                         >  4:318428/1‑148 (MQ=255)
                                                              cAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCAc                                                                                      <  8:302047/149‑1 (MQ=255)
                                                                gggCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGTCAACTGCCGTGGGCTGCCGAgg                                                                                                                                 >  2:199851/1‑104 (MQ=37)
                                                                gggCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGTCAACTGCCGTGGGCTGCCGAgg                                                                                                                                 <  1:199851/104‑1 (MQ=37)
                                                                   cAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGc                                                                                   >  3:182696/1‑147 (MQ=255)
                                                                       gAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCgg                                                                             <  2:22697/149‑1 (MQ=255)
                                                                          cACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTa                                                                          >  3:10938/1‑149 (MQ=255)
                                                                            cccATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTAcc                                                                        >  1:218872/1‑149 (MQ=255)
                                                                             ccATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCa                                                                                                   <  5:289932/121‑1 (MQ=255)
                                                                             ccATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCa                                                                                                   >  6:289932/1‑121 (MQ=255)
                                                                             ccATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACcg                                                                       >  2:99727/1‑149 (MQ=255)
                                                                               aTGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCAccc                                                                                                >  3:311708/1‑122 (MQ=255)
                                                                               aTGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCAccc                                                                                                <  4:311708/122‑1 (MQ=255)
                                                                                tGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCgcgc                                                                    >  6:357875/1‑149 (MQ=255)
                                                                                 ggCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGcc                                                                   >  6:328706/1‑149 (MQ=255)
                                                                                  gCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGc                                                                            >  2:180158/1‑139 (MQ=255)
                                                                                  gCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGc                                                                            <  1:180158/139‑1 (MQ=255)
                                                                                   cTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCgt                                                                 >  8:25435/1‑149 (MQ=255)
                                                                                    tGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCgtg                                                                >  6:52548/1‑149 (MQ=255)
                                                                                        ctcCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCAc                                                                                      >  6:59134/1‑123 (MQ=255)
                                                                                        ctcCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCAc                                                                                      <  5:59134/123‑1 (MQ=255)
                                                                                           ccGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTcgac                                                          <  5:52548/148‑1 (MQ=255)
                                                                                           ccGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTcga                                                           <  2:218872/147‑1 (MQ=255)
                                                                                              cgcCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGc                                                      <  3:57875/149‑1 (MQ=255)
                                                                                                cccTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTcgac                                                          >  1:321353/1‑143 (MQ=255)
                                                                                                cccTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTcgac                                                          <  2:321353/143‑1 (MQ=255)
                                                                                                cccTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCt                                                    >  5:358723/1‑149 (MQ=255)
                                                                                                  ctgctgGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCa                                                  >  6:179364/1‑149 (MQ=255)
                                                                                                           cTCGGAAATTTTATGGGTCAACCCGGTCGACCTCGGCACCCCACGGAAGTTTACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGAAGCACCCTGGAACGCCCCGCCAGCGGCTACAGCGCCGTGTTCTACGAGCCTCACCGAGCCg                                          >  2:326256/1‑149 (MQ=255)
                                                                                                                 aaTTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGc                                    >  3:197562/1‑149 (MQ=255)
                                                                                                                 aaTTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGc                                    <  5:357875/149‑1 (MQ=255)
                                                                                                                   ttttCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCg                                  <  7:25435/149‑1 (MQ=255)
                                                                                                                     ttCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCgg                                 <  4:371008/148‑1 (MQ=255)
                                                                                                                             cAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCgg                        >  1:36778/1‑149 (MQ=255)
                                                                                                                              aaCCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCggt                       <  2:339068/149‑1 (MQ=255)
                                                                                                                              aaCCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCggt                       >  1:339068/1‑149 (MQ=255)
                                                                                                                               aCCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCGAGACCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCCCGCCGTGTTCGCCGAGCCTCAGCGCGCCGTGCTAGCGGGCCAGTCggtg                      >  1:320733/1‑149 (MQ=255)
                                                                                                                                cccGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCggtgg                     >  7:362960/1‑149 (MQ=255)
                                                                                                                                 ccGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGaa                                                                                           <  6:109342/78‑1 (MQ=255)
                                                                                                                                 ccGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGaa                                                                                           >  5:109342/1‑78 (MQ=255)
                                                                                                                                  cGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCGGGGCCAGTCGgtggtg                   <  5:179364/149‑1 (MQ=255)
                                                                                                                                  cGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGgtggtg                   <  1:99727/149‑1 (MQ=255)
                                                                                                                                  cGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGgtggtg                   <  4:197562/149‑1 (MQ=255)
                                                                                                                                    gTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACGGCGCCGTGTGCGACTAGCCTCACCGCGCCGTGACGCCGGGCCAGTCGGTGGTGtt                 <  1:219641/149‑1 (MQ=255)
                                                                                                                                    gTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGtt                 <  8:362960/149‑1 (MQ=255)
                                                                                                                                      gcaCCTTGGCAGTCAACTGCCGTGGGCTGCCGAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTc                >  4:238261/3‑148 (MQ=38)
                                                                                                                                        aCCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTac             >  7:245560/1‑149 (MQ=255)
                                                                                                                                          cttggCAGTCAACTGCCGTGGGCTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTacga           >  8:323801/4‑149 (MQ=255)
                                                                                                                                          ctaggCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCCCCCTGGACCGCACCGCCAGCGGCTACCGCGCCGTGGTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTacgc           >  1:374554/4‑148 (MQ=255)
                                                                                                                                          cTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTacga           <  4:10938/149‑1 (MQ=255)
                                                                                                                                            cGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTacgacg         >  4:43464/1‑149 (MQ=255)
                                                                                                                                             ggCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTacgac          <  5:166293/147‑1 (MQ=255)
                                                                                                                                             ggCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACgg        <  6:9383/149‑1 (MQ=255)
                                                                                                                                              gCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGc       <  5:328706/149‑1 (MQ=255)
                                                                                                                                               cAGCCCACGGCAGTTGACTGCCAAGGTGTGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCggtg                      <  4:371832/133‑1 (MQ=255)
                                                                                                                                               cAGCCCACGGCAGTTGACTGCCAAGGTGTGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCggtg                      >  3:371832/1‑133 (MQ=255)
                                                                                                                                               cAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGc       >  7:106881/1‑148 (MQ=255)
                                                                                                                                                 gCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGaa    <  3:43464/149‑1 (MQ=255)
                                                                                                                                                  cccACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAg   >  6:151916/1‑149 (MQ=255)
                                                                                                                                                  cccACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAg   >  1:160946/1‑149 (MQ=255)
                                                                                                                                                   ccACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAgt  <  6:358723/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAGTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAGT  >  NC_002947/4525845‑4526139

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: