Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,525,990 C→T 6.1% S307S (AGC→AGT 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,9900CT 6.1% 248.1 / 3.8 81S307S (AGC→AGTmnmAtRNA 2‑thiouridine(34) synthase MnmA
Reads supporting (aligned to +/- strand):  ref base C (37/39);  new base T (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 = 9.85e-01

ATGCCAGTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAG  >  NC_002947/4525842‑4526138
                                                                                                                                                     |                                                                                                                                                    
aTGCCAGTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGc                                                                                                                                                      <  1:357839/149‑1 (MQ=255)
   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                                                                                                  >  6:289932/1‑121 (MQ=255)
                                                                                ccATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCa                                                                                                  <  5:289932/121‑1 (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                                                                           <  1:180158/139‑1 (MQ=255)
                                                                                     gCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGc                                                                           >  2:180158/1‑139 (MQ=255)
                                                                                      cTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCgt                                                                >  8:25435/1‑149 (MQ=255)
                                                                                       tGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCgtg                                                               >  6:52548/1‑149 (MQ=255)
                                                                                           ctcCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCAc                                                                                     <  5:59134/123‑1 (MQ=255)
                                                                                           ctcCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCAc                                                                                     >  6:59134/1‑123 (MQ=255)
                                                                                              ccGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTcgac                                                         <  5:52548/148‑1 (MQ=255)
                                                                                              ccGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTcga                                                          <  2:218872/147‑1 (MQ=255)
                                                                                                 cgcCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGc                                                     <  3:57875/149‑1 (MQ=255)
                                                                                                   cccTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTcgac                                                         <  2:321353/143‑1 (MQ=255)
                                                                                                   cccTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTcgac                                                         >  1:321353/1‑143 (MQ=255)
                                                                                                   cccTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCt                                                   >  5:358723/1‑149 (MQ=255)
                                                                                                     ctgctgGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCa                                                 >  6:179364/1‑149 (MQ=255)
                                                                                                              cTCGGAAATTTTATGGGTCAACCCGGTCGACCTCGGCACCCCACGGAAGTTTACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGAAGCACCCTGGAACGCCCCGCCAGCGGCTACAGCGCCGTGTTCTACGAGCCTCACCGAGCCg                                         >  2:326256/1‑149 (MQ=255)
                                                                                                                    aaTTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGc                                   <  5:357875/149‑1 (MQ=255)
                                                                                                                    aaTTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGc                                   >  3:197562/1‑149 (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                                                                                          >  5:109342/1‑78 (MQ=255)
                                                                                                                                    ccGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGaa                                                                                          <  6:109342/78‑1 (MQ=255)
                                                                                                                                     cGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCGGGGCCAGTCGgtggtg                  <  5:179364/149‑1 (MQ=255)
                                                                                                                                     cGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGgtggtg                  <  4:197562/149‑1 (MQ=255)
                                                                                                                                     cGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGgtggtg                  <  1:99727/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                     >  3:371832/1‑133 (MQ=255)
                                                                                                                                                  cAGCCCACGGCAGTTGACTGCCAAGGTGTGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCggtg                     <  4:371832/133‑1 (MQ=255)
                                                                                                                                                  cAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGc      >  7:106881/1‑148 (MQ=255)
                                                                                                                                                    gCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGaa   <  3:43464/149‑1 (MQ=255)
                                                                                                                                                     cccACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAg  >  1:160946/1‑149 (MQ=255)
                                                                                                                                                     cccACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAg  >  6:151916/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGCCAGTGACGAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAG  >  NC_002947/4525842‑4526138

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: