Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 638,742 A→T 6.8% L16Q (CTG→CAG)  PP_0550 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947638,7420AT 6.8% 217.0 / 4.9 73L16Q (CTG→CAG) PP_0550hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/34);  new base T (3/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GCCCGGCCAGTTGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGCTTCTTCGGC  >  NC_002947/638596‑638887
                                                                                                                                                  |                                                                                                                                                 
gcccgGCCAGTTGCGGAAGGAATACCGCCGGGCCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGa                                                                                                                                                 >  8:104302/1‑149 (MQ=255)
    ggCCAGTTGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCgg                                                                                                                                             >  6:83090/1‑149 (MQ=255)
     gCCAGTTGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGt                                                                                                                                            <  1:113071/149‑1 (MQ=255)
     gCCAGTTGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGt                                                                                                                                            <  5:84689/149‑1 (MQ=255)
      ccAGTTGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTg                                                                                                                                           <  3:359551/149‑1 (MQ=255)
          ttGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTc                                                                                                                                       <  8:14386/149‑1 (MQ=255)
             cGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCTGACCGGTgaacacc                                                                                                                                     <  6:194613/148‑7 (MQ=255)
             cGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCgg                                                                                                                                    >  2:224177/1‑149 (MQ=255)
              ggaaggaaTACCGCCGGGTCGGCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGc                                                                                                                                   <  2:57239/149‑1 (MQ=255)
                aaggaaTACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCaa                                                                                                                                 <  1:60636/149‑1 (MQ=255)
                 aggaaTACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGGGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAg                                                                                                                                >  4:198817/1‑149 (MQ=255)
                    aaTACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAg                                                                                                                             <  1:334247/149‑1 (MQ=255)
                     aTACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCTGACCGGTgaacacc                                                                                                                                     >  5:385566/1‑134 (MQ=255)
                     aTACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCTGACCGGTgaacacc                                                                                                                                     <  6:385566/140‑7 (MQ=255)
                      tACCGCCGGGTCACCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGa                                                                                                                           <  3:89838/149‑1 (MQ=255)
                          gccgGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATg                                                                                                                       <  5:171775/149‑1 (MQ=255)
                                      aCCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATg                                                                                                                       >  4:12434/1‑137 (MQ=255)
                                      aCCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATg                                                                                                                       <  3:12434/137‑1 (MQ=255)
                                        cGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGGGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGt                                                                                                         <  2:211580/149‑1 (MQ=255)
                                            ttttCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACat                                                                                                     <  7:104302/149‑1 (MQ=255)
                                                ccAGGTGCGGGTGCTCGGCCAATGGCCGGCGTTTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGATCAGGCTGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACatgat                                                                                                  <  7:121151/148‑1 (MQ=255)
                                                 cAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCACGACATGCTGAGGCGTCAAGTACatgatg                                                                                                 <  1:217326/148‑1 (MQ=255)
                                                        gggTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTccc                                                                                         <  2:71243/149‑1 (MQ=255)
                                                        gggTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTccc                                                                                         >  6:186590/1‑149 (MQ=255)
                                                           tGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTg                                                                                      <  7:36280/149‑1 (MQ=255)
                                                              tCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCGCCGGCAAGCAGGACCTGCTGAGGCGTCAAGTCCATGATGTCTCTCCCGTGTTa                                                                                   >  6:43920/1‑149 (MQ=255)
                                                                ggCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTa                                                                                   <  5:83090/147‑1 (MQ=255)
                                                                    aaTGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTc                                                                              >  5:89467/1‑148 (MQ=255)
                                                                      tGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAg                                                                           <  3:184472/149‑1 (MQ=255)
                                                                        gCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTc                                                                         <  1:363257/149‑1 (MQ=255)
                                                                          cGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGATCAGGCGGTCTAGCAGGTGTTCACCTGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTa                                                                       <  1:119415/149‑1 (MQ=255)
                                                                            gCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTagt                                                                     <  5:101151/149‑1 (MQ=255)
                                                                                    gcgcCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCg                                                             >  3:212871/1‑149 (MQ=255)
                                                                                      gcCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCggg                                                           >  3:235460/1‑149 (MQ=255)
                                                                                      gcCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCggg                                                           >  3:344463/1‑149 (MQ=255)
                                                                                        cAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTAAGACCTGTGGCTCCGGCAAGCAGGACATGCGGAGGCGTAAAGTACATGATGTGTCTCCCTTGTTACACTCAAGTCTAGTGACTGCCGGGCa                                                         >  5:61510/1‑149 (MQ=255)
                                                                                         aGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATg                                                                                                                       >  3:373433/1‑86 (MQ=255)
                                                                                         aGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATg                                                                                                                       <  4:373433/86‑1 (MQ=255)
                                                                                              acGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTc                                                                              <  6:52807/122‑1 (MQ=255)
                                                                                              acGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTc                                                                              >  5:52807/1‑122 (MQ=255)
                                                                                              acGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTg                                                   >  6:276905/1‑149 (MQ=255)
                                                                                                gggTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTagt                                                                     >  1:136631/1‑129 (MQ=255)
                                                                                                gggTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTagt                                                                     <  2:136631/129‑1 (MQ=255)
                                                                                                  gTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCtg                                               >  4:380990/1‑149 (MQ=255)
                                                                                                   tGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCtgc                                              <  3:190952/149‑1 (MQ=255)
                                                                                                   tGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCtgc                                              <  1:224177/149‑1 (MQ=255)
                                                                                                           ggCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCAc                                      >  1:368684/1‑149 (MQ=255)
                                                                                                            gCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGACAAGTACAAGATGTAGCTCCCTTGTTACACTCAAGACTAGTGAGTGCCGCGCATTTTTGCCTGCTTTCTCACg                                     >  5:383586/1‑149 (MQ=255)
                                                                                                                gTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTgg                                 >  4:4336/1‑149 (MQ=255)
                                                                                                                gTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTg                                  <  3:198817/148‑1 (MQ=255)
                                                                                                                     gCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCtgt                            <  6:89467/149‑1 (MQ=255)
                                                                                                                       aGGCGGTCTAGCAGGTGTTCACCGGTCAGTCCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCtgtgt                          <  6:61510/149‑1 (MQ=255)
                                                                                                                            gTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGt                     >  4:366654/1‑149 (MQ=255)
                                                                                                                             tCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGt                     >  5:51589/1‑148 (MQ=255)
                                                                                                                               tAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGGCATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCAGGCATGTTTGCCTGCTTGCTGACGTTGGCCTGTGTGGCGTTAc                  >  6:378526/1‑149 (MQ=255)
                                                                                                                                aGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACg                 >  5:85165/1‑149 (MQ=255)
                                                                                                                                 gCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTAGATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACg                 >  3:281255/1‑148 (MQ=255)
                                                                                                                                   aGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGc                                                  >  3:240954/1‑113 (MQ=255)
                                                                                                                                   aGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGc                                                  <  4:240954/113‑1 (MQ=255)
                                                                                                                                    ggTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAAc             >  2:40578/1‑149 (MQ=255)
                                                                                                                                    ggTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGGGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAAc             >  1:320865/1‑149 (MQ=255)
                                                                                                                                     gtgtTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCt                                                <  6:294860/113‑1 (MQ=255)
                                                                                                                                     gtgtTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCt                                                >  5:294860/1‑113 (MQ=255)
                                                                                                                                      gagccACCGGTCTGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGc           >  6:276273/5‑149 (MQ=255)
                                                                                                                                      gagccACCGGTCTGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGc           >  7:146671/5‑149 (MQ=255)
                                                                                                                                        ttCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGct          >  4:170983/1‑148 (MQ=255)
                                                                                                                                           aCCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGcttctt      >  5:187661/1‑149 (MQ=255)
                                                                                                                                           aCCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGcttctt      >  4:237471/1‑149 (MQ=255)
                                                                                                                                           aCCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGcttct       <  2:368684/148‑1 (MQ=255)
                                                                                                                                            ccGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGCttcttc     >  7:258450/1‑149 (MQ=255)
                                                                                                                                            ccGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGCttcttc     <  4:344463/149‑1 (MQ=255)
                                                                                                                                               gTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGCTTCTTCGGc  <  3:366654/149‑1 (MQ=255)
                                                                                                                                                tCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGc                                                  >  6:205001/1‑100 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCCCGGCCAGTTGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGCTTCTTCGGC  >  NC_002947/638596‑638887

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑