Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,893 C→A 8.7% intergenic (‑199/‑53) tktA ← / → sahR transketolase A/methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,657,8930CA 8.7% 251.5 / 10.3 80intergenic (‑199/‑53)tktA/sahRtransketolase A/methionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (37/36);  new base A (4/3);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

ATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGC  >  NC_002947/5657745‑5658039
                                                                                                                                                    |                                                                                                                                                      
aTTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAAc                                                                                                                                                        >  2:375502/1‑149 (MQ=255)
     cGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACtttt                                                                                                                                                    <  2:56617/148‑1 (MQ=255)
      gCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                                   <  2:144058/148‑1 (MQ=255)
      gCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                                   >  5:46514/1‑148 (MQ=255)
       cccACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                                   <  5:212290/147‑1 (MQ=255)
       cccACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                                   >  6:212290/1‑147 (MQ=255)
              ggggCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTTGATATTgg                                                                                                                                           <  1:373637/148‑1 (MQ=255)
                   tcTTGGGGAAAACATTTAAGCGCCTAGCGCCTGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGTGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTTCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGc                                                                                                                                     <  8:203987/148‑1 (MQ=255)
                    cTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGc                                                                                                                                     >  5:412886/1‑148 (MQ=255)
                          gCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTCGATATTgg                                                                                                                                           >  7:354042/1‑136 (MQ=255)
                          gCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTCGATATTgg                                                                                                                                           <  8:354042/136‑1 (MQ=255)
                           cAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCAt                                                                                                                             >  4:9439/1‑149 (MQ=255)
                           cAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCAt                                                                                                                             <  6:46514/149‑1 (MQ=255)
                            aaaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACtttt                                                                                                                                                    >  7:179604/1‑125 (MQ=255)
                            aaaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACtttt                                                                                                                                                    <  8:179604/125‑1 (MQ=255)
                            aaaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATAACAAAACTTTTTGATATTGGCCTTGCTAGGGCATc                                                                                                                            >  7:64887/1‑149 (MQ=255)
                               cATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGAt                                                                                                                         <  1:28594/149‑1 (MQ=255)
                                aTTTAAGCGCCAAGCGCCTGCTTTAAGCCGCAACTGAACAAGCAGAGCGAGACTCGCGGTGCCCCAAAGCGTGTAACATGAGGCTTTCGGCTTGCAGCTCATCGCCAATAACAAAACATATTGATATTGGCCTTGCTAGGGCATCGAta                                                                                                                        >  7:90011/1‑148 (MQ=255)
                                 tttAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGc                                                                                                                       <  6:412886/149‑1 (MQ=255)
                                 tttAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGc                                                                                                                       <  5:15886/149‑1 (MQ=255)
                                    aagcgccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCt                                                                                                                     <  2:218193/148‑1 (MQ=255)
                                    aagcgccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTg                                                                                                                    >  5:66449/1‑149 (MQ=255)
                                        gccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTc                                                                                                                  >  5:167269/1‑147 (MQ=255)
                                        gccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTa                                                                                                               <  8:180968/149‑1 (MQ=255)
                                               gccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTgcc                                                                                                       >  5:26090/1‑149 (MQ=255)
                                                     tttAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCa                                                                                               >  2:263812/1‑149 (MQ=255)
                                                       tAAGCAGCAACTGAAGAAGCAGAGCGCGACTCGGTGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAATCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCa                                                                                               <  4:223086/147‑1 (MQ=255)
                                                           ccGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTc                                                                                                                  >  8:361852/1‑128 (MQ=255)
                                                           ccGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTc                                                                                                                  <  7:361852/128‑1 (MQ=255)
                                                           ccGCAACTGAAGAAGCAGACCGCGACTCGGGCCGACCCAAAGCGAGTAACTTGAGCCTTTCGGCTTGCAGCTCATCGCCAAAATCAAAACTTTTAGAAAATTGCCTACCAAGGGAATCCATGCCTGTCTAGACTGCCC‑‑CCCCATGCAt                                                                                          >  4:285654/1‑148 (MQ=255)
                                                                 cTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCtgcgtg                                                                                   >  4:144232/1‑149 (MQ=255)
                                                                  tGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCtgcgtgc                                                                                  <  6:66449/149‑1 (MQ=255)
                                                                    aagaagCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCtgcgtgc                                                                                  <  1:263812/147‑1 (MQ=255)
                                                                      gaagCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCaa                                                                              <  6:167269/149‑1 (MQ=255)
                                                                            gagCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATc                                                                        >  4:408776/1‑149 (MQ=255)
                                                                            gagCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTATTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAAACTGCGTTCGCAATCGAt                                                                         >  8:320338/1‑148 (MQ=255)
                                                                               cgcgACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATcccc                                                                     >  2:338362/1‑149 (MQ=255)
                                                                                 cgACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGa                                                                   <  8:149947/149‑1 (MQ=255)
                                                                                 cgACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGa                                                                   <  3:220480/149‑1 (MQ=255)
                                                                                   aCTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAg                                                                  <  4:197062/148‑1 (MQ=255)
                                                                                   aCTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAg                                                                  >  3:197062/1‑148 (MQ=255)
                                                                                    cTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGc                                                                 >  3:144731/1‑148 (MQ=255)
                                                                                       ggggCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGCTGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAcg                                                             >  5:207325/1‑149 (MQ=255)
                                                                                         ggCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAg                                                                  >  4:373335/1‑142 (MQ=255)
                                                                                         ggCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAg                                                                  <  3:373335/142‑1 (MQ=255)
                                                                                            gcCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGaaa                                                        <  3:144232/149‑1 (MQ=255)
                                                                                              cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATg                                                                                             <  6:265587/110‑1 (MQ=255)
                                                                                              cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCt                                                                                        >  1:246659/1‑115 (MQ=255)
                                                                                              cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCt                                                                                        <  2:246659/115‑1 (MQ=255)
                                                                                              cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTGTTTGATATTTGCCTTTCTATGGCCTAGATTCCTTTCTAGCCTTGCC‑‑CCCCATGAATATTCGTTCGCAATAGATACCCCAGCAACGCGAAACa                                                      >  1:318307/1‑149 (MQ=255)
                                                                                                       tgtAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGcc                                             >  7:247166/1‑149 (MQ=255)
                                                                                                         tAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCt                                           <  4:144731/149‑1 (MQ=255)
                                                                                                            cTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgtg                                        <  3:408776/149‑1 (MQ=255)
                                                                                                               gAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCa                                      >  4:233635/1‑148 (MQ=255)
                                                                                                                aGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCATCCGGGGCACGGACTCGATNAAACTGCGTGCGCAATAGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgtg                                        >  1:337946/1‑149 (MQ=255)
                                                                                                                       cGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGt                              <  6:26090/148‑1 (MQ=255)
                                                                                                                              cAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgc                      <  3:233635/149‑1 (MQ=255)
                                                                                                                                gCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgc                    >  3:327264/1‑149 (MQ=255)
                                                                                                                                 cTCATCTCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgct                   <  4:2070/149‑1 (MQ=255)
                                                                                                                                    aTCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgctgcg                >  3:156449/1‑149 (MQ=255)
                                                                                                                                      ggCCAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGcc              >  3:171007/2‑149 (MQ=255)
                                                                                                                                       gCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgt                                         <  6:250215/121‑1 (MQ=255)
                                                                                                                                       gCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgt                                         >  5:250215/1‑121 (MQ=255)
                                                                                                                                       gCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCt             <  1:338362/149‑1 (MQ=255)
                                                                                                                                        ccAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCt             >  2:13350/1‑148 (MQ=255)
                                                                                                                                        ccAATATCAAAAAGTTTTGATATTGGCCATGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCAGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGc                                  >  1:232788/1‑127 (MQ=255)
                                                                                                                                        acAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATTCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGc                                  <  2:232788/126‑1 (MQ=255)
                                                                                                                                          aaTATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGcc              <  6:309325/145‑1 (MQ=255)
                                                                                                                                          aaTATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGcc              >  5:309325/1‑145 (MQ=255)
                                                                                                                                           atatCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTgg                                                  <  2:243103/108‑1 (MQ=255)
                                                                                                                                           atatCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTgg                                                  >  1:243103/1‑108 (MQ=255)
                                                                                                                                            tatCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTg            >  2:30345/1‑145 (MQ=255)
                                                                                                                                            tatCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTg            <  1:30345/145‑1 (MQ=255)
                                                                                                                                              tCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTa      >  2:216384/1‑149 (MQ=255)
                                                                                                                                               cAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTg                             <  2:56949/125‑1 (MQ=255)
                                                                                                                                               cAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTg                             >  1:56949/1‑125 (MQ=255)
                                                                                                                                                aaaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGTACGTAtt    <  8:64887/149‑1 (MQ=255)
                                                                                                                                                aaaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTAtt    >  6:346351/1‑149 (MQ=255)
                                                                                                                                                 aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTg   >  2:361333/1‑149 (MQ=255)
                                                                                                                                                 aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGACTGAACGTATTg   >  7:182214/1‑149 (MQ=255)
                                                                                                                                                  aaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTgc  <  1:216384/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGC  >  NC_002947/5657745‑5658039

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: