Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,896,399 T→C 12.2% intergenic (+32/‑176) PP_4304 → / → sbp‑I VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,896,3990TC12.2% 143.7 / 9.7 58intergenic (+32/‑176)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (32/18);  new base C (3/4);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

AATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCC  >  NC_002947/4896251‑4896546
                                                                                                                                                    |                                                                                                                                                         
aaTGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAt                                                                                                                                                           <  3:168594/149‑1 (MQ=255)
   gCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATCGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCg                                                                                                                                                       >  5:191998/1‑149 (MQ=255)
    ccATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGc                                                                                                                                                      <  4:129974/149‑1 (MQ=255)
     cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCt                                                                                                                                                     >  2:277120/1‑149 (MQ=255)
                           caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTTAATCCGGTGCTACAGGCGCGGACCGTgcgc                                                                                                                              >  6:210749/1‑149 (MQ=255)
                            ttgttcttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                                                             <  5:385391‑M1/50‑1 (MQ=255)
                            aTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                                                             <  3:136236/149‑1 (MQ=255)
                             tGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCACATAGCGGGTAAATCAGCTACTCAAGGCGCGGGCCCTGCGCgg                                                                                                                            >  6:229441/1‑149 (MQ=255)
                               aTTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTATTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                          >  1:224843/1‑149 (MQ=255)
                                tcttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc                                                                                                                         <  4:290098‑M1/54‑1 (MQ=255)
                                     aagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGt                                                                                                                 >  6:385391‑M1/91‑149 (MQ=255)
                                           gCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                               >  6:300240/1‑145 (MQ=255)
                                           gCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                               <  5:300240/145‑1 (MQ=255)
                                           gCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGAGGAGc                                                                                                           >  2:269043/1‑149 (MQ=255)
                                             cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                         <  6:250285/149‑1 (MQ=255)
                                              ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                         >  3:290098‑M1/82‑148 (MQ=255)
                                              ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                         >  5:373065‑M1/82‑148 (MQ=255)
                                              acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                        <  7:113116/149‑1 (MQ=255)
                                               gggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                       >  3:200255‑M1/81‑149 (MQ=255)
                                                aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGATTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGGGCGCGGGCC‑‑‑CCCGTTGGGGCcgggt                                                                                                      >  7:336250/1‑147 (MQ=255)
                                                  gaagtCCCACGAGCAGGGGGCGGCTTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCGCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTTCAg                                                                                                                                              <  6:283708/111‑1 (MQ=39)
                                                  gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAAAGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAg                                                                                                                                              >  5:283708/1‑111 (MQ=255)
                                                     gaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTAccc                                                                                                 <  1:352021/149‑1 (MQ=17)
                                                        ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCcgc                                                                                               >  3:262878/1‑148 (MQ=255)
                                                          aCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTTAATACGCGCCTACAGGCGCGGACCGGGCGCGGTGC‑‑‑CCCGGAGGAGCGGGGTTAACCGCGa                                                                                             >  1:204250/1‑148 (MQ=255)
                                                            gtgCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGATTGGGTTGTTTCCCATAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAaga                                                                                          <  6:191998/147‑1 (MQ=255)
                                                               cAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGc                                                                                       <  1:277120/149‑1 (MQ=14)
                                                               cAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATCGGTTGTTTCCCAAACCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                       <  5:210749/149‑1 (MQ=255)
                                                                   taacgggtgggtctatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                     >  5:300253‑M1/61‑147 (MQ=255)
                                                                      cggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCCTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcaca                                                                                >  6:127196/1‑149 (MQ=255)
                                                                             gtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                         <  8:6895‑M1/67‑1 (MQ=255)
                                                                             gtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                         >  7:6895‑M1/51‑117 (MQ=255)
                                                                                tgctcGCTGTGGGAATGGGGTGTTGCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGc                                                                      <  6:69904/146‑1 (MQ=255)
                                                                                   tcGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATaa                                                                  >  6:82234/1‑149 (MQ=255)
                                                                                   tcGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTAGCCGGTAAATACGCACCTACAGGCCCGGACCCTGCGCGGTCC‑‑‑CCCGGAGGAGCCGGGTTACCCGCCAAGAGGCCGGCACACGCACAACCaaa                                                                   >  4:130092/1‑146 (MQ=255)
                                                                                      cTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAAcc                                                                >  7:341228/1‑148 (MQ=255)
                                                                                             aaTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtt                                                        <  5:82234/149‑1 (MQ=255)
                                                                                                gggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAttttt                                                     >  7:216428/1‑149 (MQ=255)
                                                                                                gggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAttttt                                                     >  7:210868/1‑149 (MQ=255)
                                                                                                gggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTAGCGGGTAAATACCCTACAACAGGCCCGGACCCTGCGCGGTGC‑‑‑CCCGTAGGAGCGGGTTTACCCCCGAAGAGGCCGCCACAGGCCGAACCATaaccaaaccaaattt                                                     >  4:296158/1‑144 (MQ=25)
                                                                                                   gtttttCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACACGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGt                                                  <  8:336250/146‑1 (MQ=255)
                                                                                                          ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTAGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGCAGCATAACCAAAGCATTTTTTGTTCTTTaa                                           >  8:242349/1‑149 (MQ=255)
                                                                                                           cAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                          >  8:252183/1‑149 (MQ=255)
                                                                                                           cAAAGCGTGATAGCGCCTGCGCCGGCCGCTTAGCCGGGTAATACGCGCCGACAGGCGCGGCCCCGGGGGGGGGCC‑‑‑CCCGAGGGGCCGGGTTTACCCCGGAGAGGGCGGGCCCGGCCGAAGCCAAAACCAAGACATTTTTTTTGTTTTaa                                           >  2:135233/1‑149 (MQ=255)
                                                                                                                cGTGATCGCGCGTGCGCGGGCCTCTTGGGGGGAAATTCCGCTCCTACCGGCGGGTACCGGGCGCGGCCC‑‑‑CCCGTAGTTGGGGGTTTACCGGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGta                                     <  2:204250/149‑1 (MQ=255)
                                                                                                                 gTGATAGCGCCTGCGCCGGCCAATTAGCGGGTAAATCCGCTACAAAAGGCCCGGCCCCGGCCCCCTAC‑‑‑CCCCGAGGAGCCGGTTTAACCGCGAAGAGGCCGGCACCAGCACCAGCAAAACCCAAACATTTTTTGTTCTTTAAACCgaa                                     >  4:219008/1‑146 (MQ=255)
                                                                                                                                 CGGCCTCTTCGCGGGTAAATCCGCTCCTTCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                    >  3:130092‑M1/1‑149 (MQ=255)
                                                                                                                                 CGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                    >  5:127196‑M1/1‑149 (MQ=255)
                                                                                                                                 CGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                    >  8:341228‑M1/1‑149 (MQ=255)
                                                                                                                                       CTTCTCGGGTAAATCCGCTCCTCCAGGCGCGGACCGGGCGCGGCCC‑‑‑CCCGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt              >  8:216428‑M1/1‑149 (MQ=255)
                                                                                                                                       CTTCGCGGGTAAATCCGCGCCTCCAGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt              >  8:210868‑M1/1‑149 (MQ=255)
                                                                                                                                          gcggggTAAATGCGGCTCCTCCCGGCGGGCACCGGGCGGGGCCC‑‑‑CCGTTCGGTGGGGTTTAACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTg             <  3:14745/145‑1 (MQ=255)
                                                                                                                                          cgcgGGAAAAGCGGCTCCTACCCGCGTGCCCGGTCCGGGGTCC‑‑‑CCGTTAGGCGGGGGTTTCCCGGCGAATAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTgg            <  5:229441/148‑1 (MQ=255)
                                                                                                                                          CGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTg             >  3:296158‑M1/1‑147 (MQ=255)
                                                                                                                                            CGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCaa         <  6:42155‑M1/149‑1 (MQ=255)
                                                                                                                                               taAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt      >  1:135233‑M1/3‑149 (MQ=255)
                                                                                                                                                 AAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGc    <  3:212936‑M1/149‑1 (MQ=255)
                                                                                                                                                   ATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCACAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGccc  <  5:7430‑M1/149‑1 (MQ=255)
                                                                                                                                                   ATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGccc  <  2:4570‑M1/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                         
AATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCC  >  NC_002947/4896251‑4896546

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: