Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,525,995 G→T 6.1% R309L (CGG→CTG)  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,9950GT 6.1% 257.3 / 4.0 80R309L (CGG→CTG) mnmAtRNA 2‑thiouridine(34) synthase MnmA
Reads supporting (aligned to +/- strand):  ref base G (36/39);  new base T (3/2);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

AGTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAGT  >  NC_002947/4525847‑4526139
                                                                                                                                                     |                                                                                                                                                
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)
                                                                                                                                                     |                                                                                                                                                
AGTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAGT  >  NC_002947/4525847‑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: