Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,210,811 C→T pseudogene (1712/2001 nt) yehQ → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,210,8110CT93.8% 255.4 / 5.8 81pseudogene (1712/2001 nt)yehQpseudogene
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (38/38);  total (41/40)
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.74e-01

GTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTG  >  NC_000913/2210714‑2210907
                                                                                                 |                                                                                                
gTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGc                                                                                                <  1:279901/100‑1 (MQ=255)
gTCAGTGTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGc                                                                                                <  2:243137/100‑1 (MQ=255)
 tCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCt                                                                                               <  2:159602/100‑1 (MQ=255)
   aGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTAc                                                                                             <  1:481149/100‑1 (MQ=255)
    gTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTAcc                                                                                            >  2:1399658/1‑100 (MQ=255)
            gAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCa                                                                                    >  2:167106/1‑100 (MQ=255)
             aaGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAg                                                                                   >  2:146666/1‑100 (MQ=255)
             aaGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAg                                                                                   >  2:359954/1‑100 (MQ=255)
               gggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCg                                                                                 <  1:559872/100‑1 (MQ=255)
                ggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGt                                                                                >  2:1533908/1‑100 (MQ=255)
                  cAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGttt                                                                              >  1:473318/1‑100 (MQ=255)
                     ttGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGc                                                                           <  1:1322598/100‑1 (MQ=255)
                     ttGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGc                                                                           <  1:747781/100‑1 (MQ=255)
                      tGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGcc                                                                          >  2:647460/1‑100 (MQ=255)
                       gACACTGTAGCCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGcc                                                                          <  2:775316/99‑1 (MQ=255)
                       gACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCt                                                                         >  2:1267790/1‑100 (MQ=255)
                       gACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCt                                                                         >  1:1177654/1‑100 (MQ=255)
                        acacTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtg                                                                        <  2:1268242/100‑1 (MQ=255)
                        acacTGGAACCACGGGTGATCATTACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGCCACTGCTACCTTCTGCCAGCGTTTTGCCtg                                                                        <  2:712656/100‑1 (MQ=255)
                         cacTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgt                                                                       >  2:221603/1‑100 (MQ=255)
                           cTGGAACCATGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTTGCTCTGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                      <  2:1030253/99‑1 (MQ=255)
                           cTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGc                                                                     <  1:1273549/100‑1 (MQ=255)
                             ggAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCg                                                                   >  2:112810/1‑100 (MQ=255)
                                aCCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCt                                                                <  2:125425/100‑1 (MQ=255)
                                 ccACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTa                                                               >  2:1516780/1‑100 (MQ=255)
                                  cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                              >  1:658209/1‑100 (MQ=255)
                                  cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                              >  1:1448596/1‑100 (MQ=255)
                                  cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                              >  1:797440/1‑100 (MQ=255)
                                    cGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTg                                                            >  1:1369879/1‑100 (MQ=255)
                                    cGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTg                                                            >  2:938896/1‑100 (MQ=255)
                                     gggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGc                                                           <  1:651825/100‑1 (MQ=255)
                                      ggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCt                                                          <  1:146666/100‑1 (MQ=255)
                                       gTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtc                                                         >  2:1363716/1‑100 (MQ=255)
                                       gTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtc                                                         <  1:494188/100‑1 (MQ=255)
                                        tgatgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtca                                                        <  2:1479445/100‑1 (MQ=255)
                                        tgatgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtca                                                        <  1:919437/100‑1 (MQ=255)
                                        tgatgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtca                                                        >  1:1101239/1‑100 (MQ=255)
                                          atgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcatc                                                      <  1:325764/100‑1 (MQ=255)
                                           tgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcatca                                                     >  1:74332/1‑100 (MQ=255)
                                               gaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTg                                                 >  1:1196153/1‑100 (MQ=255)
                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                                <  1:88378/100‑1 (MQ=255)
                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                                >  2:267680/1‑100 (MQ=255)
                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               >  2:570105/1‑100 (MQ=255)
                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               >  1:570753/1‑100 (MQ=255)
                                                  aaaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTg                                              <  2:579652/100‑1 (MQ=255)
                                                     aaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATa                                           >  1:1223165/1‑100 (MQ=255)
                                                        ccGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGc                                        >  2:1231152/1‑100 (MQ=255)
                                                           tGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGc                                     <  2:847592/100‑1 (MQ=255)
                                                           tGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGc                                     <  1:1533918/100‑1 (MQ=255)
                                                            gCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGcc                                    >  1:423413/1‑100 (MQ=255)
                                                             ccTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCa                                   >  1:1136119/1‑100 (MQ=255)
                                                               tGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAgg                                 >  1:582624/1‑100 (MQ=255)
                                                                gggCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGc                                >  2:751615/1‑100 (MQ=255)
                                                                gggCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGc                                >  1:837778/1‑100 (MQ=255)
                                                                 ggCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCg                               >  1:737979/1‑100 (MQ=255)
                                                                  gcgcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGt                              <  1:647463/100‑1 (MQ=255)
                                                                   cgcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGtt                             <  1:486629/100‑1 (MQ=255)
                                                                    gcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTa                            >  1:1564960/1‑100 (MQ=255)
                                                                     cTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTAc                           <  1:1140687/100‑1 (MQ=255)
                                                                         aCGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGGTACGCTGCCAGGCGTTACTTAt                       <  2:914317/100‑1 (MQ=255)
                                                                          cGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTAtt                      >  1:233903/1‑100 (MQ=255)
                                                                           gCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCGGATACGCTGCCAGGCGTTACTTATTc                     <  1:960712/100‑1 (MQ=255)
                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                    >  1:1558897/1‑100 (MQ=255)
                                                                              gAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAAc                  <  1:1013263/100‑1 (MQ=255)
                                                                              gAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAAc                  <  2:1187263/100‑1 (MQ=255)
                                                                               aaaCTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACt                 <  2:262364/100‑1 (MQ=255)
                                                                               aaaCTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACt                 <  1:1267799/100‑1 (MQ=255)
                                                                                   tACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTc             <  1:319804/100‑1 (MQ=255)
                                                                                   tACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTc             <  2:464510/100‑1 (MQ=255)
                                                                                   tACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTc             <  2:1136112/100‑1 (MQ=255)
                                                                                    aCGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTcc            <  2:1223158/100‑1 (MQ=255)
                                                                                     cGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCa           <  2:750319/100‑1 (MQ=255)
                                                                                     cGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCa           >  1:485997/1‑100 (MQ=255)
                                                                                        cTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATaa        <  2:1369870/100‑1 (MQ=255)
                                                                                         tgtgGCATTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAAc       >  1:1398961/1‑100 (MQ=255)
                                                                                          gtgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACg      >  2:1532119/1‑100 (MQ=255)
                                                                                           tgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACgg     <  2:1196146/100‑1 (MQ=255)
                                                                                            ggCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGc    >  2:88361/1‑100 (MQ=255)
                                                                                            ggCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGc    >  2:88480/1‑100 (MQ=255)
                                                                                             gCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCt   <  2:148233/100‑1 (MQ=255)
                                                                                              cACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTg  >  2:253523/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTG  >  NC_000913/2210714‑2210907

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