Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,806,760 C→G 6.7% P459R (CCC→CGC)  PP_3364 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,806,7600CG 6.7% 213.7 / 3.9 71P459R (CCC→CGC) PP_3364two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (35/30);  new base G (3/2);  total (39/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 = 8.37e-01

GTGTGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCAGTGAACGC  >  NC_002947/3806622‑3806905
                                                                                                                                          |                                                                                                                                                 
gtgtGCTGATTATTCCGATGTCGCAACACCCTAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGTCAGCGGGGCGTCAGGCGGTCGCCCAACGcgccgctgtc                                                                                                                                           <  7:55057/147‑5 (MQ=37)
gtgtGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAaga                                                                                                                                           <  2:32171/147‑1 (MQ=255)
gtgtGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGc                                                                                                                                         >  7:52088/1‑149 (MQ=255)
   tgCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTc                                                                                                                                      <  7:3776/149‑1 (MQ=255)
           attCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTg                                                                                                                              >  7:188776/1‑149 (MQ=255)
               cGATGTCGCAACACCCAAGTGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGgaaa                                                                                                                          <  4:220087/149‑1 (MQ=255)
               cGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGgaa                                                                                                                           <  2:354751/148‑1 (MQ=255)
               cGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGga                                                                                                                            >  3:288561/1‑147 (MQ=255)
                gATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGga                                                                                                                            >  8:248242/1‑146 (MQ=255)
                gATGTCGAAACACCCAACGGATTACTAGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCAGGGGACCCGAACAAGCCCTATGACAGCGGCGGGTTGGGCGACCGCCAGACGCCCCGCAAGAGCTTCCCCATCTgaa                                                                                                                            >  4:214752/1‑144 (MQ=255)
                    tCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAAc                                                                                                                     <  4:288561/149‑1 (MQ=255)
                      gCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCg                                                                                                                   >  7:3558/1‑149 (MQ=255)
                              cAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACcag                                                                                                           >  7:33205/1‑149 (MQ=255)
                              cAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACcag                                                                                                           >  7:134926/1‑149 (MQ=255)
                               aaGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACcag                                                                                                           >  3:314310/1‑148 (MQ=255)
                               aaGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACcag                                                                                                           >  3:87833/1‑148 (MQ=255)
                                  ggATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGt                                                                                                       <  6:71309/149‑1 (MQ=255)
                                    aTTACTTTATTTTCTTCCGCTTGGAGGGGGTCGAGACACTCGATTGGGCCGGGGACCCTAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCg                                                                                                                   <  4:227769/135‑1 (MQ=255)
                                    aTTACTTGATCTTCTTCCGAAAGGAGGTGGTCGAGACACTCGATTGGGCAGGGTACACGAACAATACCTATGACAGCGGCGCGTTGGGATACCGCCTGACGACCCGCAAGAGCTTCGCAAACTGGAAAGAAACCg                                                                                                                   >  3:227769/1‑135 (MQ=255)
                                         ttGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTAcc                                                                                                >  5:10018/1‑149 (MQ=255)
                                         ttGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTAcc                                                                                                >  6:178028/1‑149 (MQ=255)
                                           gATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCt                                                                                              >  4:139770/1‑149 (MQ=255)
                                            aTCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCt                                                                                              <  7:248242/148‑1 (MQ=255)
                                              cttcttCCGCAAGGAGGTGGTCGTGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGa                                                                                           >  4:129216/1‑149 (MQ=255)
                                                tcttcCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGAcc                                                                                         >  3:278745/1‑149 (MQ=255)
                                                tcttcCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGAcc                                                                                         >  5:267320/1‑149 (MQ=255)
                                                   tcCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGa                                                                                       <  6:10018/148‑1 (MQ=255)
                                                    ccGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAg                                                                                      >  3:336189/1‑148 (MQ=255)
                                                    ccGAAATGAGGTGGTCGAGACACTCGATTGGGCCGGGTACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAATAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAg                                                                                      <  8:175944/148‑1 (MQ=255)
                                                      gCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACcag                                                                                                           <  6:211156/125‑1 (MQ=255)
                                                      gCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACcag                                                                                                           >  5:211156/1‑125 (MQ=255)
                                                      gCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAg                                                                                   >  7:99014/1‑149 (MQ=255)
                                                        aaGGAGGTGGTCGAGACACTCGATTTGGCCGGTGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGa                                                                                 <  7:97612/149‑1 (MQ=255)
                                                        aaGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCACAGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTAACCAGCAGTCGCTACCCTGGACCGAGCAGGa                                                                                 >  5:211002/1‑149 (MQ=255)
                                                            aggTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGc                                                                             <  6:267320/149‑1 (MQ=255)
                                                             ggtggtCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGc                                                                             >  6:264002/1‑148 (MQ=255)
                                                              gtggtCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCa                                                                           <  5:88624/149‑1 (MQ=255)
                                                                      gaCACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGc                                                                    <  8:3558/148‑1 (MQ=255)
                                                                         acTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATg                                                                >  7:34400/1‑149 (MQ=255)
                                                                         aatCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACAAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATg                                                                >  3:219015/3‑149 (MQ=255)
                                                                          cTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGc                                                               >  8:157999/1‑149 (MQ=255)
                                                                           tCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGc                                                               >  3:74783/1‑148 (MQ=255)
                                                                           tCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCATGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGc                                                               >  6:233468/1‑148 (MQ=255)
                                                                                      ggggACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGc                                                   <  5:233468/149‑1 (MQ=255)
                                                                                          aCCCGAACAAGACCTATGACAGCGGCGCGGTGGGCTACCTCCTGACGCCCCGAAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACAGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTAtc                                               <  6:211002/149‑1 (MQ=255)
                                                                                             cGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTAtcgtc                                            <  8:52088/149‑1 (MQ=255)
                                                                                              gAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTAtcgtc                                            >  8:64310/1‑148 (MQ=255)
                                                                                               aaCAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGa                                          <  5:264002/149‑1 (MQ=255)
                                                                                               aaCAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGa                                          <  7:64310/149‑1 (MQ=255)
                                                                                                aCAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGa                                          <  8:99014/148‑1 (MQ=255)
                                                                                                aCAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAg                                         >  3:370311/1‑149 (MQ=255)
                                                                                                aCAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAg                                         <  8:33205/149‑1 (MQ=255)
                                                                                                   aGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAggtg                                      <  3:214752/149‑1 (MQ=255)
                                                                                                     aCCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAggtg                                      >  7:114831/1‑147 (MQ=255)
                                                                                                     aCCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAggtg                                      <  8:114831/147‑1 (MQ=255)
                                                                                                       cTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAggtgg                                     <  6:136832/146‑1 (MQ=255)
                                                                                                       cTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAggtgg                                     >  5:136832/1‑146 (MQ=255)
                                                                                                        tATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGtgct                                 >  7:46288/1‑149 (MQ=255)
                                                                                                           tcttgcggggCGTCAGGCGGTCGCCCAACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGca                              >  5:255740/5‑149 (MQ=38)
                                                                                                           tcttgcggggCGTCAGGCGGTCGCCCAACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGca                              >  2:265365/5‑149 (MQ=38)
                                                                                                             ttgcggggCGTCAGGCGGTCGCCCAACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGcac                             >  5:115467/3‑148 (MQ=255)
                                                                                                                cggcgCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATa                         <  4:370311/149‑1 (MQ=255)
                                                                                                                          ggCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGc                <  5:114065/148‑1 (MQ=255)
                                                                                                                            cGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCg             <  7:157999/149‑1 (MQ=255)
                                                                                                                             gACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGACATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGAGAACTGCTctc            >  6:19124/1‑147 (MQ=255)
                                                                                                                                  ccTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCAGt        >  3:197572/1‑148 (MQ=255)
                                                                                                                                   ccaaCGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCAGTGa      <  6:255740/146‑1 (MQ=255)
                                                                                                                                    tGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCAGTGaa     >  8:5611/1‑149 (MQ=255)
                                                                                                                                       cGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCAGTGAAcgc  <  7:5611/149‑1 (MQ=255)
                                                                                                                                       cGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCAGTGAAcgc  <  4:219015/149‑1 (MQ=255)
                                                                                                                                        gCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATg                                                                >  1:227276/1‑86 (MQ=255)
                                                                                                                                        gCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATg                                                                <  2:227276/86‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
GTGTGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCAGTGAACGC  >  NC_002947/3806622‑3806905

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: