Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,608,2440TG30.2% 50.0 / 11.8 54V443G (GTG→GGG) hyfFhydrogenase 4 component F
Reads supporting (aligned to +/- strand):  ref base T (12/25);  new base G (15/1);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.57e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACG  >  NC_000913/2607999‑2608493
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
ccTGCTCAAGTACGGCCCGCGCGATCTCAACGTCGTCGGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTaa                                                                                                                                                                                                                                                        <  1:271169/250‑1 (MQ=255)
aCTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTa                                                                                                                                                                                                                                                         <  1:254041/250‑1 (MQ=255)
      cAAGTACGGCCCGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGc                                                                                                                                                                                                                                                    <  1:370508/249‑1 (MQ=255)
       aatgaCGGCCCGCGCGATCTCAACGTCTTCTGTGGTATGCGCAAAATCAGGCCATTTACCGCCGTGCTGTTGGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGaa                                                                                                                                                                                                                                                 <  2:518316/247‑1 (MQ=255)
                    gcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGGACGGGTGGCCGCGCGGGGGGTAATGGGGGAACCGCCGCAGGcc                                                                                                                                                                                                                                    >  1:99534/1‑251 (MQ=255)
                             aCGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGAGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCCCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCCCAGGCCGTTAACCgg                                                                                                                                                                                                                          >  1:306063/1‑251 (MQ=255)
                                    cTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAAACGCCGCAAGCCGTTAAACgggggggtc                                                                                                                                                                                                                   >  1:574255/1‑246 (MQ=255)
                                     tgggggaggcacAAAATCAGCCCTTTTCCCGCCGTGCTTTTGGGCGGCGTTGCGCTGGCTCTGGCAGGGTGGCCGCCCTTCAACTTTTTTCTTAGGGAATTTATGACCATTACCGCCGGACTGGCACTTAATCACCGGCTGATTATCGTCCTGCTGTTATGGCTGTTAACGCTGGTGGTGGGGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGAtc                                                                                                                                                                                                                   <  1:321734/239‑1 (MQ=255)
                                          gATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGGTGGAGCTGGCAGGGATGCCGCAATTAAACATTTTTCTTGGCGAATTTATGACCTTTACCGCCGGAATGGCCCGTAATCACCTGATGATTATCGTCCTGCTGTTATTGCTCTTAACGCCGGTCCTGGCGGGCCTGGTACGGTGGGCCGCGCGGGGGTTCATGGGGGAACCCCCGCCGGCCCTTCACCGGGGGGGTCTCGGCt                                                                                                                                                                                                             >  1:655274/1‑251 (MQ=255)
                                            tGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGGACGGGTGGGTGCGCGGGGGGTTATGGCGGAACCGCCGCAGGCCGTTAAACCGGGGGATCTCGGGTgg                                                                                                                                                                                                           >  1:242722/1‑251 (MQ=255)
                                              cTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCAGACAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGAATGGCACGCAATAACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGAAGGCTGGGCGGCCGGGGGTTTATGGCGGAACCGCCGCCGGGCGGTATCCGGGGGGAGCTCgggcgggt                                                                                                                                                                                                         >  1:79298/1‑245 (MQ=255)
                                                    aTCATGCCTTTTCCCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACTTTTTTCTTAGCGAATTTATGCCCATGCCCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGTTGATCTCGGTTGGTTGaccac                                                                                                                                                                                                   <  1:525093/251‑1 (MQ=255)
                                                              ttACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCCGGCCGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAACGCCGGACTGGCACGACCTCACCTGCTGATTATCGTTCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCTGGGGCGGAGGGCGGGGGGGGTGTTAATGGCGGAAACGCCGAAGGGGGTTTATCGGGGGGGTTTTGGCTGGTTGACCCCCCCGcccaa                                                                                                                                                                                          >  1:387423/1‑247 (MQ=255)
                                                                aCCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGATGGCTGCGCGGGGGTTTATGGCGAAACCGCCGGCGGCCGGTAACCGGGGGGGTCTCGGGCGGGTGACCCCCTCGcccaaggg                                                                                                                                                                                       >  1:455509/1‑245 (MQ=255)
                                                                      gTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCCGGGGAGGATGGCTTCGCGGGGGTTTATGGCGGAACCGCCGCCGGGCGTTTACCCGGGGGGTATCGGCTGGtttgccacccccccacaggtggttct                                                                                                                                                                                 >  2:173204/1‑227 (MQ=255)
                                                                       gggtttgtGGCGGGGGTGCGCGGGCGCGGGCGGGGAGGCCCCCCTTCAACATTTTTTTTAGCAAATTTATGCCCATCCCCGCCGGACTGGCCCTTAATCCCCGGCTGTTGATCGCCGTGCTTTTATTGCTTTTAACGCTGGTGCTGGCGGGCCGGGTCGGGAGGGCTGCGGGGGTGTTATGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTct                                                                                                                                                                                 <  2:99534/243‑1 (MQ=255)
                                                                           gttgggcggcGGTGCGCTGGCGCTGGCGGGGATGCCGCCCTTCAATATTTTTCTTACCGAATTTATGACCATTACCGCCGGATTGGCACGTAATCCCCTGCTGATTATCGTCCTGCGTTTATTGCTGTTACCGCTGGTGCTGGCGGGCCTGGAAGGGATGGCTGCGGGGTTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTgg                                                                                                                                                                            <  1:299439/247‑1 (MQ=255)
                                                                                             ggcgctggcAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGGTGACCCCCTCGCCCATGGTGATTCTGCTGGTCATGATGCTGGCGGTgg                                                                                                                                                          >  1:654743/1‑251 (MQ=255)
                                                                                                      aGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCCGGCCGGTAAACGGGGGGATCTCGGCTGGTTGACCACCTCGCCAATGGGGATTCTGCTGGGCATGATGCTGGGGCTGGGAACGCat                                                                                                                                                  >  1:127208/1‑250 (MQ=255)
                                                                                                                 ccTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGAGGGCTGCGCGGGGGTTAATGGCGGAACCGCCGCAGGCCGTTAACCGGGGGGGTCTCGGCTGGGTGACCACCCCGCCCATGGTGATTTTGCTGGGCATGTTTCTGGCGGTGGGGAGGGATATTCCACaaacc                                                                                                                                      >  2:253125/1‑248 (MQ=255)
                                                                                                                 ccTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCCGGCCGTTAACCGGGGGGATCTCCGGTGGTTGACCACCTCGCCCATGGGGATTTTGCTGGGCATGATTCTGGGGATGGGAACGCCGATTCCACAAccc                                                                                                                                      >  1:336611/1‑250 (MQ=255)
                                                                                                                              tCTTAGCAAATTTATGACCATTCCCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATGGCGGTTAACGCTGGTGCTGGCGGGCCTGGTACGGAGGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATc                                                                                                                          <  1:5503/250‑1 (MQ=255)
                                                                                                                                ttggcgatttgagacccttccccgccgGACGGGCACGTAACCCCCTGCTGATAGTCGTCCTGCTGTTATTGCTGTTACCGGTGGTGCGGGCGGCCCGGGTACGGAGGGCTGGGGGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCt                                                                                                                        <  2:446687/229‑1 (MQ=255)
                                                                                                                                 gaGCGAATTTATGCCCATTACCGCCGGACTGGCACGAAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTACCGCTGGTTCTGGCGGGCCGGGTACGGATGGCTGCGCGGTTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATcc                                                                                                                         <  1:90046/247‑1 (MQ=255)
                                                                                                                                  aGCGATTTTATGCCCATTCCCGCGGGACTGGCCCGTAATCACCTGCTGTTTAGCGTCCTGCTGTTTTTGCTGTTAACGCGGGTGCTGTCGGGCCTGTTACGGTTGGCTGCGCGGGTTTTAATGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTgg                                                                                                                      <  1:606552/250‑1 (MQ=255)
                                                                                                                                  aGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGTTGACCACCCCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTggc                                                                                                                     >  1:155839/1‑251 (MQ=255)
                                                                                                                                      attttAGGACCATTCCCCCCGTACTGGCCCGTAATCCCCTGCTGATTATCGTCCTGCTTTTATGGCTGTTAACGCTGGGGCGGGCGGGCCTGGTACGGATGGCGGCGGGGTTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGTGTGATCTCGGCTTGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGTGAACGCATATTCCACAACCTGTCATCAGGATACTggcgggc                                                                                                                 <  1:55665/249‑1 (MQ=255)
                                                                                                                                          tATGACAATTACCGCCGGACTGGCACGTACTCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGGTGGCGGCGCGGGTGGTAATGGCGGAAACGCCGCAGGCCGCTCAACGGGGGGATCTGCGCTGGGTGACCCCCCCCCCCATGGGGCTTCTCCTGGTCCAGACGCTGGCGCTGCGGCCGCCTATTCCACCACCTGTCCACAGGATGCTCGCGGGGGCtt                                                                                                             >  2:208245/1‑251 (MQ=255)
                                                                                                                                           gtGACCATGCCCGCCGGCCTGGCACGTAATCCCCTGTTATTTATCGTCCTGCTGTTTTTGCTGTAAACGCGGTTGCGGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCAAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCtt                                                                                                             <  1:600024/249‑1 (MQ=255)
                                                                                                                                              cccattcccgccggatgggccctaaaTCCCCTGCTGATTATCGTCCTGCTGTTATGCCTTTACACGCTGGTGCTGGGGGGCCTGGTACGGATGGCTGCGCGGGTGTTAAGGGCGAATCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCAc                                                                                                        <  1:281165/227‑1 (MQ=255)
                                                                                                                                               cctttcccgccgGACTGGCACGTAATCCCCTGCTGATTATCGTCCTGTTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCGGGTACGGATGGCTGCGCGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACt                                                                                                       <  1:139061/245‑1 (MQ=255)
                                                                                                                                               ccATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATAGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGCACGGGTGGCAGCGAGGGGGTTAATGGCGCAACCGCCGCAGGCCGTTAAACGGGGGGATCTCGGCTGGGTGACCCCCTCGCCAATGGGGATTCTGCTGGTCatgat                                                                                                                                                                     >  1:307194/1‑190 (MQ=255)
                                                                                                                                                         cgGACTGGCACGTAACCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCGGGCGGGCTTGGTACGTAGGCCGGCGCGGGTGTTAATGCCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                                             <  1:141775/251‑1 (MQ=255)
                                                                                                                                                             cTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGGTTATGGCGAAACCGCCGCAGGCCGGTAAACGGGGGGGTCTCGGGTGGGTGGCCACCCCGCCAATTGGGATTCTGCTGGTCATGAGGGTGGCGCTGGGAACCCATTTTCCACAAACTGTCCTCCGGGACCTGGCGGGGGCTTCCCCCATAATCCtctccgg                                                                                         >  2:166229/1‑248 (MQ=255)
                                                                                                                                                             cTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGGCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGGCCCCTCAgg                                                                                         >  1:531092/1‑251 (MQ=255)
                                                                                                                                                                    gTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATTGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCTCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGCTCCTGGCGGGCGCTTCCCCTTTATTTCTCTCAGGGACGCAc                                                                                  >  1:611572/1‑251 (MQ=255)
                                                                                                                                                                      aaTCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGGGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCCCCATAGTCctcttcggggcgcccgg                                                                                >  1:456607/1‑238 (MQ=255)
                                                                                                                                                                       aTCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAAGCTGGTGCTGGCGGGGCTGGTAAGGATGGCTGCGCGGGGGTTTATGGCGGAACACCCGACGGCCGGTAAACGGGGTGAAATCGGCTGGTTGACCCCCCCGCCAATGGTGATTCTGCTGGTCATGGTGCTGGCGAAGGGAACGCATCTTCCACAACCTGGTATCAGGACCCCGGGGGGCGCTTACCATATAGTCCTCTTAGGGGCGCActct                                                                               >  1:283191/1‑248 (MQ=255)
                                                                                                                                                                        cccccTGCTGATTATCGTCCTGCTGTTATGGCTGTAACCGCTGGTGCTGGCGGCCCTGGTACGGAGGGCGGCGCGGGTGTTAAGGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATc                                                                              <  1:661838/248‑1 (MQ=255)
                                                                                                                                                                          aCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATAGCTGCGCGGGAGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTAGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCCCCATAGTCCTCTCAGGGACGCACCATctg                                                                            >  1:248798/1‑251 (MQ=255)
                                                                                                                                                                            ctgctgATTATCGTCGTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTTCTGGTCATGATGCTGGGGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGGCGCACGATctgcc                                                                          >  1:526607/1‑251 (MQ=255)
                                                                                                                                                                              gttgtttATCGTCCTGCTGTTATGGCGGTTAACGCGGGTCCTGGGGGGCCTGGTAGGGATGGCTGCGGGGGGGTTAGTGGCAAATCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCGATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgcct                                                                         <  1:639918/245‑1 (MQ=255)
                                                                                                                                                                                 gATTATCGTCCTGCTGTTATTGCTGTTACCGCTGGTGCTGGCGGGCCTGGTACGGTTGGCTGCGCGGGTGTAAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                                                      <  1:165172/250‑1 (MQ=255)
                                                                                                                                                                                                   aTTGCTGTTACCGGTGGGGCGGGGGGGCCTGGTACGGAGGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCAt                                                    <  1:379728/250‑1 (MQ=255)
                                                                                                                                                                                                            tagcggtGGTGCGGGCGGGCGTGGTGCGGGGGGCGGTGCGGGTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGTGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCGCCTGGCATGAtttttt                                            <  2:349379/244‑1 (MQ=255)
                                                                                                                                                                                                            taccgttGGTGCGGGCGGTCCTGGTCCGTATGTCGGCGCGGCTGTTATGGGCCAATCCGCCGCAGTCCTTTAACCGGGGTGATCTCGGTTGGTTGACCCCCTCGCGAATGGTGATTCCGCTGGTCATTATGCTGGCGATGTGAGCGCAGATTCTACAGCCTGTCTTCAGGATCCTGGCGGGCGTTTCCACTATAGTGCTCGCAGGCACGCACTATCTGCCTGCACAACGTAGCACCTGGCATGAtttttt                                            <  1:293611/244‑1 (MQ=255)
                                                                                                                                                                                                                   ggTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGTTGACCACCTCGCCAATGGGGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCCCTATAGTCCTCTCAGGGACGCACGAt                                                                               >  1:3779/1‑208 (MQ=255)
                                                                                                                                                                                                                            gggCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCAt                           >  1:501259/1‑251 (MQ=255)
                                                                                                                                                                                                                              gcccTGGTACGGATGGCTGCGCGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATc                          <  1:178071/249‑1 (MQ=255)
                                                                                                                                                                                                                                ccTGGTACGGATGGCTGCGCGGGTGTTACTGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACCACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCCCCGCATCTg                        >  2:131193/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTAGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTc                    <  1:377616/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ggCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCCTGATTTTTTGCCTTCCGGCCCCCCATCTGTTTCTGAGAAac            >  2:128044/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACATCGCCAATGGTGATTATGCTGGCCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCCCTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAacacc         >  1:592796/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 gcgcGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAg        <  1:173204/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       gtgtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACg  >  2:169340/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
ACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACG  >  NC_000913/2607999‑2608493

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