Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,163,822 T→A I316N (ATC→AAC)  baeS → sensory histidine kinase BaeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,163,8220TA97.2% 233.0 / ‑1.4 71I316N (ATC→AAC) baeSsensory histidine kinase BaeS
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (28/41);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

GCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGATCCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCA  >  NC_000913/2163685‑2163950
                                                                                                                                         |                                                                                                                                 
gcgtAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCTATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAAcc                                                                                                                                 <  1:301888/140‑1 (MQ=255)
gcgtAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAAcc                                                                                                                                 <  1:399371/140‑1 (MQ=255)
 cgtAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAAccc                                                                                                                                <  1:1257050/140‑1 (MQ=255)
  gtAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCa                                                                                                                               <  2:1023458/140‑1 (MQ=255)
  gtAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCa                                                                                                                               <  1:1353249/140‑1 (MQ=255)
   tAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCAc                                                                                                                              >  2:185649/1‑140 (MQ=255)
                 agaCGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTggcg                                                                                                                <  1:317059/140‑1 (MQ=255)
                 agaCGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGCCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTggcg                                                                                                                >  2:160881/1‑140 (MQ=255)
                  gaCGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACCGCTGGAAGTggcgg                                                                                                               <  2:165891/140‑1 (MQ=255)
                  gaCGGTGGCGTCTTTACAGGCGGAGGTCGGGACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTggcgg                                                                                                               <  1:913153/140‑1 (MQ=255)
                       tGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGATCCCACTGCTGGAAGTGGCGggcggc                                                                                                          >  1:1359532/1‑140 (MQ=255)
                          cGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGgcggcggc                                                                                                       >  1:1247101/1‑138 (MQ=255)
                                 aCAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTcgcg                                                                                                <  2:382548/140‑1 (MQ=255)
                                     gCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACg                                                                                           >  1:59856/1‑140 (MQ=255)
                                               gTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGt                                                                                 <  2:315330/140‑1 (MQ=255)
                                               gTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGt                                                                                 >  1:516140/1‑140 (MQ=255)
                                                 acacTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCg                                                                               <  2:1139301/140‑1 (MQ=255)
                                                 acacTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCg                                                                               <  1:611382/140‑1 (MQ=255)
                                                     tGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGGGGGGCATTTCGCGGACGATTCCCCAGGCGTggg                                                                           >  2:718518/1‑139 (MQ=255)
                                                     tGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGGGCATTTCGCGGACGATTCGCCAGTCGTggg                                                                           >  2:1004588/1‑139 (MQ=255)
                                                     tGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCGCTTCGCGAACGATTCGCCAATCGTggg                                                                           >  2:648009/1‑139 (MQ=255)
                                                     tGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGc                                                                           >  2:896798/1‑140 (MQ=255)
                                                      gACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGcc                                                                          >  1:1285789/1‑140 (MQ=255)
                                                      gACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGcc                                                                          >  1:973718/1‑140 (MQ=255)
                                                          aaaCTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGaa                                                                      >  2:1237911/1‑140 (MQ=255)
                                                            aCTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGGAATTCGCGAACGATTTGCCAGTCGCGGCCTGgaaa                                                                    >  2:769576/1‑139 (MQ=255)
                                                               ggTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGGGCATTTTGCGCACGATTCGCCCGTCGTGGCCTGAAAATGc                                                                 >  2:500481/1‑140 (MQ=255)
                                                                 ttGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCGTTTCCCGAACGATTTGCCACTCGTGGCCCGAAAATGGaa                                                               >  2:1182655/1‑140 (MQ=255)
                                                                   gACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAAtt                                                             >  1:727910/1‑140 (MQ=255)
                                                                     cGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAAtttt                                                           <  2:825575/140‑1 (MQ=255)
                                                                       aTCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTcc                                                         <  1:475879/140‑1 (MQ=255)
                                                                         ctcCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCt                                                       >  1:1336600/1‑140 (MQ=255)
                                                                           ccATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGc                                                     <  1:1237911/140‑1 (MQ=255)
                                                                             aTCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCa                                                   <  2:1336600/140‑1 (MQ=255)
                                                                             aTCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCa                                                   <  2:1285789/140‑1 (MQ=255)
                                                                               gaGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTGGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCGAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGa                                                 <  1:674103/139‑1 (MQ=255)
                                                                               cAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGa                                                 <  1:1004588/140‑1 (MQ=255)
                                                                               cAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGa                                                 <  1:769576/140‑1 (MQ=255)
                                                                               cAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGa                                                 <  1:645614/140‑1 (MQ=255)
                                                                               cAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGa                                                 <  1:1182655/140‑1 (MQ=255)
                                                                               cAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGa                                                 <  1:185649/140‑1 (MQ=255)
                                                                               cAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGa                                                 <  1:896798/140‑1 (MQ=255)
                                                                                aGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGAc                                                <  2:347285/140‑1 (MQ=255)
                                                                                  ttGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGGGGGGCCTTTCTCGAACGATTCGCCCGTCCTGGCCTGAAAATGTAATTTTCTCTGCCCGACaa                                              >  2:1342784/1‑139 (MQ=255)
                                                                                  ttGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCCTTTCTCGGACGATTCGCCCGTCGTGGGCTGGAAATGCCATTTTTCCTGCCCGACAg                                              >  2:645614/1‑140 (MQ=255)
                                                                                  ttGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCCCGAACGGATTCCCCCTCCGGGCCTGGAAATGCAAATTTTCCTGCCAGACAg                                              >  2:1337143/1‑140 (MQ=255)
                                                                                   tGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAGCCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGt                                             <  2:656466/140‑1 (MQ=255)
                                                                                   tGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCCGTCGTGGCCTGAAAATGCAATTTTTCCTGCCCGACAGt                                             >  2:1334833/1‑140 (MQ=255)
                                                                                    gTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTa                                            <  1:1337143/140‑1 (MQ=255)
                                                                                      cGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTAtt                                          <  1:340987/140‑1 (MQ=255)
                                                                                          gTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCACTTTGCGAACGATTCGCCCGTCGTGGCCTTGAAATGGAATTTTTCCTTCCCGACAGTATTTccc                                      >  2:340987/1‑139 (MQ=255)
                                                                                             tgatgaAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGGGGCGCGTTTTTCGCACAAATACCCCGGCGGGGGCCGAAAATGCAAATTTTCCTGCCAGACAGTATTACCGTAt                                   >  2:1029585/1‑140 (MQ=255)
                                                                                                    ggCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGAGGCGCAGTTCGCGAAGGAATCCCCCTTAGCGGCCTTAAAGTTCCAATTTTCCTGCCACAGAAGATTTACGTATTTggggg                            >  2:12658/1‑137 (MQ=255)
                                                                                                     gcgcTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGGGCATTTTCCGAACGAATTACCCCGCGGGGCCTGAAAATGCCATTTTTCCTGCCAAACAGTATTACCCCATTTTGCgcc                           >  2:888218/1‑138 (MQ=255)
                                                                                                       gcTCTCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCACTTCGCCAACGATTCGCCCGTCGTGGGCCGAAACTGGAATTTTTCCTGCCCCACAGGATTACCGTATTTTGCGATcg                         >  2:1428473/1‑140 (MQ=255)
                                                                                                        ctctCGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATcgc                        <  1:160881/140‑1 (MQ=255)
                                                                                                          ctcGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGa                      <  2:973718/140‑1 (MQ=255)
                                                                                                          ctcGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGa                      <  1:710587/140‑1 (MQ=255)
                                                                                                          ctcGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGa                      <  1:12658/140‑1 (MQ=255)
                                                                                                          ctcGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCATGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGa                      <  1:505014/140‑1 (MQ=255)
                                                                                                            cGCCTATCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGAcc                    <  1:1029585/140‑1 (MQ=255)
                                                                                                                  tCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTaa              <  2:727910/140‑1 (MQ=255)
                                                                                                                  tCAAAAAGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTaa              <  1:1401339/140‑1 (MQ=255)
                                                                                                                    aaaaaGCTCCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATg            <  2:1224892/140‑1 (MQ=255)
                                                                                                                    aaaaaGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATg            <  2:223451/140‑1 (MQ=255)
                                                                                                                    aaaaaGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATg            <  1:1342784/140‑1 (MQ=255)
                                                                                                                    aaaaaGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATg            <  1:1297606/140‑1 (MQ=255)
                                                                                                                       aaGCACCGGTAGATTTGATCCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAg         <  2:1359532/140‑1 (MQ=255)
                                                                                                                        aGCACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGt        <  1:500481/140‑1 (MQ=255)
                                                                                                                          cACCGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTa      <  2:59856/140‑1 (MQ=255)
                                                                                                                             cGGTAGATTTGAACCCACTGCTGGAAGTGGCGGGCGGGGCCTTTCGCCAAAGAATCGCCAGTCGTGGGCTGAAACAGCAATTTTTCCTGCCCGACAGTAATAACGTATTTTGCGATCGCGACCGTTTTATGCAGTTAtt  >  2:196002/1‑139 (MQ=255)
                                                                                                                             cGGTAGATTTGAACCCACTGCTGGAAGTGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCa  >  2:254245/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                 
GCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATGTCTGATGAAGGCGCTCTCGCCTATCAAAAAGCACCGGTAGATTTGATCCCACTGCTGGAAGTGGCGGGCGGCGCATTTCGCGAACGATTCGCCAGTCGTGGCCTGAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCA  >  NC_000913/2163685‑2163950

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: