Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,334,995 C→A 7.3% E488* (GAA→TAA)  PP_2052 ← bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,334,9950CA 7.3% 200.9 / 4.4 64E488* (GAA→TAA) PP_2052bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base A (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

CCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGCGTGTGCTTGC  >  NC_002947/2334847‑2335139
                                                                                                                                                    |                                                                                                                                                
ccGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTc                                                                                                                                                  >  1:167511/1‑149 (MQ=255)
   aCCCAGGCATGGCCCAGCAGGCTGGCTTACCAGGCCACAATGGCGGGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGc                                                                                                                                               <  5:112595/149‑1 (MQ=255)
    ccagggCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGTGCAGATCGCTGCCGCGTTCGGCg                                                                                                                                              >  8:214120/5‑149 (MQ=255)
     ccAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCg                                                                                                                                              >  5:130083/1‑148 (MQ=255)
         gCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAg                                                                                                                                         >  5:75555/1‑149 (MQ=255)
           aTGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGTGCAGATCGCTGCCGCGTTCGGCGTAAAGCg                                                                                                                                       <  7:214120/149‑1 (MQ=255)
           aTGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCg                                                                                                                                       >  5:127215/1‑149 (MQ=255)
              gCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCTTTACGCCGAACGCGGCa                                                                                                                                    <  1:286702/149‑1 (MQ=255)
                ccAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATg                                                                                                                                  >  2:2668/1‑149 (MQ=255)
                    cagGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCTTTACGCCGAACGCGGCagcgatc                                                                                                                              <  5:100591/149‑7 (MQ=255)
                      ggCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCt                                                                                                                            <  6:130083/149‑1 (MQ=255)
                      ggCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCt                                                                                                                            <  5:147240/149‑1 (MQ=255)
                      ggCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCt                                                                                                                            <  1:2668/149‑1 (MQ=255)
                           gCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCa                                                                                                                         <  2:167511/147‑1 (MQ=255)
                           gCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCAc                                                                                                                        <  6:287224/148‑1 (MQ=255)
                           gCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTGCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCa                                                                                                                         <  3:251667/147‑1 (MQ=255)
                             gTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGt                                                                                                                     <  8:11128/149‑1 (MQ=255)
                              tACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGtt                                                                                                                    >  7:334929/1‑149 (MQ=255)
                               aCCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGc                                                                                                                                     >  1:25271/1‑131 (MQ=255)
                               aCCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGc                                                                                                                                     <  2:25271/131‑1 (MQ=255)
                                ccAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTg                                                                                                                   <  5:30529/148‑1 (MQ=255)
                                  aGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAAc                                                                                                                <  6:75555/149‑1 (MQ=255)
                                    gCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCGGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAg                                                                                                              >  6:368551/1‑149 (MQ=255)
                                         aaTGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCAc                                                                                                           >  5:11063/1‑147 (MQ=255)
                                         aaTGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCa                                                                                                         <  6:127215/149‑1 (MQ=255)
                                                     gCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGc                                                                                             <  7:118618/149‑1 (MQ=255)
                                                           cgacaGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATg                                                                                        <  8:47201/146‑1 (MQ=255)
                                                            cacaGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGgc                                                                                      >  3:108647/1‑149 (MQ=255)
                                                             acaGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGgcg                                                                                     >  3:275986/1‑149 (MQ=255)
                                                                 gCGGTTCTGGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGGATGTCGGCTTCACTGGTGAACAGCACCAGATGCAGCTGGCTCATGGCgctga                                                                                 >  8:162757/1‑148 (MQ=255)
                                                                     ttCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTggg                                                                             >  4:188729/1‑149 (MQ=255)
                                                                         tgatgaTCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCtgt                                                                         >  4:233883/1‑149 (MQ=255)
                                                                          gatgatCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGt                                                                                                                     <  1:32212/104‑1 (MQ=255)
                                                                          gatgatCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGt                                                                                                                     >  2:32212/1‑104 (MQ=255)
                                                                          gatgatCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCtgtg                                                                        >  5:89741/1‑149 (MQ=255)
                                                                             gatctggatctgGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCTTTACGCCGAACGCGGCagc                                                                                                                                  <  3:103559/88‑3 (MQ=37)
                                                                             gatctggatctgGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCTTTACGCCGAACGCGGCagc                                                                                                                                  >  4:103559/1‑86 (MQ=37)
                                                                                    atctggGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGc                                                                <  3:43421/147‑1 (MQ=255)
                                                                                    atctggGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGc                                                                >  4:43421/1‑147 (MQ=255)
                                                                                    atctggGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCt                                                              >  5:128326/1‑149 (MQ=255)
                                                                                      ctggGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTg                                                             >  6:74392/1‑148 (MQ=255)
                                                                                         ggCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGaa                                                           <  3:233883/147‑1 (MQ=255)
                                                                                         ggCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGt                                                         <  6:128326/149‑1 (MQ=255)
                                                                                               gTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGa                                                   >  5:99150/1‑149 (MQ=255)
                                                                                                  tGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGAt                                                  >  4:61672/1‑147 (MQ=255)
                                                                                                       cTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGt                                           >  7:247139/1‑149 (MQ=255)
                                                                                                       cTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGt                                           <  6:89741/149‑1 (MQ=255)
                                                                                                          tCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCg                                        <  6:99150/149‑1 (MQ=255)
                                                                                                                 cTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGATCGTACTCCg                                 >  7:21900/1‑149 (MQ=255)
                                                                                                                         ggCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAgcggt                         <  3:188729/149‑2 (MQ=255)
                                                                                                                          gCGTTCGAGCAGATCGCTGCCGCGTTAGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAAAAGCACCAGATGCATCTGGCTCAGGGCGCTGCTGGGCTGTGGAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTACGGTAgctgcg                        >  4:52648/1‑149 (MQ=255)
                                                                                                                                gAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGaa                   <  6:11063/148‑1 (MQ=255)
                                                                                                                                gAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAg                  <  4:275986/149‑1 (MQ=255)
                                                                                                                                 aGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTg                     >  7:174343/1‑145 (MQ=255)
                                                                                                                                   cAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGAtt               >  5:351865/1‑149 (MQ=255)
                                                                                                                                      aTCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGc            >  4:327437/1‑149 (MQ=255)
                                                                                                                                      aTCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCAGCGTTGAAGATTTGc            <  8:247139/149‑1 (MQ=255)
                                                                                                                                       tCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGCg           <  5:74392/149‑1 (MQ=255)
                                                                                                                                       tCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGCg           >  8:185037/1‑149 (MQ=255)
                                                                                                                                         gCTGCCGCGTTCGGCGTAAAGCTGCATGTCGGCTTCACTGTTGAACAGCACCTGATGCAGCTGGCGCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGCgt          <  8:334929/148‑1 (MQ=255)
                                                                                                                                         gCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCCGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGCgt          <  5:368551/148‑1 (MQ=255)
                                                                                                                                           tGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGCGtgtg       >  3:320680/1‑149 (MQ=255)
                                                                                                                                           tGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGCGtgtg       >  5:91824/1‑149 (MQ=255)
                                                                                                                                             ccccgTACGGCGTAAACCCGCATGACGGCTTCACTGTTGAACAACACCAGATGCAGCTGGCTAATGGCCCTGCTCGGCTGTCCAAAAGGCCTGAACACGCCGATCAGCCGTTCGTACTCCGGTAGCCGCGTAGAAGATATGCGAGTGCt     >  4:138961/4‑149 (MQ=255)
                                                                                                                                                cgTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGCGTGTGCTTGc  <  3:327437/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCGCGTTCGGCGTAAAGCGGCATGTCGGCTTCACTGTTGAACAGCACCAGATGCAGCTGGCTCATGGCGCTGCTGGGCTGTGCAAAAGGCCTGAAGTGGCCGATCAGCCGTTCGTACTCCGGTAGCGGCGTTGAAGATTTGCGTGTGCTTGC  >  NC_002947/2334847‑2335139

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: ‑