Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,5960GA42.9% 50.1 / 36.2 56F13F (TTC→TTTymgDuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (26/6);  new base A (0/24);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.54e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCAC  >  NC_000913/1222447‑1222744
                                                                                                                                                        |                                                                                                                                                    
gCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCg                                                                                                                                                      >  2:575292/1‑150 (MQ=255)
  aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGaaa                                                                                                                                                   <  2:518552/151‑1 (MQ=255)
  aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGaa                                                                                                                                                    >  2:1149127/1‑150 (MQ=255)
                 cattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGcc                                                                                                                                    >  1:1241490/1‑151 (MQ=255)
                              aTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAtt                                                                                                                       >  1:545067/1‑151 (MQ=255)
                                 tCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAttttt                                                                                                                    >  2:147717/1‑151 (MQ=255)
                                  cTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAtttttt                                                                                                                   >  1:677859/1‑151 (MQ=255)
                                      accttgtgCGCTGTGTTTTTTTGCCATTTTTTTGGCTTCTGCGCAAAAAATATGGGTTCCTGACGAAAACCCTTTTTAAATATTGTGCCCCTGCCCCCAGGTCATTGGTGGGAAAAAAAACAGGCCCGCCAGTAAAGCAAAtttttt                                                                                                                   <  2:531051/143‑1 (MQ=255)
                                       aCCTGTGCGCTGTCATTTTTTGCCATTTCGTTGGTTTCTGCGCAAGTAATATAGTTTTCGGCCGGCAACCCTTTTTTAATATTGTAGCCCTGCGCCCAGGTTTTTGGGGCAAAAAAAAACGGCCCCGCCAGTAAAGCAAATTTTTtca                                                                                                                 <  1:510909/148‑1 (MQ=255)
                                           gggcgtTGTCATTTTTTCCCTTTTCTTGGGTTTCTGCGCAAGTAATATAGGTTTCTGGCGACAAACCTTTTTTAATTTTGTAGCCCGCCCCCCAGGTCATTGGTGCGAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatcata                                                                                                            <  2:65870/147‑2 (MQ=255)
                                                 tttCATTTTTTGCTTTTTCTTTGGCTTCTCCGAAAGTAATATAGTTTTCTGACGCCAAACCTTTTTTAATATTGTAGTCCTCCCCCCGGTTTTTTGGTGCAAAAACAACCGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                      <  1:623213/147‑1 (MQ=255)
                                                 tGTCATTTTTTGCTTTTTCTTGGGTTTCTGCGAAAGTAATATAGTTTTCTGACGCCAAACCTTTTTAAATATTGTACCCCTGCCCCCAGTTTTTTGGGGCAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCAtt                                                                                                         <  1:381883/146‑1 (MQ=255)
                                                 tGTCATTTTTTCCCATTTCGTTGTCTTCTGCGCAAGTAATATAGGTTTGTGACGGCAACCCTTTTTTAATATTGTAGCCCTGCGCCCAGTTCTTTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTcc                                                                                                       <  1:157312/148‑1 (MQ=255)
                                                   tctttttttGCTATTTTTTTGTCTTCGGCAAAAAAAATATAGATTTCGAACGGAAACCCTTTTTTAATATTTAAGCCCGGCACACAGTTTATGGGTCGGAAAAAAAAAAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTa                                                                                                    <  1:517293/145‑1 (MQ=255)
                                                       tttttGCCTTTTCGTTGTCTTCTGCGCAAGAAATATAGTTTTCGGAGGGCAAACGGTTTTTAATATTGTGGCCCTGCCCCCAGGTAATGGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                                <  2:1002010/149‑1 (MQ=255)
                                                       tttttGCCATTTCGTTGGCTTCTGCGAAAGTAATATAGGTTTCTGAGGAAAACCCTTTTTTAATTTTGTAGCCCGCCCCCCAGGTCTTTGGGGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                      <  2:196354/143‑1 (MQ=255)
                                                        ttttGCCATTTCGTTGGTTTCTGCGAAAATAATATGGGTTTCTGGCGCCAACCCGTTTTTAATATTGTAGCCCTCCCCCCAGGTCATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAt                                                                                              <  2:844718/150‑1 (MQ=255)
                                                        ttttGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                             >  2:502532/1‑151 (MQ=255)
                                                         tttGCCTTTTCGTGGGTTTCTGCGCAATAAATATAGGTTTCGGACGACAACCCGTTTTAAATTTTGAGGCCCTGCGCCCAGGCTGTTGGTGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                             <  2:365621/150‑1 (MQ=255)
                                                           tGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtttt                                                                                           >  1:978481/1‑150 (MQ=255)
                                                                    gTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAAGCGTTTTTAATATTGTAGTCTTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCTTCATTCCTTATTTCATTTTACCCAGAAtt                                                                                 >  2:58085/1‑151 (MQ=255)
                                                                    gTTGGCTTCTGCACAAAGAATTTAGTTTTTTGCAGGCAACCTGTTTTTAATATTTTAATCCGGCGCCCATTCTTTTGGGGAAAAAACAAAGAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAt                                                                                  <  1:970995/150‑1 (MQ=255)
                                                                     ttGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                                >  2:1156443/1‑151 (MQ=255)
                                                                     ttGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                                >  1:132024/1‑151 (MQ=255)
                                                                     ttGGCTTCTGCGCAAGTAATATAGGTTTCTGACGACAAACTTTTTTAAATTTTGTAGCCCTGCGCCCAGGTCATTGTTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGaa                                                                                   <  1:147713/148‑1 (MQ=255)
                                                                     ttGGCTTCTGCGAAAGAAATATGGTTTTCGGCCGGCAACCTTTTTTAAATATTGTAGTCCGCCCCCCAGGTCATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAt                                                                                  <  2:615388/149‑1 (MQ=255)
                                                                       ggCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                               >  2:583582/1‑150 (MQ=255)
                                                                       ggCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                               >  1:1131994/1‑150 (MQ=255)
                                                                          ttCTGCGAAAGTAATATAGGTTTCTGACGGCAAACCTTTTTTAATATTGTAGCCCTGCGCCCAGGTCTTTGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                               <  2:243648/147‑1 (MQ=255)
                                                                          ttCTGCACAAATAATATAGGTTTCTGACGGAAACCTTTTTTTAATATTGTGCCCCTGCCCCCAGTTTTTTGGTGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGca                                                                              <  1:379778/148‑1 (MQ=255)
                                                                           tCTGCGAAAGTAATATAGTTTCCGGAGGGAAACCCTTTTTAAATTTTGTAGTCTGGCGCCCAGTTCTTTGGGGCAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcacca                                                                           <  1:685124/150‑1 (MQ=255)
                                                                             tGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCATATTTCATTTTACCCAGAATTGCACCAccc                                                                        >  2:685319/1‑151 (MQ=255)
                                                                             ggcgaAAGTAATATGGTTTTCTGACGGAAACCTTTTTTTAATATTGTAGTCCGGCGCCAAGGTCTTTGGTGAAAAAACAACAAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcaccac                                                                          <  1:1081427/148‑1 (MQ=255)
                                                                                gcAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGt                                                                      >  2:960569/1‑150 (MQ=255)
                                                                                  aaGAAATATAGGTTTCTGACGGAAACCCTTTTTAAATATTGTAGCCCGGCCCCCAGGTCATGGGTGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGt                                                                      <  2:530237/148‑1 (MQ=255)
                                                                                   ataaaTATAGTTTTCGGGAGAAAAACGTTTTTTAATGTTGTGGTCCCGCCCCAGATTTTTTGGTGAAAAAAAAACCGGTCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCg                                                                  <  1:963950/148‑1 (MQ=255)
                                                                                      aaTATAGTTTTCTGCCGCAAAACCTTTTTTAATATTGTAGCCCTGCCCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGcc                                                                <  2:77275/150‑1 (MQ=255)
                                                                                       atatAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCa                                                               >  2:706224/1‑150 (MQ=255)
                                                                                       aatttagTTTCCGGCGGGAAACCCTTTTTTAATTTTGTACCCCCGCCCCCAGTTCTTGGGGGCAAAAAAAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGcc                                                                <  2:754835/144‑1 (MQ=255)
                                                                                               ttcgtgaggCAAACCTTTTTTAAATATTGTCGCCCTACCCCCGGTTCTTTGGTCCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACgg                                                        <  1:33530/143‑1 (MQ=255)
                                                                                                  tGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCt                                                   >  2:756675/1‑151 (MQ=255)
                                                                                                       ggaaaCCTTTTTTAAATATTGTATCCCGGCCCCCATGAATTTGGTGGAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCtt                                                  <  2:484775/146‑1 (MQ=255)
                                                                                                       gCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAAt                                              >  1:761791/1‑151 (MQ=255)
                                                                                                       gCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAAt                                              >  1:137278/1‑151 (MQ=255)
                                                                                                         agaCCGGTTTTAATATTGCAGTCCTGCGCCCATCTCATTGGCGCGAAAACAAACACGCTCGCCAGTACAGCAAATTTTTTCATCATCATTCCTTATTTCATCTTACGCAGAAGTGCACCACCCGTACCGCCAATCACGCCACCTTTAAt                                              >  1:1233199/3‑149 (MQ=255)
                                                                                                                 ttttAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                     >  1:497105/1‑151 (MQ=255)
                                                                                                                 tttaAATTTTGTAGCCCGGCCCCCAGGTATTTGGAAAAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCtt                                      <  2:1139282/150‑1 (MQ=255)
                                                                                                                   ttaaaTATTGTAGCCCGGCCCCCAGGACTGTGGGGCAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                     <  2:695893/147‑1 (MQ=255)
                                                                                                                   ttAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                    >  2:934073/1‑150 (MQ=255)
                                                                                                                   ttAATATTGTAGTCCTGCCCCCAGGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                    <  2:545082/150‑1 (MQ=255)
                                                                                                                      aaTATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAgg                                 >  1:121574/1‑151 (MQ=255)
                                                                                                                        ttttGTAGTCCCCCCCCCAGGTTTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                    <  2:753666/144‑1 (MQ=255)
                                                                                                                              aGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGcc                         >  1:319294/1‑151 (MQ=255)
                                                                                                                                  cTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                     >  1:1039245/1‑151 (MQ=255)
                                                                                                                                                 tgggTGCAAAAACAACCGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg         <  2:271799/147‑1 (MQ=255)
                                                                                                                                                     tGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGc    >  1:691263/1‑150 (MQ=255)
                                                                                                                                                      gCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCAc  >  1:365767/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCAC  >  NC_000913/1222447‑1222744

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