Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,912,855 C→A 6.4% G408G (GGG→GGT)  PP_3453 ← sensor protein RstB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,912,8550CA 6.4% 202.1 / 3.6 62G408G (GGG→GGT) PP_3453sensor protein RstB
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base A (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

GGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCATCCTCCACCGT  >  NC_002947/3912731‑3912998
                                                                                                                            |                                                                                                                                                
ggCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGCTCGCCCGCCGGCTGATg                                                                                                                          >  5:27238/1‑149 (MQ=255)
               gCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGcc                                                                                                            >  3:255035/1‑148 (MQ=255)
                  ggCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCaaaa                                                                                                        >  2:287135/1‑149 (MQ=255)
                   gCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCAGGCTGATGGCCACCCCCAGGCCAAAac                                                                                                       >  3:258756/1‑149 (MQ=255)
                    cccGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAac                                                                                                       >  3:157767/1‑148 (MQ=255)
                      cGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                                    <  2:352213/149‑1 (MQ=255)
                       gTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAg                                                                                                               <  6:54042/137‑1 (MQ=255)
                       gTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAg                                                                                                               >  5:54042/1‑137 (MQ=255)
                        tctgTTCGCGGTTGCTGCGCCCCTTTCCGCCGGCGGGCAAACGATTCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCTGGCTGATGGCCAACCCCAGGCCAAAACCAC‑Cg                                                                                                 <  4:258756/149‑1 (MQ=255)
                            ttCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                                    <  1:218921/143‑1 (MQ=255)
                            ttCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                                    >  2:218921/1‑143 (MQ=255)
                            ttCGCGGTTGCTGCGCCCCTATCCGCCCGCGGGCAAACGAATCCTGAAATGCGCCCCACACAGCCCCGACTTCGCCAGGGTCAACGTCCCACCCTGCCCCACGAACGCCCCCCGGCTGATGGCCCACCCCAGGCAAAAACCAC‑CGgtgt                                                                                             >  6:236388/1‑149 (MQ=255)
                               gcgGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGc                                                                                          >  1:90191/1‑149 (MQ=255)
                                  gTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                                    >  2:144066/1‑137 (MQ=255)
                                  gTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                                    <  1:144066/137‑1 (MQ=255)
                                          gcCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGACCCACCCAGCGCGGATTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCTGCTGATGGCCAACCCCAGGCCAAAAACAC‑CGGTGTTGCGATCGCGACt                                                                                <  5:236388/148‑1 (MQ=255)
                                            cccTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCg                                                                             <  3:215655/149‑1 (MQ=255)
                                            cccTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCg                                                                             >  4:215655/1‑149 (MQ=255)
                                              cTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGt                                                                           >  1:303204/1‑149 (MQ=255)
                                                 tCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGTGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTc                                                                          >  6:96053/1‑147 (MQ=255)
                                                 tCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTc                                                                          <  2:303204/147‑1 (MQ=255)
                                                   cgccgGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                                      <  7:106563/149‑1 (MQ=255)
                                                   cgccgGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                                      <  4:255035/149‑1 (MQ=255)
                                                         cggGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                                <  6:27238/149‑1 (MQ=255)
                                                         cggGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                                <  5:255476/149‑1 (MQ=255)
                                                                  gAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTg                                                        >  8:168263/1‑148 (MQ=255)
                                                                     tCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGaa                                                     >  6:89211/1‑148 (MQ=255)
                                                                       cTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCCAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAAc                                                  <  1:287135/149‑1 (MQ=255)
                                                                            cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACTCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGc                                              <  1:96617/149‑1 (MQ=255)
                                                                            cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTa                                                                 <  2:234734/129‑1 (MQ=255)
                                                                            cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTa                                                                 >  1:234734/1‑129 (MQ=255)
                                                                               cgcCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTcgc                                          >  6:279942/1‑149 (MQ=255)
                                                                                 ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTTCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAAc                                                  <  3:70520/139‑1 (MQ=255)
                                                                                 ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGTAAAAc                                                  >  4:70520/1‑139 (MQ=255)
                                                                                 ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCgcgc                                        <  7:168263/149‑1 (MQ=255)
                                                                                  cccacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGct                                       >  1:320005/1‑149 (MQ=255)
                                                                                   ccacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctc                                      <  5:279942/149‑1 (MQ=255)
                                                                                    cacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcc                                     >  4:189842/1‑149 (MQ=255)
                                                                                       ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                                >  5:237055/1‑119 (MQ=255)
                                                                                       ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                                <  6:237055/119‑1 (MQ=255)
                                                                                       ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGACCGGCTGGAAAACCCGCTCGCGctcctc                                   >  5:189642/1‑148 (MQ=255)
                                                                                       ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCtcctcc                                  <  2:90191/149‑1 (MQ=255)
                                                                                       ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCccctcc                                  >  8:299590/1‑149 (MQ=255)
                                                                                            gcCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAgg                             <  2:320005/149‑1 (MQ=255)
                                                                                             ccGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAggg                            >  7:214756/1‑149 (MQ=255)
                                                                                                     tcgaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                                      <  6:173484/96‑1 (MQ=38)
                                                                                                     tcgaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                                      >  5:173484/4‑99 (MQ=38)
                                                                                                     gcCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAAt                    >  4:183002/1‑149 (MQ=255)
                                                                                                      ccAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATc                   <  4:157767/149‑1 (MQ=255)
                                                                                                      ccAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCCCCTCCACAGCGACGCCAATc                   <  7:299590/149‑1 (MQ=255)
                                                                                                       gaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcctc                                   >  7:37015/2‑132 (MQ=255)
                                                                                                       gaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcctc                                   <  8:37015/131‑1 (MQ=255)
                                                                                                       cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAg                              >  6:311444/1‑137 (MQ=255)
                                                                                                       cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAg                              <  5:311444/137‑1 (MQ=255)
                                                                                                       cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATcc                  >  2:118367/1‑149 (MQ=255)
                                                                                                            tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCa             <  5:96053/149‑1 (MQ=255)
                                                                                                            tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCa             <  5:89211/149‑1 (MQ=255)
                                                                                                            tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACACGGACGCCAATCCCGTCa             >  3:263468/1‑149 (MQ=255)
                                                                                                              aaGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCAtc           <  8:214756/149‑1 (MQ=255)
                                                                                                                ggTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATc                   <  6:180935/139‑1 (MQ=255)
                                                                                                                ggTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATc                   >  5:180935/1‑139 (MQ=255)
                                                                                                                        ccTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCATCCTCCACCGt  >  6:86685/1‑148 (MQ=255)
                                                                                                                            |                                                                                                                                                
GGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCATCCTCCACCGT  >  NC_002947/3912731‑3912998

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