Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 453,367 (A)5→6 coding (179/1476 nt) ampG ← muropeptide:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913453,3621.A96.4% 193.2 / ‑0.6 56coding (184/1476 nt)ampGmuropeptide:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (22/32);  total (23/33)
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

CGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATC  >  NC_000913/453223‑453492
                                                                                                                                            |                                                                                                                                  
cGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAAttt                                                                                                                                     >  1:100555/1‑140 (MQ=255)
  gAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTaa                                                                                                                                  <  1:407636/140‑1 (MQ=255)
   aGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTaaa                                                                                                                                 >  2:928399/1‑140 (MQ=255)
            ccTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAg                                                                                                                         <  2:1010539/140‑1 (MQ=255)
            ccTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGc                                                                                                                        >  2:88783/1‑140 (MQ=255)
             cTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGc                                                                                                                        >  1:1218289/1‑140 (MQ=255)
             cTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGc                                                                                                                        >  1:848744/1‑140 (MQ=255)
              tGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGcc                                                                                                                       <  2:221160/140‑1 (MQ=255)
              tGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGcc                                                                                                                       <  2:1218289/140‑1 (MQ=255)
              tGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGcc                                                                                                                       <  2:707734/140‑1 (MQ=255)
                gTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTg                                                                                                                     <  2:732816/140‑1 (MQ=255)
                gTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTg                                                                                                                     <  2:492982/140‑1 (MQ=255)
                gTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTg                                                                                                                     <  1:828433/140‑1 (MQ=255)
                gTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTg                                                                                                                     <  1:195089/140‑1 (MQ=255)
                  tCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGc                                                                                                                   >  1:1111148/1‑140 (MQ=255)
                     agaAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTa                                                                                                                <  2:1189218/140‑1 (MQ=255)
                          aCCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCaga                                                                                                           >  1:1065490/1‑140 (MQ=255)
                          aCCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCaga                                                                                                           <  2:342417/140‑1 (MQ=255)
                              aTCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAg                                                                                                                         >  2:342086/1‑122 (MQ=255)
                              aTCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAg                                                                                                                         <  1:342086/122‑1 (MQ=255)
                              aTCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAgagaga                                                                                                       >  1:416435/1‑140 (MQ=255)
                               tCGCCGCAATGGCAACTAATAACAGGATTCGCGTGGCGAGCAGCCAACGGCGCCGCCGCGCAAAAAATGGAGGCGTGTAGCGGTACATCAGCGGTGTCCGGAGGAATTTAAAAAACGTAGGCGTCGCCTACCAGAGagaa                                                                                                      <  2:374248/140‑1 (MQ=255)
                                cgccgcAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGagaag                                                                                                     >  1:56290/1‑140 (MQ=255)
                                cgccgcAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGagaag                                                                                                     >  1:706296/1‑140 (MQ=255)
                                cgccgcAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGagaag                                                                                                     >  2:73436/1‑140 (MQ=255)
                                      aaTGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCa                                                                                               <  2:706296/140‑1 (MQ=255)
                                      aaTGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCaa                                                                                              <  1:88783/140‑1 (MQ=255)
                                             aCTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCagagag                                                                                                        >  2:1206484/1‑124 (MQ=255)
                                             aCTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCagagag                                                                                                        <  1:1206484/124‑1 (MQ=255)
                                             aCTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGttt                                                                                        >  2:666767/1‑140 (MQ=255)
                                                aataaCAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTga                                                                                     <  2:1111148/140‑1 (MQ=255)
                                                 ataaCAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTgag                                                                                    <  1:701277/140‑1 (MQ=255)
                                                   aaCAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGAt                                                                                  <  1:73436/140‑1 (MQ=255)
                                                    aCAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTgag                                                                                    <  2:838172/137‑1 (MQ=255)
                                                    aCAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTgag                                                                                    >  1:838172/1‑137 (MQ=255)
                                                        gATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCaa                                                                                              >  1:412849/1‑123 (MQ=255)
                                                        gATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCaa                                                                                              <  2:412849/123‑1 (MQ=255)
                                                         aTTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGatat                                                                            >  2:999329/1‑140 (MQ=255)
                                                           ttGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGAt                                                                                  <  2:1178367/132‑1 (MQ=255)
                                                           ttGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGAt                                                                                  >  1:1178367/1‑132 (MQ=255)
                                                            tGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATtct                                                                         >  1:207191/1‑140 (MQ=255)
                                                                  gCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCg                                                                   <  2:416435/140‑1 (MQ=255)
                                                                  gCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCg                                                                   <  2:207191/140‑1 (MQ=255)
                                                                  gCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCg                                                                   <  1:666767/140‑1 (MQ=255)
                                                                  gCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCg                                                                   <  2:984947/140‑1 (MQ=255)
                                                                  gCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAAAGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCg                                                                   <  2:848744/140‑1 (MQ=255)
                                                                   cGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGActg                                                                   <  2:452349/139‑3 (MQ=255)
                                                                   cGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGActg                                                                   >  1:452349/1‑137 (MQ=255)
                                                                   cGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGt                                                                  >  1:560861/1‑140 (MQ=255)
                                                                     agcagcCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCa                                                                <  1:269292/140‑1 (MQ=255)
                                                                          ccAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAg                                                           <  1:528183/140‑1 (MQ=255)
                                                                            aaCCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGc                                                         <  2:56290/140‑1 (MQ=255)
                                                                                                                 gTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCg                    >  2:827878/1‑140 (MQ=255)
                                                                                                                 gTCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCg                    >  1:1120576/1‑140 (MQ=255)
                                                                                                                  tCCATCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCgg                   <  2:560861/140‑1 (MQ=255)
                                                                                                                      tCAGCGGTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGc               >  1:1072029/1‑140 (MQ=255)
                                                                                                                            gTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCaa             <  2:409449/136‑1 (MQ=255)
                                                                                                                            gTGACCAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCaa             >  1:409449/1‑136 (MQ=255)
                                                                                                                                ccAGAGGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAgg     <  2:1120576/140‑1 (MQ=255)
                                                                                                                                   gagGAATTTAAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATc  <  2:1072029/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                  
CGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATC  >  NC_000913/453223‑453492

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