Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,346,123 C→A 7.2% R126L (CGC→CTC)  PP_2942 ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,346,1230CA 7.2% 204.2 / 4.4 64R126L (CGC→CTC) PP_2942two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (27/32);  new base A (3/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTG  >  NC_002947/3345977‑3346271
                                                                                                                                                    |                                                                                                                                                    
cAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGgcg                                                                                                                                                     >  7:98931/1‑148 (MQ=255)
   cGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAg                                                                                                                                                 <  3:51850/149‑1 (MQ=255)
     gTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGt                                                                                                                                               <  8:98931/149‑1 (MQ=255)
      tCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTc                                                                                                                                              <  2:239489/149‑1 (MQ=255)
          ctgcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAtt                                                                                                                                           <  8:244451/148‑1 (MQ=255)
           tgcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAg                                                                                                                                         >  4:226046/1‑149 (MQ=255)
            gcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAt                                                                                                                                            <  1:324495/145‑1 (MQ=255)
            gcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAt                                                                                                                                            >  2:324495/1‑145 (MQ=255)
             cAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCTTTAGc                                                                                                                                        <  3:101993/148‑1 (MQ=255)
                cTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTACAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAAtt                                                                                                                                    >  5:130537/1‑149 (MQ=255)
                 tCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATtc                                                                                                                                   >  3:36397/1‑149 (MQ=255)
                    gcgGTACGGGCTGCAACGCGTTTCTCGAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATtcgtc                                                                                                                                <  7:262494/149‑1 (MQ=255)
                       gTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATtcgtcgtc                                                                                                                             <  3:222542/149‑1 (MQ=255)
                           gggCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTc                                                                                                                         >  4:275154/1‑149 (MQ=255)
                           gggATGCAACGCGTTTCTCAAAGTCGTTATTTAGCAGTCTCAACTGGGCTTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCTTAAAGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTc                                                                                                                         <  8:204810/149‑1 (MQ=255)
                                          tctcCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGAt                                                                                                           <  8:217948/148‑1 (MQ=255)
                                            tcCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGt                                                                                                         >  5:193947/1‑148 (MQ=255)
                                              cAAGTCGTTATTGAGCAGTCTCAAATGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAACGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAAc                                                                                                      >  7:306675/1‑149 (MQ=255)
                                                 gTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGt                                                                                                   <  6:193947/149‑1 (MQ=255)
                                                  tCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTg                                                                                                  <  2:243091/149‑1 (MQ=255)
                                                    gTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTgg                                                                                                 <  4:288839/148‑1 (MQ=255)
                                                    gTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGa                                                                                                >  5:75329/1‑149 (MQ=255)
                                                    gTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGa                                                                                                >  6:74023/1‑149 (MQ=255)
                                                     ttATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAGTCGGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTgg                                                                                                 <  1:267577/149‑1 (MQ=255)
                                                       aTTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCGGGGCTTACTGATGTAACGGTGGATCt                                                                                             >  4:56923/1‑149 (MQ=255)
                                                           aGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAt                                                                                                                                            <  3:291587/98‑1 (MQ=255)
                                                           aGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAt                                                                                                                                            >  4:291587/1‑98 (MQ=255)
                                                                tctcAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCAt                                                                                    >  1:249315/1‑149 (MQ=255)
                                                                 ctcAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCAtt                                                                                   <  5:254151/149‑1 (MQ=255)
                                                                    aaCTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGa                                                                                 <  3:275154/148‑1 (MQ=255)
                                                                        gggCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCa                                                                                                                       <  6:142432/106‑1 (MQ=255)
                                                                        gggCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCa                                                                                                                       >  5:142432/1‑106 (MQ=255)
                                                                         ggCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGcc                                                                            <  3:128321/148‑1 (MQ=255)
                                                                            gTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGCCGGCCTtga                                                                        >  8:56770/1‑149 (MQ=255)
                                                                            gTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTtga                                                                        <  4:36397/149‑1 (MQ=255)
                                                                            gTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTtga                                                                        <  3:226046/149‑1 (MQ=255)
                                                                             ttCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGCTCTGACCGTCATTGACGGCCTtgat                                                                       >  5:330810/1‑149 (MQ=255)
                                                                                          ccAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTgg                                                           <  6:75329/148‑1 (MQ=255)
                                                                                              cagcTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCgg                                                      >  7:140906/1‑149 (MQ=255)
                                                                                                  tCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTAATGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCa                                                    >  6:291054/1‑147 (MQ=255)
                                                                                                      ggCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTAACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCg                                              >  1:332700/1‑149 (MQ=255)
                                                                                                             gcgcTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCa                                       <  6:312412/149‑1 (MQ=255)
                                                                                                              ctcTCGCGTTCGGAATGCTGGTTCTCCAG‑‑GGCCTGGCGCAACGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTTCGCTTGTTTCGGCATAGCCGGTCAGCaa                                      <  5:291054/147‑1 (MQ=255)
                                                                                                               gcTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACTGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGt                                                         <  6:285230/129‑1 (MQ=255)
                                                                                                               gcTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGt                                                         >  5:285230/1‑129 (MQ=255)
                                                                                                                 tCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGc                                        >  1:193466/1‑144 (MQ=255)
                                                                                                                 tCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGc                                        <  4:146025/144‑1 (MQ=255)
                                                                                                                 tCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGc                                        >  3:146025/1‑144 (MQ=255)
                                                                                                                   gcgTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATcc                                 >  4:172049/1‑149 (MQ=255)
                                                                                                                     gTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCTCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCt                                <  8:306675/148‑1 (MQ=255)
                                                                                                                        cGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGAc                             <  6:130537/148‑1 (MQ=255)
                                                                                                                              gCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGAt                      >  8:58227/1‑149 (MQ=255)
                                                                                                                                tGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATgg                    >  3:305993/1‑149 (MQ=255)
                                                                                                                                  ttGCGCCAG‑‑GCCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGt                  >  4:263149/2‑149 (MQ=255)
                                                                                                                                  ttGCGCCAG‑‑GCCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGt                  >  4:258680/2‑149 (MQ=255)
                                                                                                                                  gTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATgg                    <  8:140906/147‑1 (MQ=255)
                                                                                                                                  gTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATgg                    <  2:193466/147‑1 (MQ=255)
                                                                                                                                  gTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGt                  >  5:45435/1‑149 (MQ=255)
                                                                                                                                  gTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGt                  <  7:58227/149‑1 (MQ=255)
                                                                                                                                     cgccAG‑‑GCCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGa                <  3:258680/146‑1 (MQ=255)
                                                                                                                                     cgccAG‑‑GCCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGa                <  3:263149/146‑1 (MQ=255)
                                                                                                                                      tCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGAt               >  7:220720/1‑148 (MQ=255)
                                                                                                                                       ccAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGc             >  8:99177/1‑149 (MQ=255)
                                                                                                                                       ccAG‑‑GGCCTGGCACAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGa                <  2:128393/146‑1 (MQ=255)
                                                                                                                                       ccAG‑‑GGCCTGGCACAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGa                >  1:128393/1‑146 (MQ=255)
                                                                                                                                                    cgcAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGTATCTGAGCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCtg  <  2:249315/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTG  >  NC_002947/3345977‑3346271

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: