Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,185,0560TC35.9% 61.0 / 64.8 64Y180C (TAC→TGC) potAspermidine preferential ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (17/24);  new base C (10/13);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGACAGTGGCGTTAAACATATTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCT  >  NC_000913/1184819‑1185306
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          
aTGACAGTGGCGTTAAACATATTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTga                                                                                                                                                                                                                                               <  2:209362/251‑1 (MQ=255)
  gACAGTGGCGTTAAACATATTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTACTGGTCGTGAGTCACAAAGACGAATGAAATGCCAAGCTTACGCTGTAACGCTTTAAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTg                                                                                                                                                                                                                                                >  1:404106/1‑248 (MQ=255)
       gggCGTTAAACATATGGATGTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGAc                                                                                                                                                                                                                                        <  2:40/250‑1 (MQ=255)
         gCGTTAAACATATTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCCTCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTc                                                                                                                                                                                                                                      >  1:94161/1‑251 (MQ=255)
         gCGTTAAACATATTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTc                                                                                                                                                                                                                                      <  2:99959/251‑1 (MQ=255)
          cGTTAAACATATTGATGTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCAGCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCa                                                                                                                                                                                                                                     <  2:359278/251‑1 (MQ=255)
           gTTAAACATATTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCAt                                                                                                                                                                                                                                    <  1:71924/251‑1 (MQ=255)
              aaaCATATTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCTGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCa                                                                                                                                                                                                                                 <  2:363223/251‑1 (MQ=255)
               aaCATATTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATcc                                                                                                                                                                                                                                  >  1:237324/1‑249 (MQ=255)
                  atatTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGcaa                                                                                                                                                                                                                             <  1:116313/251‑1 (MQ=255)
                    atTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGcaac                                                                                                                                                                                                                            >  2:116193/1‑250 (MQ=255)
                                     aGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTGGCTTACGCAGTTTGTAAACCAGCGTTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTa                                                                                                                                                                                                          <  1:278086/251‑1 (MQ=255)
                                          gCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAaccac                                                                                                                                                                                                     >  2:442382/1‑251 (MQ=255)
                                           caacaaACAGGTTTTTCGGCTCTTCGTAGATTTCGCGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTTGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAaccacc                                                                                                                                                                                                    <  2:393580/251‑1 (MQ=255)
                                           caacaaACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAaccacc                                                                                                                                                                                                    >  1:273120/1‑251 (MQ=255)
                                           caacaaACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAaccacc                                                                                                                                                                                                    <  1:94208/251‑1 (MQ=255)
                                                       ttttCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAAt                                                                                                                                                                                         <  1:168046/250‑1 (MQ=255)
                                                            ggCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGAc                                                                                                                                                                                   >  2:359503/1‑251 (MQ=255)
                                                               tctTCGTCGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGAcgcg                                                                                                                                                                                >  2:408303/1‑251 (MQ=255)
                                                                  tCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGAGACAATCGCGAATCACCACAATCATGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGCCGAATGTAATGCCAAGCATACGCTGTACCGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCCATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGATTGTTAACCACCGCGCGAGCAATGGCGACGCGTTg                                                                                                                                                                             >  1:287491/1‑251 (MQ=255)
                                                                     tAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCCATGGCGACGCGTTGCTg                                                                                                                                                                          >  2:289018/1‑251 (MQ=255)
                                                                      ggatggcACGCGGCGGGCCGGCTTGCTCAAGGCGACCATCGCGCCTCACCACAGTCCTGTCTGGCATGGTGAGTTCTTCTTCCTGGGCGGGAGTCACAAAGACGAAGGTAATGACAAGCTTACGCTGGAACGCTTGCAGCTCGTTCTGCATTTGCGTACGCAGTGTGTAATCCATCGCTTAGAGCGACTCATCCAGCAACAACAGACGAGGCGTGTAAACGACCGCGCTAGCAATGGCGACGCGTTGCTGt                                                                                                                                                                         <  2:63686/245‑1 (MQ=255)
                                                                             aCGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGaccac                                                                                                                                                                  >  2:283580/1‑251 (MQ=255)
                                                                              cgcgGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGGTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGaccacc                                                                                                                                                                 <  1:293027/251‑1 (MQ=255)
                                                                                 ggcgTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCaga                                                                                                                                                              <  1:138917/251‑1 (MQ=255)
                                                                                   cgTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGATAGCGACTCATCCGGCAACAACAGACGAGGCTTGTGAACCGCCTGGCGAGCAATGGCGACGCGTTGCTGTTGACCACCagaga                                                                                                                                                            <  1:437775/251‑1 (MQ=255)
                                                                                      gCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCt                                                                                                                                                         >  1:346672/1‑251 (MQ=255)
                                                                                        cGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGAAGCGTTGCTGTTGACCATCAGCGAGCTg                                                                                                                                                        >  2:378208/1‑250 (MQ=255)
                                                                                          tCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCACTGGCGACGCGCTGCTGGTGACCACCAGAGAGCTGagg                                                                                                                                                     >  2:255240/1‑249 (MQ=255)
                                                                                                 cAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGttt                                                                                                                                               <  1:115196/250‑1 (MQ=255)
                                                                                                 cAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTg                                                                                                                                              <  2:109079/251‑1 (MQ=255)
                                                                                                       gACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGa                                                                                                                                        <  1:278036/251‑1 (MQ=255)
                                                                                                        cccATCGCGCATCACCACAAGCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTTCGTTGAg                                                                                                                                       <  2:23175/250‑1 (MQ=255)
                                                                                                        aCCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAg                                                                                                                                       <  1:283580/251‑1 (MQ=255)
                                                                                                         ccttCGCGCATCGCCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGGCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGCCGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGc                                                                                                                                      <  2:91704/248‑1 (MQ=255)
                                                                                                             cgcgCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTTATGCGGTTTGCGTTGAGCGAAg                                                                                                                                  <  1:79847/251‑1 (MQ=255)
                                                                                                                 cATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGttt                                                                                                                              >  1:54262/1‑251 (MQ=255)
                                                                                                                  aTCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGCGCGAAGGTTTc                                                                                                                             >  2:374678/1‑251 (MQ=255)
                                                                                                                      ccacAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAAc                                                                                                                         <  1:442382/251‑1 (MQ=255)
                                                                                                                       cacaATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACt                                                                                                                        >  1:111858/1‑251 (MQ=255)
                                                                                                                                gTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACAAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATcc                                                                                                               >  2:210739/1‑251 (MQ=255)
                                                                                                                                  cTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGc                                                                                                             >  1:339583/1‑251 (MQ=255)
                                                                                                                                      cATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCATCTGCACCGTCCGCAggg                                                                                                         <  1:359685/251‑1 (MQ=255)
                                                                                                                                        tGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCt                                                                                                       >  1:118886/1‑251 (MQ=255)
                                                                                                                                          gTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGGTTCCAACTGCACCATCCGCAGGGCTTc                                                                                                     <  1:408303/251‑1 (MQ=255)
                                                                                                                                            gAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCACCCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCCCCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTcc                                                                                                    >  1:44882/1‑250 (MQ=255)
                                                                                                                                              gTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAATCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGTCTTGTTAACCGCCGCGCGAGCAATTGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGTTTTGCGTTGTGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATc                                                                                                 <  2:473332/251‑1 (MQ=255)
                                                                                                                                                 cttcttCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCAcg                                                                                              >  1:225238/1‑251 (MQ=255)
                                                                                                                                                      tcCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCgg                                                                                           <  1:201432/249‑1 (MQ=255)
                                                                                                                                                                         aGACGCATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTGCGCAGTTTGTAAGCCAGCGCTGAGAGCGACGCAGCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTGGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCGTCCACCACGCGGGGCGGAATTTCAGCAGCgggg                                                                      <  2:44882/251‑1 (MQ=255)
                                                                                                                                                                            cGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCCGCGGGGGtt                                                                   >  2:192201/1‑251 (MQ=255)
                                                                                                                                                                               aTGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGttttt                                                                <  1:289018/251‑1 (MQ=255)
                                                                                                                                                                                 gTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGtttttgc                                                              <  1:324499/251‑3 (MQ=255)
                                                                                                                                                                                              tACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGCCCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAAcccc                                                 >  2:344576/1‑249 (MQ=255)
                                                                                                                                                                                                      aaCGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGCGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGCAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCcac                                         >  2:24065/1‑250 (MQ=255)
                                                                                                                                                                                                              cAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGCCGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGa                                 <  2:34554/250‑1 (MQ=255)
                                                                                                                                                                                                                   cGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATGTTCGAACACgg                           <  1:374678/251‑1 (MQ=255)
                                                                                                                                                                                                                      tCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGGAATTTCAGCAGCGGGG‑TTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTc                         <  1:48894/250‑1 (MQ=255)
                                                                                                                                                                                                                         gCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCAt                       <  2:94161/249‑1 (MQ=255)
                                                                                                                                                                                                                              tGCTTCCGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTgggg                 <  1:116193/250‑1 (MQ=255)
                                                                                                                                                                                                                                  tACGCAGTTTGCAATCAAGCGATGAGAGAGACTCATCCAGCAACAACAGACGAGGCTTGTTAAACACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGAGGTTTGCGTTGAGCGAAGGTTTCCAACTTCACCATCCGCAGGGCTTCCATCACGCAGGTCGTCATTTCAGCAGCGGGGGGTTTT‑GCATGCGCAACACAAAGGACACAGTTTCGATCACGGTCAtggtgtggaca             >  2:393019/1‑241 (MQ=255)
                                                                                                                                                                                                                                          ttGTAATCCAGCGCTTAGAGCGACTCATCCAGCAACACCAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGGTGCTGTTGACCACCAGAGAGCTGATGCGGTTTTCGTTGAGCGAAGGTTTCCAACTGCGCCATCCGCAGGGCTTCCAGCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGGCACATTTTCGAACTCGGTCATGTGGGGGAAAAGTGCGtt     <  2:241820/250‑2 (MQ=255)
                                                                                                                                                                                                                                           tGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGAGTTGCTGTTGACCACCCGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCAGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCCGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGc   >  2:319162/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCt  <  1:69536/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          
ATGACAGTGGCGTTAAACATATTGATTTCGCCAATGAAGCCGGCAACAAACAGGTTTTTCGGCTCTTCGTAGATTTCACGCGGCGTGCCGTCTTGCTCAATGCGACCATCGCGCATCACCACAATCCTGTCTGACATGGTGAGTGCTTCTTCCTGGTCGTGAGTCACAAAGACGAATGTAATGCCAAGCTTACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCT  >  NC_000913/1184819‑1185306

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