Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 222,265 A→G 5.3% Q40Q (CAA→CAG)  PP_0170 → ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947222,2650AG 5.3% 277.6 / 3.0 93Q40Q (CAA→CAG) PP_0170ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (45/42);  new base G (2/3);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CCCGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCG  >  NC_002947/222119‑222412
                                                                                                                                                    |                                                                                                                                                   
cccGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAAtt                                                                                                                                                   >  5:349034/1‑149 (MQ=255)
   gTTCTGACGACTGAGGAAACCATCCTGCGCAACTCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCGGGGCGCAAGCCTGGTCGTCGG‑CAGCCGAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCgg                                                                                                                                                <  4:21596/149‑1 (MQ=255)
   gTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCgg                                                                                                                                                <  6:357631/149‑1 (MQ=255)
    ttCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGc                                                                                                                                               <  5:143280/149‑1 (MQ=255)
         acgacTGAGGAAACCATCATGCGCAAATCCATCAGCCGGCTGGCGGCAAGCATC‑GGCCTGTGCGCAAGCCTGGTCGTCGT‑GAGCCTAGCCGCCCCCGCCGTTGCCGAGGCCTAAGGCAAGATCCGTATTGCCGAACAATTCGGCAtg                                                                                                                                            <  4:331903/147‑2 (MQ=255)
         acgacTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAATTGTTCGGCaata                                                                                                                                           <  5:363921/148‑4 (MQ=255)
         acgacTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAATTGTTCGGCaata                                                                                                                                           <  2:311929/148‑4 (MQ=255)
         acgacTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAATTGTTCGGCaa                                                                                                                                             <  3:379120/146‑2 (MQ=255)
             cTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGc                                                                                                                                               <  5:191535/140‑1 (MQ=255)
             cTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGc                                                                                                                                               >  6:191535/1‑140 (MQ=255)
                aGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACctg                                                                                                                                   >  5:356602/1‑149 (MQ=255)
                aGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACctg                                                                                                                                   >  1:141192/1‑149 (MQ=255)
                 ggAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACctgc                                                                                                                                  <  4:18944/149‑1 (MQ=255)
                  gAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGGCCTGGGCGCAAGCCTGGACGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGAATTGCCGAACAATTCGGCATCGTCTACctgc                                                                                                                                  >  1:85488/1‑149 (MQ=255)
                  gAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCACAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACctgct                                                                                                                                 >  6:10904/1‑149 (MQ=255)
                  gAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTCGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACcagct                                                                                                                                 >  5:165377/1‑149 (MQ=255)
                     aCCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACctgctgct                                                                                                                              >  4:160713/1‑149 (MQ=255)
                       catcatGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACctgctgctg                                                                                                                             >  3:379809/1‑148 (MQ=255)
                           atGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACgt                                                                                                                        >  4:101163/1‑149 (MQ=255)
                              cgcAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACgtggt                                                                                                                     <  1:110756/149‑1 (MQ=255)
                               gcAAATCCATCAGCCGACTTGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑AAGCCTAGCCGCCCCCGCCGTGGCAAAGGCAGAAGGCAATATACCTATTGCCGAACAATTCGGCATCGTATACCTTCTGCTGAACgtggtg                                                                                                                    >  4:51015/1‑149 (MQ=255)
                                cAAATCCATCAGCCGCCTGGGGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGc                                                                                                                   <  5:56496/149‑1 (MQ=255)
                                    tCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAATGTGGTGCGTGa                                                                                                               >  5:278594/1‑149 (MQ=255)
                                    tCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGa                                                                                                               <  2:85488/149‑1 (MQ=255)
                                      cATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTCCTGCTGAACGTGGTGCGTGAc                                                                                                              <  1:94145/148‑1 (MQ=255)
                                       aTCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGAc                                                                                                              >  2:6726/1‑147 (MQ=255)
                                       aTCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGAc                                                                                                              >  3:294003/1‑147 (MQ=255)
                                           gccgccTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCa                                                                                                            <  6:199124/145‑1 (MQ=255)
                                           gccgccTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCa                                                                                                            >  5:199124/1‑145 (MQ=255)
                                           gcagccTGCCGGCAAGCCTC‑GGCCTTGGCGCAATCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCAc                                                                                                        <  8:346930/146‑1 (MQ=255)
                                            ccgccTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACACTTTGGCATCGTCTACCTGCTGGTGCACGTGGTGCGTGACCAGCAcc                                                                                                       >  4:70560/1‑149 (MQ=255)
                                             cgccTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTg                                                                                                                >  5:143641/1‑139 (MQ=255)
                                             cgccTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTg                                                                                                                <  6:143641/139‑1 (MQ=255)
                                               ccTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGa                                                                                                    <  6:349034/149‑1 (MQ=255)
                                               ccTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGa                                                                                                    <  3:116541/149‑1 (MQ=255)
                                               ccTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGACATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGa                                                                                                    >  5:134359/1‑149 (MQ=255)
                                                    cggcAAGCATC‑GGCCTGGTCTCAAGCCTGGTCGTCGG‑CAGACTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAAGAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTgag                                                                                               <  6:301903/149‑1 (MQ=255)
                                                    cggcAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGAtt                                                                                                  >  1:243771/1‑146 (MQ=255)
                                                      gcAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGGTCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTgag                                                                                               >  1:58357/1‑148 (MQ=255)
                                                      gcAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGaa                                                                                             <  3:284941/149‑1 (MQ=255)
                                                       cAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTTCGGAACAGTTCGGCATCGTCTCACTGCTGGGGACCGTGGTGCGTGACCAGCAGCTGCTTggga                                                                                              >  4:293719/1‑147 (MQ=255)
                                                        aaGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGaaa                                                                                            >  2:213089/1‑148 (MQ=255)
                                                          gCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAacac                                                                                         <  4:128536/149‑1 (MQ=255)
                                                            aTC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACgg                                                                                       <  8:327029/149‑1 (MQ=255)
                                                            aTC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACgg                                                                                       >  1:326981/1‑149 (MQ=255)
                                                            aTC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCACGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCGACCTGCTGGTGAACGTGGTGCGTGACCAGCACCTGATTGAGATACACgg                                                                                       >  7:204919/1‑149 (MQ=255)
                                                                ggCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGc                                                                                      >  4:81750/1‑147 (MQ=255)
                                                                 gCCTGGGCTCAAGCCTGGTCGTCTG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTTACCAGCACCTGATTGAGAAACACTGCAAg                                                                                   <  3:51015/149‑1 (MQ=255)
                                                                      ggCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTga                                                                                                <  8:33981/131‑1 (MQ=255)
                                                                      ggCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTga                                                                                                >  7:33981/1‑131 (MQ=255)
                                                                        cgcAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAgg                                                                            <  4:379809/149‑1 (MQ=255)
                                                                         gcAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAggg                                                                           >  2:93355/1‑149 (MQ=255)
                                                                         gcAAGCCTGGTCGTCGG‑CAGCCTAGCCGCACCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCAGAACAATTCGGCATCGTCTTCCTGCTGCTGAACGTGGTGCGTGACCAGAACCTGATTGAGAAACACGGCAAGGAGCAggg                                                                           >  2:311365/1‑149 (MQ=255)
                                                                            aGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCGAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGcat                                                                        >  4:172356/1‑149 (MQ=255)
                                                                                tGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGcatcga                                                                     >  2:293304/1‑148 (MQ=255)
                                                                                   tcgtcgG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGcatcgacatc                                                                 >  1:78941/1‑149 (MQ=255)
                                                                                     gtcgG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGcatcgacatcg                                                                <  1:213089/148‑1 (MQ=255)
                                                                                      tcgG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGaa                                                              <  3:81750/149‑1 (MQ=255)
                                                                                      tcgG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGaa                                                              <  6:356602/149‑1 (MQ=255)
                                                                                      tcgG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGaa                                                              <  3:160713/149‑1 (MQ=255)
                                                                                      gcgG‑CAGCCTAGCCGTCCCCGCCGTGGCCCAGGCCGAAGGCAAGTTCCGTATTGCCGAACAATTGGGCATCGTCTACCTGCTGCTGTACGTGGTGCGTGTCGTGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGaa                                                              <  5:321721/148‑1 (MQ=255)
                                                                                              ccTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGAtt                                                       <  4:294003/149‑1 (MQ=255)
                                                                                              ccTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGAtt                                                       <  2:141192/149‑1 (MQ=255)
                                                                                              ccTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGAtt                                                       <  2:58357/149‑1 (MQ=255)
                                                                                              ccTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGAtt                                                       <  5:10904/149‑1 (MQ=255)
                                                                                               cTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTg                                                      >  6:331579/1‑149 (MQ=255)
                                                                                                tAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAAATCGGCATCGTCTACCTGCTGCTGAACGTGGTGCCTGACCCGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTgg                                                     >  7:284474/1‑149 (MQ=255)
                                                                                                  gccgccCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGgc                                                   >  5:195148/1‑149 (MQ=255)
                                                                                                     gccCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCa                                                <  2:326981/149‑1 (MQ=255)
                                                                                                      cccccGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAATGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAg                                               <  6:278594/149‑1 (MQ=255)
                                                                                                        cccGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCt                                             >  4:134608/1‑149 (MQ=255)
                                                                                                               tGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAATCACGGGAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGt                                      <  6:165377/149‑1 (MQ=255)
                                                                                                               tGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGt                                                            <  1:390007/127‑1 (MQ=255)
                                                                                                               tGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGt                                      <  3:101163/149‑1 (MQ=255)
                                                                                                               tGGCCCAGGACGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGt                                                            >  2:390007/1‑127 (MQ=255)
                                                                                                               tGGCCAAGGCCGAAGGCAAGAGCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGCACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGt                                      <  3:293719/149‑1 (MQ=255)
                                                                                                                 gCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTgg                                    <  6:195148/149‑1 (MQ=255)
                                                                                                                  cccAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTgg                                    >  8:273939/1‑148 (MQ=255)
                                                                                                                   ccAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCt                                  >  1:197913/1‑149 (MQ=255)
                                                                                                                   ccAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCt                                  >  8:293654/1‑149 (MQ=255)
                                                                                                                       gCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGc                              >  1:99231/1‑149 (MQ=255)
                                                                                                                              gCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAAc                       >  2:24312/1‑149 (MQ=255)
                                                                                                                                aaGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGa                     >  3:309535/1‑149 (MQ=255)
                                                                                                                                      cGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGTTCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCAt                <  1:311365/148‑1 (MQ=255)
                                                                                                                                         gatGCCGAATTGTTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCAttgtt            >  4:102170/3‑149 (MQ=255)
                                                                                                                                          ttGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCAttgttg           >  3:282073/1‑149 (MQ=255)
                                                                                                                                          ttGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCAttgtt            <  1:6726/148‑1 (MQ=255)
                                                                                                                                             ccGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTga                                                                                                >  5:53672/1‑61 (MQ=255)
                                                                                                                                             ccGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTga                                                                                                <  6:53672/61‑1 (MQ=255)
                                                                                                                                             ccGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTcc        >  3:227938/1‑149 (MQ=255)
                                                                                                                                              cGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGTAtt               >  4:267644/1‑141 (MQ=255)
                                                                                                                                              cGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGTAtt               <  3:267644/141‑1 (MQ=255)
                                                                                                                                               gAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCgg      >  1:12047/1‑149 (MQ=255)
                                                                                                                                                   aaTTCGGCATCGTCTACCTGCTGCTTAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCg  <  2:99231/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATC‑GGCCTGGGCGCAAGCCTGGTCGTCGG‑CAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCG  >  NC_002947/222119‑222412

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