Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,591,716 A→C 8.6% C60G (TGC→GGC)  PP_1395 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,591,7160AC 8.6% 167.2 / 6.9 58C60G (TGC→GGC) PP_1395AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (33/20);  new base C (3/2);  total (36/22)
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.89e-01

CGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTCCAGG  >  NC_002947/1591568‑1591864
                                                                                                                                                    |                                                                                                                                                    
cgcAGGCCCAGGCTTGGGGTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGca                                                                                                                                                      <  8:169864/149‑1 (MQ=255)
       ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAggg                                                                                                                                                 >  6:203034/1‑147 (MQ=255)
       ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAggg                                                                                                                                                 >  2:26723/1‑147 (MQ=255)
                 gTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTgcgc                                                                                                                                       >  3:189990/1‑147 (MQ=255)
                     tcGGGGCTGGCATCTAGCCTCAAGGCGGTGAGTGGCAATGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggt                                                                                                                                 <  8:340289/149‑1 (MQ=255)
                     tcGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggt                                                                                                                                 >  6:281251/1‑149 (MQ=255)
                       ggggCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggt                                                                                                                                 <  3:232783/147‑1 (MQ=255)
                        gggCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggtggg                                                                                                                              <  1:26723/149‑1 (MQ=255)
                           cTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTa                                                                                                                            >  8:63667/1‑148 (MQ=255)
                              gCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACa                                                                                                                          <  5:203034/147‑1 (MQ=255)
                                           cGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCg                                                                                                           <  6:402526/149‑1 (MQ=255)
                                             gcgcTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                                         <  7:63667/149‑1 (MQ=255)
                                             gcgcTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                                         <  7:145248/149‑1 (MQ=255)
                                                 tGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTaa                                                                                                     >  7:356218/1‑149 (MQ=255)
                                                       gCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCAc                                                                                                >  1:402237/1‑148 (MQ=255)
                                                        cAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCa                                                                                              <  4:24119/149‑1 (MQ=255)
                                                               aCCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAg                                                                                       >  7:308672/1‑149 (MQ=255)
                                                                 cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAgg                                                                                            <  4:224822/142‑1 (MQ=255)
                                                                 cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAgg                                                                                            >  3:224822/1‑142 (MQ=255)
                                                                 cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAgg                                                                                      >  5:404030/1‑148 (MQ=255)
                                                                 cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGt                                                                                     >  4:198308/1‑149 (MQ=255)
                                                                         ccATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                                             >  1:168515/1‑149 (MQ=255)
                                                                         ccATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                                             <  5:281251/149‑1 (MQ=255)
                                                                                      tGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGAc                                                                <  8:356218/149‑1 (MQ=255)
                                                                                      tGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTG‑CCTTGGCAAAGAAGCACAAGGCAGGCAGCGCCAGGGTTGGTACAGAAAGTACATTTTAGGTGTTACGCACCAGGGACAAGTCATAGAATGCCCAAGCTTGAc                                                                >  2:136490/1‑148 (MQ=255)
                                                                                          tCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCa                                                              >  6:277304/1‑147 (MQ=255)
                                                                                          tCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCa                                                              >  2:93873/1‑147 (MQ=255)
                                                                                              aGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGc                                                        >  8:327011/1‑149 (MQ=255)
                                                                                                gCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTc                                                      >  2:229806/1‑149 (MQ=255)
                                                                                                     cTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgg                                                                                                                                   >  4:265967/1‑67 (MQ=255)
                                                                                                     cTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgg                                                                                                                                   <  3:265967/67‑1 (MQ=255)
                                                                                                         tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                                             <  2:278900/117‑1 (MQ=255)
                                                                                                         tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                                             >  1:278900/1‑117 (MQ=255)
                                                                                                         tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATg                                             >  2:32615/1‑149 (MQ=255)
                                                                                                         tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATg                                             >  3:95082/1‑149 (MQ=255)
                                                                                                          cGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATgg                                            >  1:358339/1‑149 (MQ=255)
                                                                                                           gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTAc                                                                                                                           >  2:328632/1‑69 (MQ=255)
                                                                                                           gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTAc                                                                                                                           <  1:328632/69‑1 (MQ=255)
                                                                                                           gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGCCCTAAGGGCTCGCCCGCAAGGc                                           >  8:384333/1‑149 (MQ=255)
                                                                                                           gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGc                                           <  5:213334/149‑1 (MQ=255)
                                                                                                              aCCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCgg                                         <  4:189990/148‑1 (MQ=255)
                                                                                                                ccAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCACCCGCATGGCGGtt                                       >  1:236927/1‑148 (MQ=255)
                                                                                                                    gAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCa                                  <  2:168515/149‑1 (MQ=255)
                                                                                                                        aGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTc                              >  6:204831/1‑149 (MQ=255)
                                                                                                                          gcgcagcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCgg                            >  6:335432/8‑149 (MQ=38)
                                                                                                                          gCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCgg                            >  4:1846/1‑149 (MQ=255)
                                                                                                                           cgcagcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGc                           <  5:335432/143‑1 (MQ=255)
                                                                                                                            ttttGCTGCCCTGGGCAAGGATGCACAGGGCAGGCAGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGc                           >  1:268745/1‑148 (MQ=255)
                                                                                                                              agcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAg                        >  6:172323/4‑149 (MQ=255)
                                                                                                                                gctgcCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCt                      >  7:29436/1‑149 (MQ=255)
                                                                                                                                cctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTg                                                                   >  7:200816/2‑104 (MQ=38)
                                                                                                                                cctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTg                                                                   <  8:200816/103‑1 (MQ=38)
                                                                                                                                  tgcCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTc                     >  6:75856/1‑148 (MQ=255)
                                                                                                                                           cAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGCTCAGCTCGGCCAGCTCCAGGCGCTg            >  8:92260/1‑148 (MQ=255)
                                                                                                                                           cAAGGATGCACAGGGCAGGCTGCACCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCaaa                                                            >  5:357521/1‑100 (MQ=255)
                                                                                                                                           cAAGGATGCACAGGGCAGGCTGCACCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCaaa                                                            <  6:357521/100‑1 (MQ=255)
                                                                                                                                            aaGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTg            <  4:196091/147‑1 (MQ=255)
                                                                                                                                            aaGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTg            >  3:196091/1‑147 (MQ=255)
                                                                                                                                              ggATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTCATCAGCTCGGCCAGCTCCAGGCGCTGCAGt        <  1:32615/149‑1 (MQ=255)
                                                                                                                                                   cacaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTCCAg   <  2:358339/149‑1 (MQ=255)
                                                                                                                                                   cacaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCACCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTCCAg   <  2:236927/149‑1 (MQ=255)
                                                                                                                                                    acaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTCCAgg  <  2:402237/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTCCAGG  >  NC_002947/1591568‑1591864

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