Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,917,603 T→G V849V (GTT→GTG barA → sensory histidine kinase BarA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,917,6030TG98.3% 173.4 / ‑5.9 59V849V (GTT→GTGbarAsensory histidine kinase BarA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (28/30);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TGGTGAACCCGAATGCGACCCTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTTGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAGAAGAAGATTT  >  NC_000913/2917464‑2917740
                                                                                                                                           |                                                                                                                                         
tgttaaaCCCGAATGCGACCCTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAACAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGga                                                                                                                                        <  2:100687/138‑1 (MQ=255)
    gcaCCCGTAGGCGACCCTCGCCTGGTATCTCGCACTACGCCAGGCACCCGGAAAAACCTACTTAGCGCGCGATATGCTGCAAATGTTTCTCGATTTCATGCCTTTAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGaaa                                                                                                                                    <  2:378140/141‑1 (MQ=255)
     aaCCCAAATGCGACCCTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGaaaa                                                                                                                                   >  2:1029761/1‑143 (MQ=255)
        ccGAATGCGACCCTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGACTGAAGTTCGCAACAAGGTTGAGGAACATCTGGTGGGAGAAAAccc                                                                                                                                >  2:667729/1‑143 (MQ=255)
         cGAATGCGACCCTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCg                                                                                                                               >  2:33092/1‑143 (MQ=255)
          gAATGCGACCCTCGACTGGCAACTGGCACTGCGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCgg                                                                                                                              <  1:441875/143‑1 (MQ=255)
            gtGCGACCCTCGGCTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGaa                                                                                                                            <  2:175302/142‑1 (MQ=255)
            aTGCGACGCGCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCCGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGaa                                                                                                                            <  1:84993/143‑1 (MQ=255)
                gACCCTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAGCCCGGAAGGcc                                                                                                                        >  1:678905/1‑143 (MQ=255)
                   ccTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTgg                                                                                                                     >  2:525111/1‑143 (MQ=255)
                   ccTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGAGATGCTGCCCACGTTACGCTAGTTCCTGCCTGAAGTTCGCAACAAGGTTGAGGAACAGCTGGTGGGAGAAAGCACGGCAGGCCTgg                                                                                                                     >  2:149170/1‑143 (MQ=255)
                    cTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGt                                                                                                                    >  2:339328/1‑143 (MQ=255)
                        actggcaactggcaCTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGCGGAACAGGTGGTGGGAGAAAACCAGGAAGGCCTGGttgat                                                                                                                >  2:597316/1‑143 (MQ=255)
                           ggcaactggcaCTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGGCTGGttgatgtg                                                                                                             >  1:183151/1‑143 (MQ=255)
                            gcaactggcaCTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGgtgatttga                                                                                                            >  1:1193569/1‑143 (MQ=255)
                                  ggcgCTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATa                                                                                                      <  1:107642/143‑1 (MQ=255)
                                             aGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAAAACCACAGATTATGTTTGGGGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGc                                                                                           >  1:283304/1‑143 (MQ=255)
                                               gcagcaGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAg                                                                                         >  2:897389/1‑143 (MQ=255)
                                               gcagcaGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAg                                                                                         <  1:1137920/143‑1 (MQ=255)
                                                 agcagGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGtt                                                                                       >  1:620966/1‑143 (MQ=255)
                                                   ctggAAAATCCTATTTAGCGCGCGATATGCTGCAAATGTGACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGAACCGGCAGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGc                                                                                     <  1:339441/141‑1 (MQ=255)
                                                   cagGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGc                                                                                     >  2:580915/1‑143 (MQ=255)
                                                     ggAAAAACGCATCTAGCGCGCGGTATGCTGCAAAGGGTACGCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCgg                                                                                   <  1:1179246/143‑1 (MQ=255)
                                                      gagagaCCGATTTGGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGc                                                                                  <  2:636648/138‑1 (MQ=255)
                                                       aaaaaCCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGCACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCt                                                                                 >  1:269387/1‑143 (MQ=255)
                                                       aaaaaCCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATGTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGTTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCt                                                                                 <  1:748902/143‑1 (MQ=255)
                                                        aaaaCCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAGGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCta                                                                                <  2:860583/143‑1 (MQ=255)
                                                          aaCCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATACACTGCATGGCAGTTGCGGCtata                                                                              >  1:246668/1‑143 (MQ=255)
                                                           aCCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGGTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAg                                                                             >  1:1025107/1‑143 (MQ=255)
                                                              gATTTAGCGCGCGATATGCTGCAACTGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCgg                                                                          >  2:372635/1‑143 (MQ=255)
                                                               aTTTAGCGCGCGATATGCGGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGgt                                                                         <  1:737155/143‑1 (MQ=255)
                                                               aTTGAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGgt                                                                         <  1:745187/143‑1 (MQ=255)
                                                                 ttAGCGTGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGgtgt                                                                       <  2:541013/143‑1 (MQ=255)
                                                                   aGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCGTGGCAGTTGCGGCTAGAGAGGTGTGc                                                                     >  2:458839/1‑143 (MQ=255)
                                                                     tgcgcgATATGCTGCAAATGTTACTCGATTGCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCa                                                                   <  2:714182/142‑1 (MQ=255)
                                                                         ctagatGCTGCAAATGTTACGCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGCGGAACAGCTGGTGGGAGAACACCCGGAAGGCCTGGTTGATTTGATTCATAAACTTCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTa                                                               <  2:445549/139‑1 (MQ=255)
                                                                            tatGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATga                                                            <  2:290536/143‑1 (MQ=255)
                                                                            tatGCTGCAAATGTTACTCGATTTCCTGCCTGAAAACCTCAACAAAGTTGAGGAACAGCTGGTGGGAGGAAACCCGGAAGGCCTGGGTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATga                                                            >  2:210960/1‑143 (MQ=255)
                                                                                 tgcAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAAt                                                       >  1:1080240/1‑143 (MQ=255)
                                                                                  gcAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTACGAAGAAtc                                                      >  2:699114/1‑143 (MQ=255)
                                                                                   cAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGCAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAAtct                                                     >  2:1040336/1‑143 (MQ=255)
                                                                                    aaaTGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAAtctc                                                    <  1:1029761/143‑1 (MQ=255)
                                                                                          tACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCaa                                              <  1:33092/143‑1 (MQ=255)
                                                                                          tACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTAGATTCATAAACTGCATGGAAGTTGCGGCTATAGAGGTGTACAACGAATGAAGCACCTCGGCCaa                                              >  2:1001235/1‑143 (MQ=255)
                                                                                            cTCGATTTCCTGCCAGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACt                                            <  1:251545/143‑1 (MQ=255)
                                                                                                 tttCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCg                                       <  1:848288/143‑1 (MQ=255)
                                                                                                            aaGTTCGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTa                            <  1:922717/143‑1 (MQ=255)
                                                                                                                 cGCAACAAAGTTGAGGAACAGCTGGTGGTAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAg                       <  2:126280/143‑1 (MQ=255)
                                                                                                                 cGCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACCACAGCTACGtcg                       >  1:1196356/1‑141 (MQ=255)
                                                                                                                  gCAACAAAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACATCTACGTAGc                      >  2:151790/1‑143 (MQ=255)
                                                                                                                    aacaaAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCCCGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCgg                    <  1:580915/143‑1 (MQ=255)
                                                                                                                      caaAGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTa                  >  1:665604/1‑143 (MQ=255)
                                                                                                                       aaaGTTGAGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTAc                 <  1:923937/143‑1 (MQ=255)
                                                                                                                       aaaGGTGAGGAGCAGCGGGGGGGAGAAAGCCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGGTGGGGCTATCGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTAc                 <  1:424499/143‑1 (MQ=255)
                                                                                                                            tGAGGAATATGTGGTGGCAGAAAACCCGGAAGGCTTTGTTGATGTGATTCATAAACTGCATGGCAGTTGCGGCGATAGCGGTGGGCCACGCATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaag            <  2:129011/143‑1 (MQ=255)
                                                                                                                              aGGAACAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCGACTGATCGAACAACAGCTGCGTAGCGGAAcagaagga          >  1:34341/1‑136 (MQ=255)
                                                                                                                                   cAGCTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAagaagaaga     >  1:705771/1‑143 (MQ=255)
                                                                                                                                     gCTGGTGGGAGAAAGGCGGCGAGGCCTGGTTGATTTGATTCATAAACTGCGTGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAGAAGAAGAtt   <  2:602085/143‑1 (MQ=255)
                                                                                                                                      cTGGTGGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAGAAGAAGAttt  <  2:620966/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TGGTGAACCCGAATGCGACCCTCGACTGGCAACTGGCACTACGCCAGGCAGCAGGAAAAACCGATTTAGCGCGCGATATGCTGCAAATGTTACTCGATTTCCTGCCTGAAGTTCGCAACAAAGTTGAGGAACAGCTGGTTGGAGAAAACCCGGAAGGCCTGGTTGATTTGATTCATAAACTGCATGGCAGTTGCGGCTATAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAGAAGAAGATTT  >  NC_000913/2917464‑2917740

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: