Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,448,200 T→C 5.2% intergenic (‑51/+114) rplN ← / ← rpsQ 50S ribosomal subunit protein L14/30S ribosomal subunit protein S17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,448,2000TC 5.2% 299.6 / 3.2 90intergenic (‑51/+114)rplN/rpsQ50S ribosomal subunit protein L14/30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base T (43/41);  new base C (2/3);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGTGCGCACCACTACCGCCTTCAGCACATCACCTTTTTTGACCTTACCACGCGGAATTGCTTCTTTGATGGTGATCTTGATGATGTCGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGAT  >  NC_000913/3447955‑3448457
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                 
ggTGCGCACCACTACCGCCTTCAGCACATCACCTTTTTTGACCTTACCACGCGGAATTGCTTCTTTGATGGTGATCTTGATGATGTCGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATAACGAGGgcgc                                                                                                                                                                                                                                           >  1:78356/1‑270 (MQ=255)
   gcgcACCACTACCGCCTTCAGCACATCACCTTTTTTGACCTTACCACGCGGAATTGCTTCTTTGATGGTGATCTTGATGATGTCGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGACCTTTTTAGTGCTCCGCTAATGTCAACCACTACTGAGACCCGAAAATCAGGTCGGTAAAAAACCCCAGATCCAGGGGGCgggcg                                                                                                                                                                                                                                      >  2:994455/1‑269 (MQ=255)
        ccacTACCGCCTTCACCACATCCCCTTTTTTGCCCTTACCACGCGGAATTGCTTCTTTGATGGTGATCTTGATGATGTCCCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATaa                                                                                                                                                                                                                                 <  1:61116/272‑1 (MQ=255)
                    tCAGCACATCACCTTTTTTGACCTTACCACGCGGAATTGCTTCTTTGTTGGTGATCTTGATGATGTCGCTTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTATTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTCAGACCCGACACTCAGGTCGTTAAAAATCCCAATATCAAGGGCGCGACATTATAACAcctgtatctt                                                                                                                                                                                                                     >  1:568119/1‑264 (MQ=255)
                    tCAGCACATCACCTTTTTTGACCTTACCACGCGGAATTGCTTCTTTGATGGTGATCTTGATGATGTCGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCAAATACAGGGCCCGGCAATAAAacaacacccccaa                                                                                                                                                                                                                     >  2:917762/1‑262 (MQ=255)
                                 ctttttGACCTTACCACGCGGAATTGCTTCTTTGATGGTGATCTTGATGATGTCGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCAAATCGAgg                                                                                                                                                                                                                                               >  2:654357/2‑233 (MQ=255)
                                  tttttGACCTTACCACGCGGAATTGCTTCTTTGATGGTGATCTTGATGATGTCGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAggg                                                                                                                                                                                                                                              <  1:654358/233‑1 (MQ=255)
                                                      aaTTGCTTCTTTGATGGTGATCTTGATGATGTCGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGgc                                                                                                                                                                                                                                             >  2:729569/1‑214 (MQ=255)
                                                       aTTGCTTCTTTGATGGTGATCTTGATGATGTCGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGgcg                                                                                                                                                                                                                                            <  1:729570/214‑1 (MQ=255)
                                                              tgttGATTGTGATCTTGGTGATGTCGCCTACGCCTGCGTAGCGACGGTGCGCGCCACCCAGAACCTTGATACACATTACGCGCCGTGCGCCGGAGTTGTCAGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                                                           <  1:994456/270‑1 (MQ=255)
                                                                        gATTTTGATGATGTCGCCTCCGCCTGCGTAGCGGCGGTGCGAGCCACCCAGAACCTTGATAAACAGTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATAAGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                                                                 <  1:1029525/272‑1 (MQ=255)
                                                                               atgatgTCGCCGTCGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGCTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgt                                                                                                                                                          <  1:745836/272‑1 (MQ=255)
                                                                                      cGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCCTCAAGGATATGGGTAGAAAAAATAAACGGCCACTTTCTGAGCCGTTTATTCGTATTtagagctgtttctt                                                                                                                                                    >  1:694059/1‑262 (MQ=255)
                                                                                       gCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGCTAAAAATCCCCATACCGAGGGCGCGGCATTATAACACCCCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCCGAGCCGTTTATTCGTAATTAGAAAGTGTACTGTa                                                                                                                                                  >  2:318128/1‑272 (MQ=255)
                                                                                                     aGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGgc                                                                                                                                                                                                                                             >  2:474200/1‑167 (MQ=255)
                                                                                                     aGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATgg                                                                                                                                                                                                              >  1:320475/1‑198 (MQ=255)
                                                                                                     aGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTga                                                                                                                                                                >  1:517006/1‑244 (MQ=255)
                                                                                                     aGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCt                                                                                                                                     >  1:1008390/1‑271 (MQ=255)
                                                                                                     aGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTACTCGTATTGTGAGAGTGTACTTTATTACAGAACCGcct                                                                                                                                    >  2:818711/1‑270 (MQ=255)
                                                                                                     aGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGATCATTTCCAAGCCGTTTATTGGTTTTGAGAGAGTGTACTGTAGTAAAAAACCgact                                                                                                                                    >  2:883405/1‑269 (MQ=255)
                                                                                                      gCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGgcg                                                                                                                                                                                                                                            <  1:474201/167‑1 (MQ=255)
                                                                                                      gCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATggg                                                                                                                                                                                                             <  2:320474/198‑1 (MQ=255)
                                                                                                      gCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTgag                                                                                                                                                               <  1:57068/244‑1 (MQ=255)
                                                                                                          cGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCGGAGCCGTTTATTCGTATTgagag                                                                                                                                                             >  1:273118/1‑242 (MQ=255)
                                                                                                                          gAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGCCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGAAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGt                                                                                                                                                                     >  1:455426/1‑218 (MQ=255)
                                                                                                                          gAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGt                                                                                                                                                                     >  2:904275/1‑218 (MQ=255)
                                                                                                                          gAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAg                                                                                                                 >  1:1026798/1‑270 (MQ=255)
                                                                                                                          gAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAAAGTGTACTGTATTACAGAACCGCTTTCTCAACAACACAAACCAg                                                                                                                 >  1:927187/1‑270 (MQ=255)
                                                                                                                          gAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGCTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCg                                                                                                                                                                                                                            >  1:876912/1‑163 (MQ=255)
                                                                                                                           aaCCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTa                                                                                                                                                                    <  1:904276/218‑1 (MQ=255)
                                                                                                                           aaCCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGCTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGc                                                                                                                                                                                                                           <  2:876911/163‑1 (MQ=255)
                                                                                                                                  aTACACATTACGCGACGTGCACCGGAGTTGTCTGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAgg                                                                                                         <  1:818712/270‑1 (MQ=255)
                                                                                                                                   tACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTaaaaa                                                                                                                                                                                                                                                             >  1:989721/1‑121 (MQ=255)
                                                                                                                                   tACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTa                                                                                                                                                                    <  1:1125740/210‑1 (MQ=255)
                                                                                                                                   tACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTa                                                                                                                                                                    >  2:1125739/1‑210 (MQ=255)
                                                                                                                                   tACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCAAGGAtt                                                                                                      >  1:198888/1‑272 (MQ=255)
                                                                                                                                   tACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACGGTATTAACGAAACGCTTTCTCAAAAAAGAGAAACAGCgatcagaaat                                                                                                      >  2:315964/1‑263 (MQ=255)
                                                                                                                                    acacATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAAt                                                                                                                                                                                                                                                            <  2:989720/121‑1 (MQ=255)
                                                                                                                                    acacATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAgg                                                                                                         <  1:315965/268‑1 (MQ=255)
                                                                                                                                     cacaTTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTTCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACTCCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATGTCTGAGCCGTTTATTCGTATTGAGAGAgt                                                                                                                                                          <  1:1062428/218‑1 (MQ=255)
                                                                                                                                           aCGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGt                                                                                                                                                                     >  2:851950/1‑201 (MQ=255)
                                                                                                                                            cgcgACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTa                                                                                                                                                                    <  1:851951/201‑1 (MQ=255)
                                                                                                                                                     cACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTa                                                                                                                                                                                                                                                                 >  1:426226/1‑99 (MQ=255)
                                                                                                                                                     cACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGc                                                                                                                                                                                                                           >  1:414563/1‑137 (MQ=255)
                                                                                                                                                      aCCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTaa                                                                                                                                                                                                                                                                <  2:426225/99‑1 (MQ=255)
                                                                                                                                                      aCCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCt                                                                                                                                                                                                                          <  2:414562/137‑1 (MQ=255)
                                                                                                                                                        cGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAgg                                                                                                                                                                                                                                               >  2:751563/1‑114 (MQ=255)
                                                                                                                                                         ggAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAggg                                                                                                                                                                                                                                              <  1:751564/114‑1 (MQ=255)
                                                                                                                                                          gAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCAt                                                                               <  1:534876/272‑1 (MQ=255)
                                                                                                                                                                tCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGgc                                                                                                                                                                                                                                             >  2:654090/1‑108 (MQ=255)
                                                                                                                                                                tCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAtt                                                                                                                                                                  >  2:431331/1‑183 (MQ=255)
                                                                                                                                                                tCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTga                                                                                                                                                                >  1:437161/1‑185 (MQ=255)
                                                                                                                                                                 cGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGgcg                                                                                                                                                                                                                                            <  1:654091/108‑1 (MQ=255)
                                                                                                                                                                 cGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTgag                                                                                                                                                               <  2:437160/185‑1 (MQ=255)
                                                                                                                                                                 cGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                                                                 <  1:431332/183‑1 (MQ=255)
                                                                                                                                                                   gCGACGTTCAGCATAGTCTGTTCTTAGATCATTTCAGTCCGCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCTCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGTAGATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTc                                                                                                                                                                       <  1:397044/175‑1 (MQ=255)
                                                                                                                                                                                 aGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTga                                                                                                                                                                >  2:1059947/1‑168 (MQ=255)
                                                                                                                                                                                  gTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTgag                                                                                                                                                               <  1:1059948/168‑1 (MQ=255)
                                                                                                                                                                                      gTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGagagag                                                                                                                                                           >  2:844877/1‑168 (MQ=255)
                                                                                                                                                                                      gTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAAACACGTCAccaaa                                                     >  1:417693/1‑267 (MQ=255)
                                                                                                                                                                                      gTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGCCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCCCATTTCTGAGCCGTTTCTTCGTAtt                                                                                                                                                                  >  2:402264/1‑161 (MQ=255)
                                                                                                                                                                                       ttCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTCCTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACCCCGCTTCAAGGATATTGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                                                                 <  1:402265/161‑1 (MQ=255)
                                                                                                                                                                                       ttCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgt                                                                                                                                                          <  1:844878/168‑1 (MQ=255)
                                                                                                                                                                                           tggttaaTTTTAGTGTCCCGCTAATGTCAACTCCTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCAt                                              <  2:316229/266‑1 (MQ=255)
                                                                                                                                                                                           tGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAACAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTGTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCAt                                              <  1:506538/272‑1 (MQ=255)
                                                                                                                                                                                           tGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGtt                                                                                                                                                                            >  2:357210/1‑146 (MQ=255)
                                                                                                                                                                                           tGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAcgc                                                                                                                        >  1:278184/1‑198 (MQ=255)
                                                                                                                                                                                            ggATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                                                           <  1:357211/146‑1 (MQ=255)
                                                                                                                                                                                            ggATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAcgcg                                                                                                                       <  2:278183/198‑1 (MQ=255)
                                                                                                                                                                                                cATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                                                                 <  1:760689/152‑1 (MQ=255)
                                                                                                                                                                                                cATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAg                                                                                                          >  2:161972/1‑207 (MQ=255)
                                                                                                                                                                                                 aTTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAgg                                                                                                         <  1:161972/207‑1 (MQ=255)
                                                                                                                                                                                                    ttAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGATCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACGGCTTTCCCTACCATCCTAACCAGCGTCCAGGATTTAGTCTTGGTCAGCGGACGGGATTCGCGGAACTCACACACGTCCCCGATAACGCaacggtttggg                                     >  2:506537/1‑262 (MQ=255)
                                                                                                                                                                                                     actgccccggCTAATTTCAATTACTACTGGAACCCGAAAATCAGGGCGTTAAAAATCCCCATATCGAGTGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtct                                    <  2:417692/263‑1 (MQ=255)
                                                                                                                                                                                                          ccccGCAAATTTAAACTATTACGGAGACCTGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGTc                                <  2:78356/269‑1 (MQ=255)
                                                                                                                                                                                                           tCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGagagag                                                                                                                                                           >  1:616050/1‑147 (MQ=255)
                                                                                                                                                                                                            ccGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgt                                                                                                                                                          <  2:616049/147‑1 (MQ=255)
                                                                                                                                                                                                             cataagttcaaaCTACTACGGAGACCCAAAAAGCAGGTCGTTAAAAATCCCCATTTCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGTCATGt                            <  2:1026797/262‑1 (MQ=255)
                                                                                                                                                                                                               cTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGCTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGTCATGTAc                          <  1:318129/272‑1 (MQ=255)
                                                                                                                                                                                                                                  agaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATtata                                                                                                                                                                                                                                  >  1:249902/1‑53 (MQ=255)
                                                                                                                                                                                                                                  agaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGagagag                                                                                                                                                           >  1:1121185/1‑124 (MQ=255)
                                                                                                                                                                                                                                   gaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATaa                                                                                                                                                                                                                                 <  2:249902/53‑1 (MQ=255)
                                                                                                                                                                                                                                   gaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgt                                                                                                                                                          <  2:1121184/124‑1 (MQ=255)
                                                                                                                                                                                                                                       cGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGagagag                                                                                                                                                           >  2:269172/1‑119 (MQ=255)
                                                                                                                                                                                                                                       cGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATGTAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCCTTCGTTTTTCTCGTCATGTACGTTCAGTTTGGTCGTACGCTtgat  >  1:36212/1‑272 (MQ=255)
                                                                                                                                                                                                                                        gAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgt                                                                                                                                                          <  1:269173/119‑1 (MQ=255)
                                                                                                                                                                                                                                          aaaTCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAtt                                                                                                                                                                  >  2:60493/1‑109 (MQ=255)
                                                                                                                                                                                                                                           aaTCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                                                                 <  1:60493/109‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGagagag                                                                                                                                                           >  1:1085343/1‑109 (MQ=255)
                                                                                                                                                                                                                                                  tCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgt                                                                                                                                                          <  2:1085342/109‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                 
GGTGCGCACCACTACCGCCTTCAGCACATCACCTTTTTTGACCTTACCACGCGGAATTGCTTCTTTGATGGTGATCTTGATGATGTCGCCTACGCCTGCGTAGCGACGGTGCGAGCCACCCAGAACCTTGATACACATTACGCGACGTGCACCGGAGTTGTCGGCGACGTTCAGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGAT  >  NC_000913/3447955‑3448457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: