Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,842,455 C→T *445* (TAG→TAA adeQ ← adenine permease, high affinity, adenine:H+ symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,842,4550CT98.7% 239.4 / NA 75*445* (TAG→TAAadeQadenine permease, high affinity, adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (40/34);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGT  >  NC_000913/3842320‑3842586
                                                                                                                                       |                                                                                                                                   
aTCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTaa                                                                                                                                  <  2:1457903/139‑1 (MQ=255)
       ttCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATccacca                                                                                                                           <  1:1389106/139‑1 (MQ=255)
       ttCTGGCCACGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATccacca                                                                                                                           <  1:452710/139‑1 (MQ=255)
                cGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCt                                                                                                                  <  1:580094/139‑1 (MQ=255)
                 ggATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCtt                                                                                                                 <  1:1238362/139‑1 (MQ=255)
                  gATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTc                                                                                                                <  1:81118/139‑1 (MQ=255)
                       cAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGc                                                                                                           >  1:687863/1‑139 (MQ=255)
                            cccGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACa                                                                                                      >  1:1075786/1‑139 (MQ=255)
                                 aaaCCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCt                                                                                                 <  2:1075786/139‑1 (MQ=255)
                                               cacTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACcaccaccac                                                                                   >  2:222404/1‑139 (MQ=255)
                                                acTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcaca                                                                                  <  2:192596/139‑1 (MQ=255)
                                                  tACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                <  1:822565/139‑1 (MQ=255)
                                                  tACCGCCAGCGTGGCGAGCAGTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATTATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                <  1:1227390/139‑1 (MQ=255)
                                                   aCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGaccacca                                                                                       >  1:1412519/1‑131 (MQ=255)
                                                   aCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGaccacca                                                                                       <  2:1412519/131‑1 (MQ=255)
                                                   aCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacacac                                                                               >  1:1354006/1‑139 (MQ=255)
                                                           cGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTaa                                                                                                                                  >  2:396793/1‑80 (MQ=255)
                                                           cGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTaa                                                                                                                                  <  1:396793/80‑1 (MQ=255)
                                                           cGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAAGCAGAACTGATGATATCTGGAAACACTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCAACCACCACAACACACAGGTTCAg                                                                       >  2:441233/1‑139 (MQ=255)
                                                              ggCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATc                                                                    >  1:1371511/1‑139 (MQ=255)
                                                               gCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATccac                                                                                                                             <  2:452690/81‑1 (MQ=255)
                                                               gCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATccac                                                                                                                             >  1:452690/1‑81 (MQ=255)
                                                                    cATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCa                                                              <  1:490412/139‑1 (MQ=255)
                                                                         tCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGcc                                                         <  1:222404/139‑1 (MQ=255)
                                                                              aTCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGaata                                                                                                                      <  2:467153/73‑1 (MQ=255)
                                                                              aTCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGaata                                                                                                                      >  1:467153/1‑73 (MQ=255)
                                                                              aTCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGa                                                                                                                         <  1:808407/70‑1 (MQ=255)
                                                                              aTCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGa                                                                                                                         >  2:808407/1‑70 (MQ=255)
                                                                                  tGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATa                                                >  2:728408/1‑139 (MQ=255)
                                                                                   gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGac                                                                                            <  2:1365145/94‑1 (MQ=255)
                                                                                   gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGac                                                                                            >  1:1365145/1‑94 (MQ=255)
                                                                                   gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATAc                                               >  1:1027553/1‑139 (MQ=255)
                                                                                   gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATAc                                               >  1:124845/1‑139 (MQ=255)
                                                                                     gTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGaata                                                                                                                      <  2:386714/66‑1 (MQ=255)
                                                                                     gTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGaata                                                                                                                      >  1:386714/1‑66 (MQ=255)
                                                                                          ccTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCatg                                        >  1:1397252/1‑139 (MQ=255)
                                                                                           cTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCg                                                       >  2:473554/1‑123 (MQ=255)
                                                                                           cTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCg                                                       <  1:473554/123‑1 (MQ=255)
                                                                                              cGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATccac                                                                                                                             <  1:792982/50‑1 (MQ=255)
                                                                                              cGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATccac                                                                                                                             >  2:792982/1‑50 (MQ=255)
                                                                                                tAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCa                                          >  1:116902/1‑131 (MQ=255)
                                                                                                tAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCa                                          <  2:116902/131‑1 (MQ=255)
                                                                                                tAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAg                                  >  1:1208510/1‑139 (MQ=255)
                                                                                                   tCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGcc                                                         >  1:623234/1‑113 (MQ=255)
                                                                                                   tCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGcc                                                         <  2:623234/113‑1 (MQ=255)
                                                                                                   tCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTAc                               >  2:986320/1‑139 (MQ=255)
                                                                                                        aCTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCg                                                                                              <  1:1086382/71‑1 (MQ=255)
                                                                                                        aCTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCg                                                                                              >  2:1086382/1‑71 (MQ=255)
                                                                                                        aCTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                >  1:1259408/1‑85 (MQ=255)
                                                                                                        aCTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                <  2:1259408/85‑1 (MQ=255)
                                                                                                        aCTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCGCACACAGGTTCAGATCGCGCCAgcgc                                                          <  1:536293/107‑1 (MQ=255)
                                                                                                        aCTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAgcgc                                                          >  2:536293/1‑107 (MQ=255)
                                                                                                           cATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATaaa                       <  2:1354006/139‑1 (MQ=255)
                                                                                                            aTGATATCTGGAAACCCTCGCCAGTTATTAATccacca                                                                                                                           <  2:1325564/38‑1 (MQ=255)
                                                                                                            aTGATATCTGGAAACCCTCGCCAGTTATTAATccacca                                                                                                                           >  1:1325564/1‑38 (MQ=255)
                                                                                                            aTGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCGAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAg                      >  1:882018/1‑139 (MQ=255)
                                                                                                            aTGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAg                      >  2:1380099/1‑139 (MQ=255)
                                                                                                             tGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGc                     >  2:1219429/1‑139 (MQ=255)
                                                                                                                 atCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAg                 >  1:62684/1‑139 (MQ=255)
                                                                                                                    tGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGc              <  2:173326/139‑1 (MQ=255)
                                                                                                                     ggAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCa             >  1:928784/1‑139 (MQ=255)
                                                                                                                      gAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCaa            >  1:190084/1‑139 (MQ=255)
                                                                                                                        aaCCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATc          >  2:564200/1‑139 (MQ=255)
                                                                                                                          cccTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                <  2:386370/67‑1 (MQ=255)
                                                                                                                          cccTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                >  1:386370/1‑67 (MQ=255)
                                                                                                                           ccTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCa                                                              <  2:441713/84‑1 (MQ=255)
                                                                                                                           ccTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCa                                                              >  1:441713/1‑84 (MQ=255)
                                                                                                                             tCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTc     >  2:571572/1‑139 (MQ=255)
                                                                                                                              cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGaccacc                                                                                        >  1:752247/1‑55 (MQ=255)
                                                                                                                              cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGaccacc                                                                                        <  2:752247/55‑1 (MQ=255)
                                                                                                                              cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGacc                                                                                           <  1:895940/52‑1 (MQ=255)
                                                                                                                              cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGacc                                                                                           >  2:895940/1‑52 (MQ=255)
                                                                                                                              cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGc              <  2:1249290/129‑1 (MQ=255)
                                                                                                                              cGCCAATTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGc              >  1:1249290/1‑129 (MQ=255)
                                                                                                                                ccAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGt  >  1:590531/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                   
ATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGT  >  NC_000913/3842320‑3842586

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