Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,571,0360TG45.7% 93.6 / 15.0 71N239T (AAC→ACC) eutEaldehyde oxidoreductase, ethanolamine utilization protein
Reads supporting (aligned to +/- strand):  ref base T (12/26);  new base G (32/0);  total (45/26)
Fisher's exact test for biased strand distribution p-value = 2.54e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAACCA  >  NC_000913/2570899‑2571186
                                                                                                                                         |                                                                                                                                                             
gTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCCGCGGGGTGCCAGCGcccg                                                                                                                                                   >  2:231042/1‑148 (MQ=255)
  tCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCgg                                                                                                                                               <  2:1083275/151‑1 (MQ=255)
      cGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTgcccgcgcccgccggcat                                                                                                                                           >  1:1633303/1‑145 (MQ=255)
      cGGCCACGCTATCAACAACAATCCGTACCTTTTCGTCGGTACAAATGATGTTGTTATCGAAAGACGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCCCTACCCGCGGGCTGCCTGCGCCTGCGGcctt                                                                                                                                          >  1:310619/1‑148 (MQ=255)
                aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTCCCGCGCCTTCGGCagctcgaccgttt                                                                                                                                >  1:290282/1‑139 (MQ=255)
                aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCACGCCCGCGGCAATCAGaaggttt                                                                                                                               >  2:651180/1‑145 (MQ=255)
                aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACCACCCGCGGGTTTCCAGCGCCTCCgggcaaccgacggtta                                                                                                                                 >  2:1699057/1‑137 (MQ=255)
                 tCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTAt                                                                                                                             <  1:832640/151‑1 (MQ=255)
                    acaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTgg                                                                                                                          <  1:1713001/151‑1 (MQ=255)
                    acaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTgg                                                                                                                          <  1:1448797/151‑1 (MQ=255)
                    acaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCCGCCATCAGAAGCTTTTTgg                                                                                                                          >  1:1146930/1‑151 (MQ=255)
                         aaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGGCAGCGCCTGGGGGAAACAGACGGTTTTTTgtgggt                                                                                                                     >  2:1626021/1‑151 (MQ=255)
                         aaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCCATCAGACGTTTTTTggggggg                                                                                                                     >  2:826802/1‑146 (MQ=255)
                         aaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCAATCAGACGGTTTTTtgtgtgg                                                                                                                     >  1:1706115/1‑150 (MQ=255)
                          aTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCCCTACCGGCGGGTTGCCCCCGCCCCCGGGAAATAAACGTTTTTTtgtggggt                                                                                                                    >  2:1135963/1‑147 (MQ=255)
                          aTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCAATCAGACGGTTTTTTGTGTGtt                                                                                                                    >  1:592824/1‑151 (MQ=255)
                              gtcacTTTTCGTAGGCACAAATGAAGTTGTTATCGAAAGAAGCGCCTTTGACGAAGGACTGAGCGGCACGGGCGAGGGCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCACGTCAGATGTGTTTTGTGGTGtttc                                                                                                                  >  1:3198/5‑149 (MQ=255)
                              gTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATTCACCACTACCGGCGGGGTGCCAACGCCCGCGGCAAACAGACGCTTTTTGGGGGGTTTTcg                                                                                                                >  1:213820/1‑151 (MQ=255)
                              gTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCCGCGCCTGCGGCAATCAGACGGTTTTTTGGGGGTTTTcg                                                                                                                >  2:430709/1‑151 (MQ=255)
                              gTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGGAATCAGACGTTTTTTTGGGTGTTTAcg                                                                                                                >  1:556006/1‑151 (MQ=255)
                                aCCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTAcgcg                                                                                                              <  2:1706995/151‑1 (MQ=255)
                                aCCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGGGGCAATCCAAAGGTTTTTTGTGTGTTTGcggg                                                                                                              >  2:133490/1‑149 (MQ=255)
                                 ccTTTTCGTCGGCACCAATGATGTTGTTATCGATAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGCTTCATCCACCACTACCCGCGGGTTGCCCGCCGCCGCGGCAAACACACCTTTTTTTGTGTGgtttctcc                                                                                                              >  1:1423642/1‑145 (MQ=255)
                                  cTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGGCAGCGCCTGGGGGAATCAGACGTTTATTGGGGGGGTTACccgc                                                                                                             >  2:737227/1‑150 (MQ=255)
                                            gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTCCCGCCCCCGCGGCCAACAGACGGTTTTTGGGGGGGTTTCGCGCCGCCTcttcta                                                                                                  >  1:1281800/1‑151 (MQ=255)
                                            gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGCCAATCGACGGTTATTGGTGTGTTTACGCGCCGCTTctcctt                                                                                                  >  1:449767/1‑150 (MQ=255)
                                            gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGCGATCCGACGTTTATTTGGGTGGTTTCGCGCCGCCTcttcta                                                                                                  >  1:965793/1‑151 (MQ=255)
                                            gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGCCATCAGACGGTTTTTTGGGGGGTTTCGCGCCGCTTctacta                                                                                                  >  1:36596/1‑151 (MQ=255)
                                             cacaAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCtactac                                                                                                 <  2:416628/151‑1 (MQ=255)
                                             cacaAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCtactac                                                                                                 <  1:139632/151‑1 (MQ=255)
                                               caAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCCGCGCCCGCGGCCATCCGACGGTTTTTTGGGGGTTTAAGCGCCGCTTcttctaac                                                                                                >  1:130412/1‑148 (MQ=255)
                                               caAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTCCAGCGCCTTCGGCAATCAGACGGTTTTTTGTGTGTTTTCGCGCCGCCTCTTCTAcc                                                                                                >  2:42408/1‑150 (MQ=255)
                                                 aaTGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCCCCACTACCGGCGGGGTGGCCGCGCCCGCGGCACTCAGACGGTTATTGGTGTGTTTTCGCGCCGCTTCTTCTAccc                                                                                               >  1:612051/1‑148 (MQ=255)
                                                  atgatgTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGGAAACCGGCGTTTTTTGGGGTGTTTTAGCGCCCCCTTTACTACCGCtt                                                                                            >  1:17542/1‑151 (MQ=255)
                                                  atgatgTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAACGCCTGCCGGAATCAAACGTTTATTTGGGTGTTTACGCGCCCCTTTTTCTACCGc                                                                                              >  2:422634/1‑149 (MQ=255)
                                                          gttATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCCATCAGACGTTTTTTTGGGTGTTTTCGGGCCGCTTCTACTACCGCCTTGcccccg                                                                                  >  1:1592272/1‑151 (MQ=255)
                                                          gttATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCAATCAGGCGTTTTTTGGTGTGTTTACGCGCCGCCTCTTCCACCGCTTCGccgccg                                                                                  >  2:948501/1‑151 (MQ=255)
                                                               cGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGGAATCAGACGGTTATTGGGGTGTTTTCGCGCCCCTTCTACTACCGCTTCGCCGCCgggtt                                                                              >  1:598948/1‑147 (MQ=255)
                                                                 aaaGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTacca                                                                           <  1:56340/151‑1 (MQ=255)
                                                                 aaaGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCCGCAATCAGAAGGTTTTTTGGGTAGTTTCCCGCCGCCTTTACTACCGCCTCGCCGCCGGGTatc                                                                            >  2:682373/1‑148 (MQ=255)
                                                                 aaaGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCAATACCGGCGGGTTGGCCCCGCCTTGCGGCAAACGGCGGTTCTTTGGGGGTTTACACGCCGCTTTTACTACCCCCTCGCCGCCCGTTacc                                                                            >  2:1378335/1‑150 (MQ=255)
                                                                  aagaagCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTaccac                                                                          <  2:720899/151‑1 (MQ=255)
                                                                    gaagCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTCCAGCGCCTGCGGCGATCAGACGGTTTTTTGTGTGGTTTCGCCCCGCTTCTACTACCGCCTCCCCGCCGGGTAccacca                                                                        >  1:815716/1‑151 (MQ=255)
                                                                          ccTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCCCTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGc                                                                  <  2:36133/151‑1 (MQ=255)
                                                                          ccTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTTTTGGGGGGTTTTCGCGCCCCCTCTACTACCCCTTCCCCGCCCGGTTCcccccaccggg                                                                   >  2:729748/1‑146 (MQ=255)
                                                                          ccTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGc                                                                  <  1:509252/151‑1 (MQ=255)
                                                                               gACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCAATCAGACGTTTTTTTGTGTGTTTTCGCGCCCCTTCTACTACCCCCTCGCCCCCGGGTACCCCCAGCAGGCCGCTa                                                             >  2:1161471/1‑151 (MQ=255)
                                                                                aCGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCCAACAGACGTTTTTTGGTGTGTTTTCCCGCCGCTTCTACTACCCCTTCGCCGCCCGGTTCCCCCAGCAGGCCGATaa                                                            >  2:424446/1‑150 (MQ=255)
                                                                                aCGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCCCCCGCGGGCAACAGACGGTTTTTTGTGGGTTTACGCGCCGCCTCTACTACCCCTTCCCCGCCGGGTatcacccacacggcgcta                                                             >  1:1313941/1‑138 (MQ=255)
                                                                                     ggACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCCGCGGCAATCAGACGTTTTTTGGGGTGTTTTCGCGCCCCTTCTTCTACCGCCTCGCCGCCGGGTTCCACCAGCAGGCCGATACCCggg                                                       >  2:943147/1‑150 (MQ=255)
                                                                                       aCTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGaaa                                                     <  2:159946/151‑1 (MQ=255)
                                                                                       aCTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGaaa                                                     <  1:708827/151‑1 (MQ=255)
                                                                                               gCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCAATCAGAAGGTTTTTGGGGTGGTTTCGCGCCGCTTCTACTAACCCCTCCCCGCCGGTTAACACCAGCAGGCCGGTACCCCGGAACCTGGACa                                            >  1:970228/1‑151 (MQ=255)
                                                                                                cACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGATTACGCGCCGCTTCTACTACCGCTGCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACa                                            <  1:1175635/150‑1 (MQ=255)
                                                                                                 aCGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAg                                          <  2:636349/151‑1 (MQ=255)
                                                                                                 aCGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAg                                          <  1:836218/151‑1 (MQ=255)
                                                                                                    ggCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTCCAGCGCCCGCGGCCATTAGACGGTTATTTGGGGGGTTTCGCGCCCCTTCTCCTACCGCCTCGCCGCCGGGTTACCCCCGCAGGCCGGTACCCGGGAACTTGGAACAACCg                                        >  1:1034406/1‑151 (MQ=255)
                                                                                                          gTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTgcg                                   <  2:1531528/149‑1 (MQ=255)
                                                                                                                 gTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAAAAAGCGGTGCGCGGGTTCGa                          >  2:22381/1‑151 (MQ=255)
                                                                                                                   ttCATCCGCCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGata                        <  1:1682694/151‑1 (MQ=255)
                                                                                                                   ttCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGata                        <  2:36617/151‑1 (MQ=255)
                                                                                                                   ttCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGAAAACTTGAACAAGCGTTGCGCGGTTTCGata                        <  2:599209/151‑1 (MQ=255)
                                                                                                                   ttCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGGGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGata                        <  1:1426628/151‑1 (MQ=255)
                                                                                                                      aTCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACgcgc                                                                                                             <  1:124234/66‑1 (MQ=255)
                                                                                                                      aTCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGGGTTTACgcgc                                                                                                             >  2:124269/1‑66 (MQ=255)
                                                                                                                           ccacTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGAtt                <  2:128515/151‑1 (MQ=255)
                                                                                                                             acTACCGGCGGGGTGCCAGCGCCCGCGGCAATCAGACGTTTATTGGTGGGGTTACGCGCCGCCTCCACCAACGCTTCGCCCCCGGGTAACACCAGCAGGCCGATACCCGGGAAACTGAACTAACGTTGCGCGGTTTCGCTATTCGGAtttt              >  1:1201942/1‑150 (MQ=255)
                                                                                                                               tACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTAGTGGTGTGTTTACGCGCCGCTGCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGcc            <  1:1279746/151‑1 (MQ=255)
                                                                                                                               tACCGGCGGGTTGCCAGCGCCGGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGcc            <  2:1642767/151‑1 (MQ=255)
                                                                                                                                     cggGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTa      >  1:881919/1‑151 (MQ=255)
                                                                                                                                     cggGTTGCCAGCGCCTGCGGCAATCAGACGGTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGt       <  1:327420/150‑1 (MQ=255)
                                                                                                                                         ttGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAAcc   <  1:1625212/150‑1 (MQ=255)
                                                                                                                                         ttGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAAcc   <  2:593082/150‑1 (MQ=255)
                                                                                                                                         ttGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAAcc   <  1:829237/150‑1 (MQ=255)
                                                                                                                                         ttGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAACCa  >  1:28543/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                                             
GTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAACCA  >  NC_000913/2570899‑2571186

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