Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,478,3240TG43.9% 91.1 / 15.3 57G33G (GGT→GGGynbBputative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (6/26);  new base G (25/0);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 7.52e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTT  >  NC_000913/1478175‑1478473
                                                                                                                                                       |                                                                                                                                                     
gACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTg                                                                                                                                                      <  1:347407/151‑1 (MQ=255)
 aCTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGa                                                                                                                                                     <  2:639051/151‑1 (MQ=255)
     cGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGtt                                                                                                                                                 <  1:361288/151‑1 (MQ=255)
              gTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGt                                                                                                                                        <  1:785255/151‑1 (MQ=255)
                tGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAAAGGGGGGGAAAAGGGGGGAGGGTTCACTACGtct                                                                                                                                      >  1:139054/1‑149 (MQ=255)
                 gTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTAtt                                                                                                                                     <  1:872428/151‑1 (MQ=255)
                  tCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGGGAGAAAAAGGGGGGGGGGTACATTTCGTAtt                                                                                                                                     >  2:932279/1‑150 (MQ=255)
                     gAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGGGAGAAAAAGGGGGGAGGGTACATTTCGTATTCg                                                                                                                                   >  1:330438/1‑149 (MQ=255)
                     gAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTa                                                                                                                                 <  1:705519/151‑1 (MQ=255)
                       gtgtGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTAcc                                                                                                                               <  2:915179/151‑1 (MQ=255)
                       gtgtGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTAcc                                                                                                                               <  2:201290/151‑1 (MQ=255)
                           gCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCTGGGGGAAAAGGGGGGGGGGGTTTCTTTAGTATTCGTACCt                                                                                                                              >  1:185467/1‑148 (MQ=255)
                              tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggtgg                                                                                                                         <  2:887477/150‑1 (MQ=255)
                              tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggtg                                                                                                                          <  1:408309/149‑1 (MQ=255)
                                 tATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACCGGGAGAAAAGGGGGGGGGGGTTCATTTCGTATTTTTACCTTGTggggg                                                                                                                      >  1:818585/1‑149 (MQ=255)
                                  aTGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCAGCGGGGAAAAGGGGTGGGGTTTTAATATCTATTTGTACCtgggtggg                                                                                                                       >  1:116348/1‑143 (MQ=255)
                                     aGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGt                                                                                                                 <  2:677257/151‑1 (MQ=255)
                                             aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAAGGGGGGGGGGTATATTTCGTATTCGTACCTGGTGGGGCGTCGTCATTTGtt                                                                                                         >  1:494707/1‑151 (MQ=255)
                                             aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGtt                                                                                                         <  2:835393/151‑1 (MQ=255)
                                             aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGAGAGGAAAAGGGGGGGGGGTACATTTCGTATTTGTACCTTGGGGGGGGTCCTCAtttttt                                                                                                         >  1:658040/1‑151 (MQ=255)
                                               tataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGAGAGAAAAGGGGGGGGGGGTTAATTTAGTTTTTGTAACTTGGGGGGGATCGGCATTTGtttt                                                                                                       >  1:109592/1‑151 (MQ=255)
                                                ataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACCGGGAGAAAAGGGGGGGGGGGTAAACTTAACATTTGTATCTTGTGGGGGATTGTCCTTTGttttt                                                                                                      >  2:768944/1‑151 (MQ=255)
                                                  aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAAGGGGGGGGGGTACATTACGTATTCGTACCTTGGGGGGGATCGTCATtttttttt                                                                                                      >  2:160607/1‑149 (MQ=255)
                                                    tGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAAGGGGGGGGGGTTCATTTAGTATTTGTACCTGGTGGGGGATCGTCCTttttttttt                                                                                                     >  2:395772/1‑147 (MQ=255)
                                                    tGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGAGGGTACATTTCGTATTTGTACCTGGGGGGGGATCGTCCTTTTTTTTTTaa                                                                                                   >  2:409832/1‑150 (MQ=255)
                                                      cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGGGAAAAAGGGGGGAGGGTTAATTACGTTTTCGTACCTTGTGGGGGGTTGTCATttttttttt                                                                                                     >  1:380636/1‑145 (MQ=255)
                                                          aaaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGATTACGGGGGGAAAAAAGGGGGAGGATAACATAAATAATTGTTCCTGGGGGGGGATCGGTATTTGttctttcttt                                                                                                 >  1:1613/1‑141 (MQ=255)
                                                                ctctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAAGGGGGGGGGGTACATTACGTATTCGTACCCGGTGGGGGATCGTCATTTTTTTTTTTATTGGGAtttttgg                                                                                      >  2:6603/1‑147 (MQ=255)
                                                                  ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAgg                                                                                      <  1:641299/149‑1 (MQ=255)
                                                                       aaCACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAAGGGGGGGGGGTACATTACGTATTTGTACCTGGGGGGGCATCGTCATTTGTTTTTTAATGGTGGTTTTAGGTCCACGc                                                                               >  2:653170/1‑151 (MQ=255)
                                                                             tttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGa                                                                           <  2:250103/149‑1 (MQ=255)
                                                                             tttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGa                                                                           <  1:545545/149‑1 (MQ=255)
                                                                              ttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTTCTGTGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTAACCGGGGGGGCGTCGGCATTTGTTTTTTAATGGGGGTTTTTGGTCCACGCTgggt                                                                          >  2:527426/1‑147 (MQ=255)
                                                                                cgcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgac                                                                       <  1:953665/150‑1 (MQ=255)
                                                                                 gcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgt                                                                     <  2:669211/151‑1 (MQ=255)
                                                                                 gcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGAGGGTACATTACGTATTCGTACCTGGGGGGGGATCGTCCTTTTTTTTTTAATGGTGATTTTAGGTCCACGCTGGGtgacgt                                                                     >  2:623281/1‑151 (MQ=255)
                                                                                  cgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgtt                                                                    <  2:330442/151‑1 (MQ=255)
                                                                                   gttgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAAGGGGGGAGGGTAAATTTCCTATTCGGTCCTGGGGGGGCGTCGTCATTTGTTTTTTTATGGTGGTTTTTGGTCCACGCTGGAtgacgtt                                                                    >  1:319861/1‑150 (MQ=255)
                                                                                      gttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttgacg                                                                <  1:34999/151‑1 (MQ=255)
                                                                                       ttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAGGGGGGGAGGGTAAATTACGGATTTGGACCTGGTGGGGGCTCGGCATTTTTTTTTCAATGGGGATTTTAGGGCCACCCTGGGtgacgttttcg                                                                >  2:796031/1‑146 (MQ=255)
                                                                                         aaTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATgacgttgacgtt                                                              <  1:108963/150‑1 (MQ=255)
                                                                                          aTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGTGAAAAAGGGGGGGGGTTACAATACGTATTTGTACCTGGTGGGGGCTCCGCATTTTTTTTTTTATGGTGGTTTTAGGGCCCCGCTTGATTACGTTGGCGgttt                                                            >  2:776242/1‑151 (MQ=255)
                                                                                           ttttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAGGGGGGGAGGGTACATTTTGTATTTGTAACCGGTGGGGGATTGTCCTTTGTTTTTTAATGGTGGTTTTTGGTCCCCGCTGGATGACGTTGAAGttttt                                                           >  1:322811/1‑151 (MQ=255)
                                                                                                   aCGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGGAAAGGGGGGGGGGGTATATTACGTATTCGTACCTGGGGGGGGATCGGCATTTTTTTTTTAATGGGGGTTTTTGGGCCACGCTGGGTGACGTTGACGtttttt                                                          >  1:175836/1‑144 (MQ=255)
                                                                                                         aTAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGAGGGTACAATACTTTTTTGTTCCTTGGGGGGGATCGTCCGTTTTTTTTTTATTGGGGTTTTAGGGCCACGCTGGGTGAAGTTTACGTTTTTTTTGCTggt                                                  >  1:411672/1‑144 (MQ=255)
                                                                                                              ccGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTg                                        <  1:726512/151‑1 (MQ=255)
                                                                                                               cGTTTTCTGCTATGGCGGTTTACCGAGGGAAAAGGGGGGGAGGGTAAATTACGTATTCGTACCCGGTGGGGGCTCGTCATTTTTTTTTTTATGGGGGTTTTAGGGCCCCGCTGGATGACGTTGTCGTTTTTTTCGCTGATAAGtttt                                           >  1:830387/1‑147 (MQ=255)
                                                                                                                            ggCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAat                           <  2:139054/150‑1 (MQ=255)
                                                                                                                                    aCCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACg                   <  2:844748/151‑1 (MQ=255)
                                                                                                                                         agagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtat              <  2:185468/151‑1 (MQ=255)
                                                                                                                                          gagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtata             <  2:700970/151‑1 (MQ=255)
                                                                                                                                          gagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtat              >  1:224484/1‑150 (MQ=255)
                                                                                                                                             aaaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCt          >  1:218114/1‑151 (MQ=255)
                                                                                                                                              aaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTg         >  2:512666/1‑151 (MQ=255)
                                                                                                                                                aGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGt        >  2:736523/1‑150 (MQ=255)
                                                                                                                                                    gggTGAGGGTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAt    >  1:903070/1‑151 (MQ=255)
                                                                                                                                                      gTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAttt  <  1:532185/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTT  >  NC_000913/1478175‑1478473

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