Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,525,998 A→C 6.0% Q310P (CAG→CCG)  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,9980AC 6.0% 257.2 / 4.2 80Q310P (CAG→CCG) mnmAtRNA 2‑thiouridine(34) synthase MnmA
Reads supporting (aligned to +/- strand):  ref base A (36/39);  new base C (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 = 1.00e+00

GAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAGTGTGCCTG  >  NC_002947/4525853‑4526146
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                                        <  1:199851/104‑1 (MQ=37)
                                                        gggCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGTCAACTGCCGTGGGCTGCCGAgg                                                                                                                                        >  2:199851/1‑104 (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                                                                                                       <  4:311708/122‑1 (MQ=255)
                                                                       aTGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCAccc                                                                                                       >  3:311708/1‑122 (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                                                                                             >  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                                                                 <  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                                           >  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                              >  1:339068/1‑149 (MQ=255)
                                                                                                                      aaCCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCggt                              <  2:339068/149‑1 (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          >  1:160946/1‑149 (MQ=255)
                                                                                                                                          cccACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAg          >  6:151916/1‑149 (MQ=255)
                                                                                                                                           ccACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAgt         <  6:358723/149‑1 (MQ=255)
                                                                                                                                                  aGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAGTGTGCCTg  >  2:75051/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAGCCGTGGTACGTGCTGCACAAGGACCTGGCCCGCAACGTGCTGGTGGTCGGCCAGGGCAACGAGCACCCATGGCTGTTCTCCCGCGCCCTGCTGGCCTCGGAAATTTTCTGGGTCAACCCGGTCGACCTCGGCAGCCCACGGCAGTTGACTGCCAAGGTGCGCTACCGCCAGAGCGACCAGCGATGCACCCTGGAACGCACCGCCAGCGGCTACCGCGCCGTGTTCGACGAGCCTCAGCGCGCCGTGACGCCGGGCCAGTCGGTGGTGTTCTACGACGGCGAAGTGTGCCTG  >  NC_002947/4525853‑4526146

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: