Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,288,8520CT27.7% 193.0 / 34.5 95pseudogene (63/2525 nt)yejOpseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (20/48);  new base T (26/0);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 5.16e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGT  >  NC_000913/2288632‑2289089
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
ccATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCCGCCGGCTAGCTTTTTCTAGGTCAACACCATCa                                                                                                                                                                                                                                   >  2:193667/1‑237 (MQ=255)
      ttGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAATAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGacaac                                                                                                                                                                                                                                          <  1:512647/224‑1 (MQ=255)
       tGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCa                                                                                                                                                                                                                          <  1:373695/239‑1 (MQ=255)
       tGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACGGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCa                                                                                                                                                                                                                          <  2:1149966/239‑1 (MQ=255)
         aGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCCGCTGGATATCTTTTTTCAGGACAACAACATCCACTGTCCCAcc                                                                                                                                                                                                                        >  1:53032/1‑239 (MQ=255)
              aGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCCTAAACCAGCTGGATATCTTTTTTCTGGACAACAACATCAACTTTCTCACCAAAc                                                                                                                                                                                                                    >  2:447715/1‑238 (MQ=255)
                tCTGGTTTGCCGCGCCACCCATTTGAATAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCCTGATATCTTTTTTCACGACCACACCACCCACt                                                                                                                                                                                                                                >  2:512647/1‑224 (MQ=255)
                   ggTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACa                                                                                                                                                                                                              <  2:1312330/239‑1 (MQ=255)
                    gTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCa                                                                                                                                                                                                                          <  2:298790/226‑1 (MQ=255)
                    gTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCa                                                                                                                                                                                                                          >  1:298790/1‑226 (MQ=255)
                      ttGCCGCTCCACCCATTTGAACAATTTGATAACCGTAATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGtt                                                                                                                                                                                                           <  2:526794/239‑1 (MQ=255)
                                   cATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTa                                                                                                                                                                                              <  1:447715/239‑1 (MQ=255)
                                    aTTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGcc                                                                                                                                                                                                     <  1:1048643/231‑1 (MQ=255)
                                        gAACAATTTGATAACCGTCATTGAATACTGAGTATTCTGCGGCTCCATTCATTTCAATGTACTGATACCCTCCTTTAATTTCAGTATTATTACTTACACCAAATTCTTTTAGATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAGCCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAAGGCTATGGGt                                                                                                                                                                                         <  2:1155245/239‑1 (MQ=255)
                                              tttGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAATCAGCTTGGTATCTTTTTCCAGTACCACACACTCAACTGTCTCACCATACTCAACAATTTATACCCATGCTCTgt                                                                                                                                                                                           >  2:1048643/1‑230 (MQ=255)
                                              tttGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGtt                                                                                                                                                                                                           <  1:844573/215‑1 (MQ=255)
                                                        gTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCCTAAACCCACTGGAGATCTTTTTTCAGGACAACACCATCAACCTTCTCCCCCTACTCAACAGttttttccaaat                                                                                                                                                                                                 >  2:844573/1‑205 (MQ=255)
                                                            tttAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTg                                                                                                                                                                                                                               <  1:401467/181‑1 (MQ=255)
                                                            tttAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGtc                                                                                                                                                                                                                             <  1:434749/183‑1 (MQ=255)
                                                            tttAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGcc                                                                                                                                                                                                     <  1:936850/207‑1 (MQ=255)
                                                            ttgaaTACTGAGTATTCTGCGGCCCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAACCCAGCTGGATATCTTTTTCCAGGACAACGCCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTa                                                                                                                                                                                              <  1:134973/211‑1 (MQ=255)
                                                                 tACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGacaac                                                                                                                                                                                                                                          <  1:1056492/165‑1 (MQ=255)
                                                                  aCTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATTTTTTTCCAGGACCACAACATTAACTGTCTCACCATACTCAACAGGTGATGCCAATGCTATGGGGGAATATAAGGAAGGTGCGAATAAGCg                                                                                                                                                              >  1:192353/1‑239 (MQ=255)
                                                                      aGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCTGTCCCATCAACCACCTGGATATCTTTTTCCCGAACAACACCATCAACTCTCTCACCACACTcc                                                                                                                                                                                                                 >  2:434749/1‑184 (MQ=255)
                                                                      aGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCCTAAACCAACTGGGTATCTTTTTTCAGGACAACACCATCAACTGTCTCACCCTa                                                                                                                                                                                                                     >  2:401467/1‑181 (MQ=255)
                                                                      aGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAGCCAGCTGGGTATTTTTTTTCCGGACAACACCATCAACCGCCTTACCAAACCTAACAGGTGATGCCAATGATCTGGGGGAATATAAGGAAGGGGCGAATAAGAGGGGa                                                                                                                                                          >  2:1097827/1‑239 (MQ=255)
                                                                      aGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATTAACCAGCTGGGTATCTTTTTTCAGGACAACACCATCAACTTTCTCACCATACATCACACTTGATGCCCATGCTATgg                                                                                                                                                                                           >  2:936850/1‑207 (MQ=255)
                                                                      aGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGGTATCTTTTTTCTGGACAACACCATCAACTGTCTCCCCCTACTCAACAGTTGTTGCCCATGCTATGGGGGAATATAAGGAAGGTGCGCATAAGCGGGGa                                                                                                                                                          >  2:193633/1‑239 (MQ=255)
                                                                      aGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGt                                                                                                                                                                                                                              <  1:386826/172‑1 (MQ=255)
                                                                        tATTCTGCGGGGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGAGGCCAATGCTATggg                                                                                                                                                                                          <  1:528280/206‑1 (MQ=255)
                                                                           tCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCCGCTGGATATCTTTTTTCAGGACAACACCATCAACt                                                                                                                                                                                                                                >  2:1056492/1‑165 (MQ=255)
                                                                                ggcgCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCATTAACCCGCTGGGTATTTTTTTTCAGGACAAAACCATCAACTGTATCACCATACTCAACAGTTGATGTCAATGCTATGGGTGAATATTAAGAAGGGGCGAATACGCGGGGGAATTCAtctc                                                                                                                                                >  2:1325344/1‑239 (MQ=255)
                                                                                ggcgCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTGCCCTAAACCAGCTGGATATCTTTTTTCAGGACAAAAACCTCCACTGTCCCAcccatc                                                                                                                                                                                                                    >  2:386826/1‑168 (MQ=255)
                                                                                       ttcattTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGTATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCaa                                                                                                                                                                                                                                  <  1:1323443/151‑1 (MQ=255)
                                                                                       ttcattTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGt                                                                                                                                                                                         <  1:378711/192‑1 (MQ=255)
                                                                                              cAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGtt                                                                                                                                                                                                           <  1:370014/167‑1 (MQ=255)
                                                                                                 tGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGTATTGATTTTCGTATTATTGGCGGTTCCATTAACCCGCTTGATATCTTTTTTCAGGAAAAAACCATCAACTGTCTCAcc                                                                                                                                                                                                                        >  2:1323443/1‑151 (MQ=255)
                                                                                                 tGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCATAAACCAGCTGGATATTTTTTTCCAGGACAAAACCATCAACTGTCTCACCATACCCAAAAATTTATGCCAATGCTATGGGGGAAAATaagg                                                                                                                                                                               >  2:378711/1‑192 (MQ=255)
                                                                                                   tACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGacaac                                                                                                                                                                                                                                          <  1:614227/131‑1 (MQ=255)
                                                                                                      tGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGAGAATAAACGGGGAAATTCTTCTCGGCTGACTCAGTactttcattt                                                                                                                        >  2:889614/1‑239 (MQ=255)
                                                                                                        aTACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATCTTTTTTCAGGACAACACCCTCAACTGTCTCACCATACTCCACACttgttgtccatt                                                                                                                                                                                                 >  2:370014/1‑161 (MQ=255)
                                                                                                        aTACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAATTGATTCCAAAGCTATGGGTGAATATAAGGAAGGTGGGAATAAGAGGGGAAATTCTTCTTGGCTGACTCAGTcattttatttat                                                                                                                      >  2:178418/1‑237 (MQ=255)
                                                                                                             cTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCATAAACCAGATGGGTATCTTTTTTCTGGACAACACCATCAACTGTTCCCCCCTACTCAAAAGGTGATGCCAATGCTATTGGGGAATATAAGGAAAGTGCGAATAAACGGGGAAATTCTTCTCGGCCGGCTCAGttatgtccattttttatt                                                                                                                 >  2:1071913/1‑228 (MQ=255)
                                                                                                             cTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATCTTTTTTCAGGACAACACCATCAACTGTCTCACCATACTCAAAAATTGATGCCCAAGCTATGGGTGAATATAAGGAAGGTGCGAAAAAGCGGGGTAATTCTTCTCGGCTGGTTCAGGTcttttctttttt                                                                                                                      >  2:751547/1‑233 (MQ=255)
                                                                                                             cTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGacaa                                                                                                                                                                                                                                           <  1:128498/120‑1 (MQ=255)
                                                                                                             cTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGacaa                                                                                                                                                                                                                                           >  2:128498/1‑120 (MQ=255)
                                                                                                             cTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTa                                                                                                                                                                                              >  2:680213/1‑165 (MQ=255)
                                                                                                             cTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTa                                                                                                                                                                                              <  1:680213/165‑1 (MQ=255)
                                                                                                                cGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGGTATCTTTTTCCAGGACAACACCATCAACTGTCTCACCCTACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGCACCAGTCATTTCATTTCTTCATGttt                                                                                                             >  1:299366/1‑239 (MQ=255)
                                                                                                                cGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAg                                                                                                                                                                                                                                                <  1:611772/112‑1 (MQ=255)
                                                                                                                cGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATa                                                                                                                                                                                                                     >  1:1202263/1‑139 (MQ=255)
                                                                                                                cGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATa                                                                                                                                                                                                                     <  2:1202263/139‑1 (MQ=255)
                                                                                                                   tAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                                                                                                          <  1:1149965/239‑1 (MQ=255)
                                                                                                                          aGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCCTAAACCAGCTGGGTATCTTTTTTCAGGACAACACCa                                                                                                                                                                                                                                      >  2:611772/1‑112 (MQ=255)
                                                                                                                            tattatTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCCGCTTGATATTTTTTTTCAGGACAACACCATCAACTGTCTCACCCTACTCAAAAGTTGTTTCCAATGCTATGGGGGAATATAAGGAAGGTGCGAAAAATCGGGGAAATTCTTCTCTGCGGACTCATTAACTTCCTTTTTTCCTGTTTGAGCCGCttttt                                                                                                 >  2:1421385/1‑239 (MQ=255)
                                                                                                                            tattatTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATaagga                                                                                                                                                                              <  1:200052/166‑1 (MQ=255)
                                                                                                                             attattACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                                                <  1:1421385/239‑1 (MQ=255)
                                                                                                                               tattaCTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtct                                                                                              <  1:178418/239‑1 (MQ=255)
                                                                                                                                      tacACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATCTTTTTTCCGGACAACACCATCAACTGTCTCACCAAACTCTACCGTTGATGCCAATGCTATGGGGGAATATTAGGAAGGTGCGCAt                                                                                                                                                                   >  2:200052/1‑166 (MQ=255)
                                                                                                                                         acCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCCACTGGATATCTTTTTTCAGGACAAAACCATCAACTGTTTCACCATACTCAACAATTTATGCCAATGCTATGGGGGAATATAAGGAAGGTGCGAAAAAACGGGGGAATTCTTTTCGGCTGGCTCCGTCATTTCATTTCTTCATGTTTTAGCCGCTTTTTTTTCCCGTaaaa                                                                                     >  2:267721/1‑237 (MQ=37)
                                                                                                                                           cAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAAACCACTGGGTATCTTTTTTCCGGACAAAACCCTCCAATTTCTCACCCTACCCCACAGTTTATGTCCATTCTATTGGGGAATATAAGGGAAGGGCGGATAAACGGGGGAATTCTTCTCTGCGGACTCAGTCAATTTCATTCTTCATGTTTTGTCCGGATTTTTTTTCCGTAAATg                                                                                    >  2:1353934/1‑237 (MQ=255)
                                                                                                                                                  ttttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa                                                                                                                                      <  1:1091920/183‑1 (MQ=255)
                                                                                                                                                    ttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTg                                                                                                                                                                                                                               <  1:820710/93‑1 (MQ=255)
                                                                                                                                                    ttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                                                         <  1:193633/239‑1 (MQ=255)
                                                                                                                                                    ttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                                                         <  1:159595/239‑1 (MQ=255)
                                                                                                                                                    ttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTAACCTAAACCACCTTGGTATCTTTTTTCAGGACAACACCCTCAACAGTCCCACCCTACTCAACAGTTTTTGTCAAATCTATTGGGGAAGATAAGGAAGGGGCGCAACAACGGGGGAAATTTTTTTGGCTGAACTATGCATTTACTTTCTTCTTGTTTGTGGCGGTTTTTTCTTCCGTAAATGCCTTGGATCa                                                                          >  2:1402323/1‑238 (MQ=255)
                                                                                                                                                     tataTGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCACCTGGATATCTTTTTACCGGACAACACCATCCACTGGCTCACCATACTCAACAGTTGTTGCCCAAGCTATGGGGGGATATAAAGAAGGTACGACCAAACGGGGAAATTTTTTTGGGATGGCTCACGCATTTTATTTTTCTCTTGTTGAGGCGGTTTTTTTCTCCGTAAATGCAtttcatcca                                                                         >  2:989417/1‑231 (MQ=255)
                                                                                                                                                      atatGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGcc                                                                       <  1:889614/239‑1 (MQ=255)
                                                                                                                                                         tGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCaa                                                                                                                                                                                                                                  <  1:583055/85‑1 (MQ=255)
                                                                                                                                                          gCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTTTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                                                                   <  1:121433/239‑1 (MQ=255)
                                                                                                                                                            tGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCACCTGGGTATCTTTTTTCAGGACAACACCATCAACCGTCTCACCCTACTCAACAGTTTATTGCAATGCTATGGGGGAATATAAGGAAGGGGCGAATAAGCGGGGGAAATCTTCTCGGCGGGCTCAGTCATTTTATTTCTTCATGTTTGTGCCGCTTTTTTTTCCCCTAAATACCCTGAATTAGCCCAtttt                                                                 >  2:372674/1‑238 (MQ=255)
                                                                                                                                                              ttCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCCGCTGGGTATTTTTTTTCAGGACAACACCATCAACTGTCTCACCCTa                                                                                                                                                                                                                     >  2:820710/1‑93 (MQ=255)
                                                                                                                                                              ttCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGcc                                                                                                                                                                                                     <  2:847492/109‑1 (MQ=255)
                                                                                                                                                              ttCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGcc                                                                                                                                                                                                     >  1:847492/1‑109 (MQ=255)
                                                                                                                                                              ttCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGaaa                                                                                                                                                        <  1:452484/153‑1 (MQ=255)
                                                                                                                                                                   cGCCAGGATTGATTTTCGTATTATTGGCGGTTCCCCAACCCAACAGGGTAACTTTTTTTTGAACAACCCCAGCAACTTTCTTACCAAACTCCAAAGTTGAGGCCAACGCTAATGGCGTATATATGGTAGGGGGGAAAAAATGGGGAAAATCTTTTCGGCGGGCTCCTCGATTTTATTTCTTCATGTTTTGGCCGAGTTTTTCTTCCGCAAAAGCCTTGAATCAGCCTTATTAGTCCGtt                                                          >  2:1053173/1‑239 (MQ=255)
                                                                                                                                                                   cGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGGTATCTTTTTTCAGGACAACACCATCAACTGTCTCAcc                                                                                                                                                                                                                        >  2:583055/1‑85 (MQ=255)
                                                                                                                                                                        ggATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGGTATCTTTTTTCAGGGCAAAACCCTCAACTGTCTCACCCTACCCAACAGTTTGTGCCAAAGCTATGGGTGAATATAAGGAAGGGGCGGATAAACGGGGGAAtttttttagg                                                                                                                                              >  2:452484/1‑148 (MQ=255)
                                                                                                                                                                               ttttCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt                                              <  1:372674/239‑1 (MQ=255)
                                                                                                                                                                                                       aaaCCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGt                      <  1:1097827/239‑1 (MQ=255)
                                                                                                                                                                                                          cccgCTGGATTTCTTTTTCCAGGACACCACCATCAACTGTCTCACCATACTCACCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTgag                   <  1:729623/236‑1 (MQ=255)
                                                                                                                                                                                                             gCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGt                                                                                                                                                                                         >  2:865367/1‑74 (MQ=255)
                                                                                                                                                                                                             gCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGt                                                                                                                                                                                         <  1:865367/74‑1 (MQ=255)
                                                                                                                                                                                                             gCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                                                                   <  2:820810/111‑1 (MQ=255)
                                                                                                                                                                                                             gCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                                                                   >  1:820810/1‑111 (MQ=255)
                                                                                                                                                                                                             gCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc                                <  1:158363/223‑1 (MQ=255)
                                                                                                                                                                                                             gCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc                                >  2:158363/1‑223 (MQ=255)
                                                                                                                                                                                                             gCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTgaga                  >  2:1292671/1‑237 (MQ=255)
                                                                                                                                                                                                              cTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                                                                  <  1:1379106/172‑1 (MQ=255)
                                                                                                                                                                                                              cTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                                                                  >  2:1379106/1‑172 (MQ=255)
                                                                                                                                                                                                              cTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGAtc                >  1:894467/1‑238 (MQ=255)
                                                                                                                                                                                                               tGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAc                                                              >  2:333652/1‑191 (MQ=255)
                                                                                                                                                                                                               tGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAc                                                              <  1:333652/191‑1 (MQ=255)
                                                                                                                                                                                                                         tttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGAc    <  1:1353934/239‑1 (MQ=255)
                                                                                                                                                                                                                           tCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGt  <  1:267721/239‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
CCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGT  >  NC_000913/2288632‑2289089

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