Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,390,0540GA21.4% 143.5 / 10.3 70E179K (GAA→AAA) aaeRDNA‑binding transcriptional activator AaeR
Reads supporting (aligned to +/- strand):  ref base G (39/16);  new base A (0/15);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 4.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATTGGCTGTTCTTCAACTATGGCACAAAATGTTCTCGCCGGGCTGACAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGTGCCGCTGATGTGGGTGATCAACGAGATCAATCGTGGGGAGCTGGAGATCCTGCTGCCGCGTTACCAGT  >  NC_000913/3389804‑3390303
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                           
cATTGGCTGTTCTTCAACTATGGCACAAAATGTTCTCGCCGGGCTGACAGCCAAAATGCTGAAAGAATAACCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGAACCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCCGGGCGCGATGCCAATGGTGGCGTGCGCCGCGAAAAGCTATCTCCCACAATACGGCATACCgg                                                                                                                                                                                                                                                             >  2:381808/1‑251 (MQ=255)
     gCTGTTCTTCAACTAGGGCACAAAATGTTCTCGCCGGGCTGACAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGaaaaa                                                                                                                                                                                                                                                        <  1:427771/251‑1 (MQ=255)
        gTTCTTCAACTATGGCACAAAATGTTCTCGCCGGGCTGACAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAAc                                                                                                                                                                                                                                                      <  1:327492/249‑1 (MQ=255)
             tcAACTTTGGCACAAATGTTTCTTGCCGGGCTGACAGCCAAAATGTTGAAAAAATACCCAGGTTTGAGCTTCTATCGGGTTACCGGATATCCAGCCCCCTCCCTGATTGCCGACGGTCTGGATTTGGTGATCCGCGTCGGCTCGTTGCAGGATTCCAGCCGGTTTTCCCGCCGTCTGGGCGCGATGCCAATGTTGTTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGAAAAAACccgcc                                                                                                                                                                                                                                                <  2:476853/249‑1 (MQ=255)
                aCTATGGCACAAAATGTTCTCGCCGGGCTGACAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGAttt                                                                                                                                                                                                                                           >  2:176229/1‑251 (MQ=255)
                      gCACAAAATGTTCTCGCCGGGCTGACAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGagt                                                                                                                                                                                                                                       >  1:130816/1‑249 (MQ=255)
                      gCACAAAATGTTCTCGCCGGGCTGACAGCAAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGCCCGGATTGCCGAGGGTTTGGATGTGGTGACCCGCGTCGGCGCTTTGCAGGATTCCAGCCGGTTTTCCCGCCGTCGGGGGGCGATGCAAATGGGGGTGTGCCCCGGGAAAAGCTACCTCACACAATCCGGCATACCGAAAAAACCCGCCGATTTGagtag                                                                                                                                                                                                                                     <  1:339446/251‑1 (MQ=255)
                         caAAATGTTCTCGCCGGGCTGACAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGtc                                                                                                                                                                                                                                   >  1:88202/1‑250 (MQ=255)
                                    ccGCCGGGCTGACAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTAAATCTGGTGACCGGAATTCCGGCCCCCGCCCTGATTGCCGACGGTCTGGTTGTGGTGACCCGCGTCGGGGCGTTGCAGGATTCCAGCTTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGGGGTGTGCGCCGCGAAAAGCTATCTCACACAAGCCGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCtt                                                                                                                                                                                                                        <  1:111410/250‑1 (MQ=255)
                                    ccGCCGGGCTGACAGCAAAAACGCCGAAAGAATCCCCAGTTTTGCGCGTCAATTTGTTTACCGGATTTCCAGCCCCGGCCCGTATTGCCGGGGGTCTGGATGTGGTGATCCGCGTCGGGGCGTTGCGGTATTCCAGCCTTTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACACTCCGGCATACCGAAAAACCCCGCCGATTTGAGTAGTCATTCATGGCtt                                                                                                                                                                                                                        <  1:350/250‑1 (MQ=255)
                                     cGCCGGGCTGACACCCAAAATGCTGAAAAAATCCCCAGGTTTGAGCGTCAATTTGGTTCCCGAAATCCCAGCCCCCGACCTGTTTGCCGAGGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATCCCAGCTTTTTTCCCCGCCGTCTGGGCGCGAGCCCAATGGTGGTGTGCGCCGCGAAAAGTTATTTCACACAATCCGGCATACCGGAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTg                                                                                                                                                                                                                       <  1:239248/251‑1 (MQ=255)
                                          ggCTGACAGCAAAAATGCGAAAAGAAACCCCAGGTTGGGGCGTAATGTGGGTTGCCGGAAGTCCCGCCCCCGCCGTGTTGGCGGACGGTCTGGAGGTGGTGGCCCGCGTCGGGGGTTTGCAGAATTCCAGCCTGTTTTCCCGCCGTCTGGGCGGGATGCGAATGGGGTTGTGCGCCGCGAAAAGATATCTCACACAATACGTCTTACCGAAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATAc                                                                                                                                                                                                                  <  2:422310/251‑1 (MQ=255)
                                             tGACAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATATGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGCATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGCTTTGAGTATACATTCATGGCTTGAATACAGc                                                                                                                                                                                                               >  2:35930/1‑251 (MQ=255)
                                                cAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTg                                                                                                                                                                                                            >  1:413699/1‑251 (MQ=255)
                                                   ccAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGCTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGTTGGTGTGCGCCGCGAAAAGCTCTCTCACACAATACGGCATCCCGGAAAACCCCGACGGTTTGAGTAGTCATTCATAGCTTGAATAAAGCGAGCgg                                                                                                                                                                                                         >  2:321992/1‑251 (MQ=255)
                                                       aaTGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCg                                                                                                                                                                                                     >  2:543018/1‑251 (MQ=255)
                                                          gCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGCTTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGATGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGCTTTGAGTAGTCATTCATGGCTTGACAACCGCGTGCGGCCCgccc                                                                                                                                                                                                  >  2:76772/1‑248 (MQ=255)
                                                           cTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACa                                                                                                                                                                                                  >  2:49412/1‑250 (MQ=255)
                                                            tGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAAt                                                                                                                                                                                                >  2:228710/1‑251 (MQ=255)
                                                               aaGAATCCCCAGGTTTGAGCGTCAATCTGGTTACCGGATTTCCAGCCCCCGCCCTGATTGCCGACGGTCTGGATTGGGTGATCCGCGTCGGGGCTTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGAGGCCAATGGTGGTGTGCGCCGCAAAAAGCTATCTCACACAATACGGCATACCGAAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGaa                                                                                                                                                                                             <  1:47649/251‑1 (MQ=255)
                                                                aGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATgcat                                                                                                                                                                                            >  1:75962/1‑248 (MQ=255)
                                                                   aTACCAAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGTCGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTAATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTg                                                                                                                                                                                         >  2:108580/1‑251 (MQ=255)
                                                                         aGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTg                                                                                                                                                                                    >  1:206201/1‑250 (MQ=255)
                                                                              tGAGCGTCAATTTGGTTGCCGGAATCCCCGCCCCCGCCCTGATTGCAGCGGGTCTGTTGTGGTTGTCCCGCGTGGGCGCGTTGCAGGTTCCCAGCCTGTTGTCCCGCCTCCTGGGCGCGAGGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATCCGGCATACCGGAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGc                                                                                                                                                                               <  2:348459/250‑1 (MQ=255)
                                                                               gAGCGTAAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTGGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATCCGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCa                                                                                                                                                                              <  1:347888/250‑1 (MQ=255)
                                                                                aGCGTCAATCTGGTTACCGGAATTCCAGCCTCCGACCTGATTGCCGACGTTCTGGATGGGGTGATCCGCGTCGGCGCGTTGCAGGATGCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGAGTGCGCCGCGACAAGCTCTCTCACACAATACGGCCTACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTTAATACAGCGTGCGGCCCGACACAGAATTTGAACTGAACGCAcc                                                                                                                                                                            >  2:501336/1‑251 (MQ=255)
                                                                                 gCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCg                                                                                                                                                                           >  1:236899/1‑251 (MQ=255)
                                                                                       aTCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCACGATGCAGGATTCAAGCCTGTTTTCCCGCCGTCTGGCCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGCTTTGAGTAGTCATTCATGGCATGAATACAACGGGCGGCCCGACAATGAATTTGCACTGATCGCACCGGAAggg                                                                                                                                                                     >  1:449749/1‑251 (MQ=255)
                                                                                           ggTTCCCGGATTTCCGGCCCCCGCCCGGTTTGCCGTCGGTTGGGGTGTGGTGTTCCGCGTCGGCGCGTTGAAGGATTCCAGCCTTTTTCCCCGCCGTCTGGGCGCGATGCAAAGGGTGGTGTGCGCCGCGAAAAGCTATCTCACAAAATCCGGCATACCGAAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGAtct                                                                                                                                                                 <  2:358181/251‑1 (MQ=255)
                                                                                                 cgggattcccAGCCCCCAACCTGTTGGGCGACGGTGTGGATGTGGAGACCGGCGTCGGCGTTGGGAGGGATACCGGCCTGTTTCCCCGCCGTCGGGGCGCGATGCAAGTGGGGGGGGGCGCCGCAAAAAGTTATCTCACAAAATACGGCATCCGAAAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACt                                                                                                                                                            <  2:398120/243‑1 (MQ=255)
                                                                                                     aaTTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGCTCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGcc                                                                                                                                                        >  1:375674/1‑251 (MQ=255)
                                                                                                              ccccGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCcac                                                                                                                                               >  1:533876/1‑251 (MQ=255)
                                                                                                                     cttaTTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGCAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGGACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCAAAAGGa                                                                                                                                          >  2:23367/4‑249 (MQ=255)
                                                                                                                       gATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGa                                                                                                                                          >  1:476198/1‑247 (MQ=255)
                                                                                                                                 tgacTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGACAACCCCGCCGATTTGAGTAGTCATTCATGGCGTGAATACAGCATGCGGTCCGTCAATGTATTTGAACTGATCGCACCGGAAGGGATCTTGACTCTCCTGATCTCACAGGGAGGATTTGTGACTa                                                                                                                             >  2:42949/4‑250 (MQ=255)
                                                                                                                                 ggtgtgtttttggtgATCCGCTTCGGCGCGTCGCAGGTTTCCACCCTGTTTCCCCGCCTTTGGGGCGCGATGCAAATGGTGTTGTGCGCCGCAAAAACCTATCTCACACATTACGGCGTCCCGAAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTaa                                                                                                                            <  1:109915/241‑1 (MQ=255)
                                                                                                                                  gTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTCCCCGCCGTCTGGGCGCGATGCCAATGGGGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAAt                                                                                                                           <  1:116886/251‑1 (MQ=255)
                                                                                                                                    cTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATg                                                                                                                          >  1:194124/1‑250 (MQ=255)
                                                                                                                                           gtggtactcgggtgggccgggTTGAAGGATGCCAGCCTGTTTTCCGCCCTTCGGGGCGCGATCCAAGTGGTGGTGGGCGCCTGAAAAAGCTATCTCACACAATCCGGCATCCAGAAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGAt                                                                                                                   <  2:258457/231‑1 (MQ=255)
                                                                                                                                           gggtgATCCGCGTCGGCGCGTTGAAGGATCCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGGGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATg                                                                                                                  <  1:144652/250‑1 (MQ=255)
                                                                                                                                                     gcgTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTg                                                                                                         >  1:204828/1‑251 (MQ=255)
                                                                                                                                                            cgcgTTGCGGGATTCCAGCCTGTTTTCCCGCCGTCTGGTCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATCCGGCATACCGGAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCgctggc                                                                                                  <  1:86641/251‑1 (MQ=255)
                                                                                                                                                             gcgTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGGGACTAATGATCCGATGACGCTGGTGCgctggct                                                                                                 >  2:438243/1‑251 (MQ=255)
                                                                                                                                                                 tGCAGGTTCCAGCCCTGTTTCCCCGCCGTCTGGGCGCGATGCAAATGGTGGTGTGCGCCGCGAAAAGCTATCTCAAACAATACGGAATACCGAAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGAc                                                                                              <  1:319645/250‑1 (MQ=255)
                                                                                                                                                                 ggCAGGATGCCAGGCTGTTTCCCCCCTGGCGGGGCGCGAGGCGAAGGGTGGTGTGCGCCAGAAAAAGATATCTCAAACAGTACGGCATCCCGGAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGAc                                                                                              <  2:164056/249‑1 (MQ=255)
                                                                                                                                                                 cgaaGGATTCCGGCCTGTTTTCCCGCCGGCTGGGCGCGGACCCGAGGGGGGTGTCCCCGGGGAAAGCCTATCTCCCACAATGGGGCATACAAAAAAAGCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGGCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGAc                                                                                              <  1:180535/247‑1 (MQ=255)
                                                                                                                                                                   cAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGAcggc                                                                                           >  1:527838/1‑251 (MQ=255)
                                                                                                                                                                              gCCTGTTTACCCGCGTTCTGGGCGCGATGCCAATGTGGTTGTGCCCCGCGAAAAGCTATCCACCACAATACGCCTTACGGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCggg                                                                                 <  1:381808/251‑1 (MQ=255)
                                                                                                                                                                                cTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAGACCCGCCGATTTGAGTAGTCAGTCATGGCTTGAATACGGCGTGCGGCCCGACAATGAATTTGATCTGATCGCACCGGAAGGGAGCTCGACGCGCCTGATCCCACAAAGAAGGTTTGTGACTAATGATCCGAGGACGCTAGAGCGCAGGCTGActgctgggtccaggg                                                                                >  1:271377/1‑236 (MQ=255)
                                                                                                                                                                                cTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGAt                                                                               >  1:115967/1‑251 (MQ=255)
                                                                                                                                                                                  tGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGGGCCGGGATc                                                                              >  1:276720/1‑251 (MQ=255)
                                                                                                                                                                                  ggttgtccctcctgcggggagcgATGCAAATGGTGGTGGGCGCCGCAAAAATCTATCTCACACAATACGGCATACAGAAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATc                                                                              <  2:386533/231‑1 (MQ=255)
                                                                                                                                                                                       tCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGGCggggttgcctc                                                                         >  2:253614/1‑243 (MQ=255)
                                                                                                                                                                                       ccccGCCGTCGGGGCGGGATGCAAAGGGTGGTGGGCCCCGCGAAAAGCTATCTCACACAAACCGGCATCCCGAAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCt                                                                          <  2:194124/249‑1 (MQ=255)
                                                                                                                                                                                              gTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGGGGCGGGAttgactaccggccg                                                                  >  2:105466/1‑238 (MQ=255)
                                                                                                                                                                                                        gATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCTCCTAAGCGCCGCTGATGTggg                                                        >  1:53209/1‑251 (MQ=255)
                                                                                                                                                                                                         aTGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGAc                                                                                              >  2:438082/1‑212 (MQ=255)
                                                                                                                                                                                                         aTGCCAAGGGGGGTGTGCCCCGCGAAAAGTTATCTCCCACAATCCGGCATACCGAAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGAc                                                                                              <  1:438082/212‑1 (MQ=255)
                                                                                                                                                                                                           gCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATTGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGAATTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGGTCGCCtccgaagctctgatgtgggggg                                                     >  1:77087/1‑230 (MQ=255)
                                                                                                                                                                                                                     ggtgTGCGCCGCGAAAAGCTATCTCACACAATACGGCAGACCGGAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGTGCCGCTGATGTGGGTGATCAACgaga                                            <  1:294993/251‑1 (MQ=255)
                                                                                                                                                                                                                          tgCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCCACGTGCCGCTGATGGGGGGGATCAACGAGATCa                                         >  1:79404/1‑250 (MQ=255)
                                                                                                                                                                                                                            ggccgcAAAAAGCTATCCCCCACAATCGGGCAAACCGAAAAACCCCGCCGATTTGAGTATTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGTGCCGCTGATGTGGGTGATCAACGAGATCAATc                                      <  2:103128/250‑1 (MQ=255)
                                                                                                                                                                                                                                 cgAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGACGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCCGGGTCGCCTCCGTGCCGCTGCTGTGGGGGTTCAACCAGATCAATTGTggg                                 >  2:131499/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaaGCTATCTCACACAATACGGAATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGTGCCGCTGATGTGGGTGATCAACGAGATCAATCGTGGGGAg                              <  1:380520/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaaaaccgccttccGGAAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGTGCCGCTGATGTGGGTGATCAACGAGATCAATCGTGGGGAGCTGGAGATCctgc                 <  2:400477/239‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaaTACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGTGCCGCTGATGTGGGTGATCAACGAGATCAATCGTGGGGAGCTGGAGATCctgct                <  1:228710/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGGGCCGCCGATGTGGGGGATCAACGAGATCAATCGTGGGGAGCTGGAGATCCtgctgc              >  1:480038/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      aTACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGCGCCGGGATCGCCTACGGGCCGCTGGTGTGGGGGGTCAACGAGATCAATCGTGGGGGGCTgggggtgctcctgc              >  1:456213/1‑239 (MQ=255)
                                                                                                                                                                                                                                                         ggcaaaaCCGAAAAACCCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGTGCCGCTGATGTGGGTGATCAACGAGATCAATCGTGGGGAGCTGGAGATCCTGCTGCcgc           <  2:456434/245‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGTGCCGCTGATGTGGGGGATCAACAAGATCAAACGTGGGGGGCTGGAGATCCTGCTGCcgcg          >  1:483108/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            cATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGGGCCGGGATCGCCTACGTGCCGCTGATGTGGGGGATCAACGAGATCAATCGTGGGGAGCTGGAGATCCTGCTGCCGCGtt        >  1:457771/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                  ggAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCCACGTGCCGCTGATGTGGGTGATCAACGAGATCAATCGTGGGGGGCTGGGGATCCTGCTGCCGCGTTACCAGt  >  1:80244/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                           
CATTGGCTGTTCTTCAACTATGGCACAAAATGTTCTCGCCGGGCTGACAGCCAAAATGCTGAAAGAATACCCAGGTTTGAGCGTCAATCTGGTTACCGGAATTCCAGCCCCCGACCTGATTGCCGACGGTCTGGATGTGGTGATCCGCGTCGGCGCGTTGCAGGATTCCAGCCTGTTTTCCCGCCGTCTGGGCGCGATGCCAATGGTGGTGTGCGCCGCGAAAAGCTATCTCACACAATACGGCATACCGGAAAAACCCGCCGATTTGAGTAGTCATTCATGGCTTGAATACAGCGTGCGGCCCGACAATGAATTTGAACTGATCGCACCGGAAGGGATCTCGACTCGCCTGATCCCACAAGGAAGATTTGTGACTAATGATCCGATGACGCTGGTGCGCTGGCTGACGGCGGGTGCCGGGATCGCCTACGTGCCGCTGATGTGGGTGATCAACGAGATCAATCGTGGGGAGCTGGAGATCCTGCTGCCGCGTTACCAGT  >  NC_000913/3389804‑3390303

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: