Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 568,254 G→C 7.2% R414R (CGC→CGG uvrA ← uvrABC excinuclease protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947568,2540GC 7.2% 213.1 / 6.5 69R414R (CGC→CGGuvrAuvrABC excinuclease protein A
Reads supporting (aligned to +/- strand):  ref base G (37/27);  new base C (3/2);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

GCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGACGAT  >  NC_002947/568107‑568400
                                                                                                                                                   |                                                                                                                                                   
gCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCgcggc                                                                                                                                                    >  7:225675/1‑149 (MQ=255)
gCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCgcggc                                                                                                                                                    <  6:130768/149‑1 (MQ=255)
  aGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGgcgc                                                                                                                                                  <  5:303657/149‑1 (MQ=255)
  aGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTGGCTGGCTTCACCGATCGGCAAATTGGTTACCGCAGGCAGGGTTTACTCGCCCCCCCACACATGGCGCGCCTCgcggcg                                                                                                                                                   >  6:179213/1‑148 (MQ=255)
       gccTTCAGCATTTTTGCCGCCATTCCGCCACGCCGGCCAGTAAGGGCCAGTGCCCCAAAGTAGTGGCGGCCTTCACCGATCTGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGgc                                                                                                                                             <  1:313964/147‑1 (MQ=255)
         cTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTGCGCCGCGAGGcgcg                                                                                                                                          <  3:137515/149‑1 (MQ=255)
          ttGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGcgcg                                                                                                                                          >  7:243307/1‑148 (MQ=255)
           tGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTCCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTa                                                                                                                                        >  6:148476/1‑149 (MQ=255)
             aGGATTTTTGCCGCGATTTCGCCCCGCCGGCCATTAAGGGGCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGACGCGCCTGCGCCGCGAGGcgcgcca                                                                                                                                       <  1:304495/148‑4 (MQ=255)
                                cGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACgg                                                                                                                    <  3:292089/148‑1 (MQ=255)
                                 gCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGtt                                                                                                                  >  3:388409/1‑149 (MQ=255)
                                   cACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCTCGGTTGc                                                                                                                <  8:154323/149‑1 (MQ=255)
                                           ccAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAg                                                                                                        <  2:161891/149‑1 (MQ=255)
                                           ccAGTAAGGGGCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAg                                                                                                        <  7:134035/149‑1 (MQ=255)
                                              gTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGaa                                                                                                     <  2:249078/149‑1 (MQ=255)
                                              gTAAGCGGCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGaa                                                                                                     <  4:282462/149‑1 (MQ=255)
                                               tAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAAt                                                                                                    >  5:38814/1‑149 (MQ=255)
                                               tAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAAt                                                                                                    >  7:221455/1‑149 (MQ=255)
                                               tAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAAt                                                                                                    >  7:48439/1‑149 (MQ=255)
                                                  gggTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTg                                                                                                 <  4:388409/149‑1 (MQ=255)
                                                  gggTCAGTTGACCAAAGTAGTTCCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGACACGGTTGCGTGCCAAGGAATTTg                                                                                                 <  5:148476/149‑1 (MQ=255)
                                                           gACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGttct                                                                                        >  2:35696/1‑149 (MQ=255)
                                                           gACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGttc                                                                                         <  3:56467/148‑1 (MQ=255)
                                                            aCCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGttctt                                                                                       >  1:187423/1‑149 (MQ=255)
                                                              cAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTca                                                                                     >  8:221581/1‑149 (MQ=255)
                                                               aaaGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTcac                                                                                    >  2:365661/1‑149 (MQ=255)
                                                                aaGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTcacg                                                                                   >  1:354163/1‑149 (MQ=255)
                                                                aaGTAGTTGCTGGCTTCACAGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATCGAGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCAAGCTTTCGTGCAAAGGAATTTGGCCAGTTCTTcacg                                                                                   >  7:346362/1‑149 (MQ=255)
                                                                           ggCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGACAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCg                                                                        >  8:82087/1‑149 (MQ=255)
                                                                             cTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGAt                                                                      <  8:225675/149‑1 (MQ=255)
                                                                             cTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGAt                                                                      >  5:385641/1‑149 (MQ=255)
                                                                                cACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTc                                                                    <  2:354163/148‑1 (MQ=255)
                                                                                 aCCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCgg                                                                  >  5:280987/1‑149 (MQ=255)
                                                                                  ccGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGt                                                                 >  3:57919/1‑149 (MQ=255)
                                                                                   cGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGc                                                                          <  6:242776/139‑1 (MQ=255)
                                                                                   cGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGc                                                                          >  5:242776/1‑139 (MQ=255)
                                                                                     aTCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTc                                                              >  1:59993/1‑149 (MQ=255)
                                                                                     aTCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTc                                                              >  5:25636/1‑149 (MQ=255)
                                                                                        ggCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGGCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTcgcg                                                           >  3:356465/1‑149 (MQ=255)
                                                                                         gCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCgg                                                          >  2:296574/1‑149 (MQ=255)
                                                                                           aaaTTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCgg                                                          >  6:180253/1‑147 (MQ=255)
                                                                                           aaaTTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTa                                                        <  2:187423/149‑1 (MQ=255)
                                                                                            aaTTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAg                                                       >  5:71442/1‑149 (MQ=255)
                                                                                                     ccgccgGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTCGCCGATTCGGTTTCGCGGTAGCGGCGCTCa                                              <  1:35696/149‑1 (MQ=255)
                                                                                                          ggCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGtt                                         <  5:180253/149‑1 (MQ=255)
                                                                                                            cAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCg                                       <  2:59993/149‑1 (MQ=255)
                                                                                                             aGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCgg                                      <  6:38814/149‑1 (MQ=255)
                                                                                                               ggTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGTACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGca                                    >  1:315381/1‑149 (MQ=255)
                                                                                                                   ttCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAAt                                >  6:198509/1‑149 (MQ=255)
                                                                                                                   ttCTCGCCCACCCACACATGGCGCGCCTCGCGGAGCAGGCGCGTACCCCGGCAATCCGCGACCGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACTCACGGTCGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAAt                                >  3:298445/1‑149 (MQ=255)
                                                                                                                    tctcGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAAt                                <  2:164876/148‑1 (MQ=255)
                                                                                                                        gCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCCCGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACaa                                 <  7:206578/143‑1 (MQ=255)
                                                                                                                        gCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACaa                                 >  8:206578/1‑143 (MQ=255)
                                                                                                                        gCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCtt                           >  5:82715/1‑149 (MQ=255)
                                                                                                                        gCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCtt                           >  3:72387/1‑149 (MQ=255)
                                                                                                                         cccacccaCACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTc                          <  8:243307/149‑1 (MQ=255)
                                                                                                                            acccacACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGaa                       >  4:387760/1‑149 (MQ=255)
                                                                                                                            acccacACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGCCCGATACGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGaa                       >  3:312779/1‑149 (MQ=255)
                                                                                                                              ccacACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCCCGGTTGCGTGCCAAGGAAGTTGCCCAGTTCTACCCGCACGGTGGCCCATTCGGTTTCCCGGTAGCGCCGCTCAAGGTTCGGCACAATCCCTGCGATCg                     >  5:311972/1‑149 (MQ=255)
                                                                                                                              ccacACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACg                     >  3:350436/1‑149 (MQ=255)
                                                                                                                              ccacACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAAACCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACg                     >  2:280730/1‑149 (MQ=255)
                                                                                                                               cacacaTGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACgg                    <  6:280987/149‑1 (MQ=255)
                                                                                                                                acacaTGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACggg                   >  5:59347/1‑149 (MQ=255)
                                                                                                                                  acaTGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTg                 >  5:289134/1‑149 (MQ=255)
                                                                                                                                  acaTGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTg                 <  1:365661/149‑1 (MQ=255)
                                                                                                                                   caTGGCGCGGCTCGCGGCGAAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGTCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGc                <  8:346362/149‑1 (MQ=255)
                                                                                                                                      tacgcgcCTGCGCCGCGAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGAt             >  2:99957/3‑149 (MQ=255)
                                                                                                                                      tacgcgcCTGCGCCGCGAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGAt             >  5:228307/3‑149 (MQ=255)
                                                                                                                                      ggCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGAt             >  8:96998/1‑149 (MQ=255)
                                                                                                                                       gcgcgcCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGAt             <  1:296574/148‑1 (MQ=255)
                                                                                                                                       acgcgcCTGCGCCGCGAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATc            >  7:238289/2‑149 (MQ=255)
                                                                                                                                         gcgcCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGc          <  6:385641/149‑1 (MQ=255)
                                                                                                                                             cTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGa      <  3:387760/149‑1 (MQ=255)
                                                                                                                                                 cggcgCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGACGAt  <  6:25636/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGCAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGACGAT  >  NC_002947/568107‑568400

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: