Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,954,658 A→G 11.6% intergenic (‑161/‑120) PP_2586 ← / → PP_2587 putative amino acid transporter/LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,954,6580AG11.6% 152.4 / 9.3 61intergenic (‑161/‑120)PP_2586/PP_2587putative amino acid transporter/LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base G (2/5);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATG  >  NC_002947/2954521‑2954806
                                                                                                                                          |                                                                                                                                                     
aGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCtt                                                                                                                                             >  8:104632/1‑148 (MQ=255)
aGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCtt                                                                                                                                             <  7:104632/148‑1 (MQ=255)
aGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTg                                                                                                                                            >  6:349542/1‑149 (MQ=255)
aGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTg                                                                                                                                            >  3:153370/1‑149 (MQ=255)
 ggACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTg                                                                                                                                            <  6:368755/148‑1 (MQ=255)
  gACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTTGGGTGAGCGGGCACAGCCCATGcc                                                                                                                                          >  6:243418/1‑149 (MQ=255)
    cacCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCg                                                                                                                                        <  5:349542/149‑1 (MQ=255)
       cgcgGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCt                                                                                                                                     >  4:76398/1‑149 (MQ=255)
               gCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCa                                                                                                                              >  2:243402/1‑148 (MQ=255)
                 gAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGGCAAGGGCTGTGTCCGCTCACCCCAAAgg                                                                                                                          <  4:46864/149‑1 (MQ=37)
                   gCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCGGAATCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGt                                                                                                                        <  7:57266/149‑1 (MQ=18)
                    ccggCTCGCGACCGTCAGGTGTTGGGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGCCACCCGGAATCCCCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGt                                                                                                                        <  5:363244/147‑1 (MQ=21)
                           gcgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCg                                                                                                                <  8:142766/149‑1 (MQ=255)
                           gcgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCg                                                                                                                <  1:243402/149‑1 (MQ=255)
                           gcgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCg                                                                                                                <  5:149016/149‑1 (MQ=255)
                            cgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCgg                                                                                                               <  5:243418/149‑1 (MQ=255)
                               ccGTCAGGTGTTGTGCAAGGATGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGc                                                                                                            <  8:260303/149‑1 (MQ=255)
                                cGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGATGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTGATTCCgggtgc                                                                                                           <  1:295354/149‑5 (MQ=18)
                                  tCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCtt                                                                                                                                             <  7:90744/114‑1 (MQ=255)
                                  tCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCtt                                                                                                                                             >  8:90744/1‑114 (MQ=255)
                                  tCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTACGGTGCCTg                                                                                                         >  8:227616/1‑149 (MQ=255)
                                      gtgtTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCGGAATCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTGATTCCGGGTGCCTGCAgg                                                                                                     <  6:289162/149‑2 (MQ=14)
                                       tgttgtGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAgg                                                                                                                          <  8:124691/128‑1 (MQ=255)
                                       tgttgtGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAgg                                                                                                                          >  7:124691/1‑128 (MQ=255)
                                       tgttgtGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAgcg                                                                                                    >  6:165292/1‑149 (MQ=255)
                                           gtgCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGt                                                                                                >  1:195214/1‑149 (MQ=255)
                                            tgCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTg                                                                                               <  3:76398/149‑1 (MQ=255)
                                              cAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTgg                                                                                              >  3:83188/1‑148 (MQ=255)
                                               aaGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTTATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGAcaca                                                                                                                                                   <  5:99613/95‑1 (MQ=255)
                                               aaGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGAcaca                                                                                                                                                   >  6:99613/1‑95 (MQ=255)
                                               aaGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGt                                                                                                <  6:80007/145‑1 (MQ=255)
                                               aaGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGt                                                                                                >  5:80007/1‑145 (MQ=255)
                                                              aaaGCTCGCCGCGGTGGGGCGGTTTTATGCCCCCACCCCCGGCCGCCCCCCGTAATCCCCCTTGGGGGTTAGCGGACACCGCCCTTGCCGGGCTCACCCCAATGGGTTTTTCGGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAacga                                                                             <  4:153370/149‑1 (MQ=255)
                                                              aaaGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAacga                                                                             >  4:350370/1‑149 (MQ=255)
                                                                         cgcgGGGGCGGGTTCATTCTCCCACACCCTGCAGCCACCCGTAATCACCCTTTGGGGTGAGCGGACACAGCCCTTCCCCGGCTCACCCCAAAGTTTGATTCCGGTGCCTGCAGCGCCGTGGTAGTAACTGCTCAACGACGCACCGCACt                                                                  <  7:227616/149‑1 (MQ=255)
                                                                                 cGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGcggcgg                                                          <  4:83188/149‑1 (MQ=255)
                                                                                    gTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGCCACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGc                                                       >  3:358192/1‑149 (MQ=255)
                                                                                       gATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCAc                                                    >  4:280682/1‑149 (MQ=255)
                                                                                           cTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCCCCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGGGGGAGGATCTGGGGAACGACGCACCGCCCTGGCGGCGGTGCCCCGgag                                                >  2:288712/1‑149 (MQ=38)
                                                                                            tCCCCCCCCCTCCAGCCCCCCGCTAGCCCCCTTTGGGGGTTGCGGACACCCCCCTTGCCCGCCTCACCCCAATGGTTTTTTCGGGTGCGTGCAGCGCCGTGGGAGTATCTTCTCAACGACGCACTGCACTGGCGGCGGTGCCACTga                                                 <  7:55315/147‑1 (MQ=255)
                                                                                            tCCCACACCCTGCAGGCACCCGGAATCACCCATTGGGGTGAGCGGACACAGCCCTTGCCCGGCACACCCCAAAGGGTGATTGCGGGGCCGGCAGCGCCGGGGGGGTAGCGGCGCAACGGCGCACCGCCCTGGCGGCGGGGccccggg                                                 >  8:55315/1‑142 (MQ=255)
                                                                                              ccACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCaa                                                                                                                             <  3:248576/70‑1 (MQ=255)
                                                                                              ccACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCaa                                                                                                                             >  4:248576/1‑70 (MQ=255)
                                                                                                acccccTCCAGGCACCCGTAATCACCCTTTGGGGTGAGCGGACACCGCCCTGGCCCGGCTCACCCAAATGGGTTATTCGGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTga                                                 <  2:155272/140‑1 (MQ=255)
                                                                                                acacCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTga                                                 >  1:155272/1‑143 (MQ=255)
                                                                                                      tGCAGGCACCCGGAATCACCCTTTGGGGGGAGCGGACACAGCCCTTGCCGGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGCCGCACCGGACTGGCGGCGGTGCGACTGAGAAAACGGCAgg                                     >  1:313816/1‑149 (MQ=39)
                                                                                                           gCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTa                                 >  1:67669/1‑148 (MQ=255)
                                                                                                             aCCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAg                              <  5:165292/149‑1 (MQ=255)
                                                                                                             aCCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAg                              >  6:311416/1‑149 (MQ=255)
                                                                                                               ccGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGt                             >  8:56137/1‑148 (MQ=255)
                                                                                                                   aaTCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGc                                                       <  4:346564/118‑1 (MQ=255)
                                                                                                                   aaTCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGc                                                       >  3:346564/1‑118 (MQ=255)
                                                                                                                    aTCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACcaa                       >  7:276274/1‑149 (MQ=21)
                                                                                                                      cACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCGACTGAGAACACGGCAGGCGTACAGTGCCcaat                      >  5:55358/1‑148 (MQ=255)
                                                                                                                          cTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAAtt                 <  5:311416/149‑1 (MQ=255)
                                                                                                                           tttGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCa                                       <  7:144186/126‑1 (MQ=255)
                                                                                                                           tttGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCa                                       >  8:144186/1‑126 (MQ=255)
                                                                                                                           tttGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGt                             >  3:24981/1‑136 (MQ=255)
                                                                                                                           tttGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGt                             <  4:24981/136‑1 (MQ=255)
                                                                                                                               gggTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCttt            >  7:1827/1‑149 (MQ=255)
                                                                                                                               gggTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCtt             >  8:197863/1‑148 (MQ=255)
                                                                                                                                          acacAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATg  >  3:110062/1‑149 (MQ=255)
                                                                                                                                          acacAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATg  <  7:197863/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
AGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATG  >  NC_002947/2954521‑2954806

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: