Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913954,3880CT25.0% 79.3 / 12.2 57A27T (GCA→ACA) focAformate channel
Reads supporting (aligned to +/- strand):  ref base C (36/6);  new base T (1/13);  total (38/19)
Fisher's exact test for biased strand distribution p-value = 1.75e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGA  >  NC_000913/954153‑954605
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                           
aTGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGGAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCttt                                                                                                                                                                                                                          >  1:1754521/1‑239 (MQ=255)
 tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATAGCCAGATAGAAAGTCTTAAGCGGAAGTTTTGTTGCTTta                                                                                                                                                                                                                         >  1:1969021/1‑239 (MQ=255)
 tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGAGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCCGATAGAAAAACTTAAACGGATGTTTCGTAGCTTta                                                                                                                                                                                                                         >  1:193864/1‑239 (MQ=255)
  gcccccACTCGACTTAGCAACAACAACCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGTTAGAAAGTCTTAAGCGGATGTGTCGTTTCTTtat                                                                                                                                                                                                                        <  2:946491/236‑1 (MQ=255)
        aCTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGAGAGCAAGTCTTAAGCGGGTGTGTCCGTGCTTTATAGaca                                                                                                                                                                                                                   >  2:615952/1‑238 (MQ=255)
                                 acGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAAACTTAAGCGGATGTTTCGTTGCTTTATAGACACCCCCCCCTTCGGTCCCTTTGGCCAt                                                                                                                                                                                        >  1:1428089/1‑239 (MQ=255)
                                    gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTTTAGACACCCGCCTCTTCGGCCCCTTTGGCCATTgc                                                                                                                                                                                     >  1:1547935/1‑239 (MQ=255)
                                    gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGAAAGAAAGGATTAAGCGGATGTTTCGGTGATTTATAGACACCCGCCTCCTCGGGCACCTTTGGCAAtgg                                                                                                                                                                                      >  1:1391442/1‑237 (MQ=255)
                                     tGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACCTTGGCCATTgca                                                                                                                                                                                    >  1:4783/1‑239 (MQ=255)
                                              agagATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAACGGATGTTTCGTTGCTTTTTAGACACCCGCCCCTTCGGCCCCTTTTGCCATTGCAGCAGGaaag                                                                                                                                                                          >  1:1635893/1‑237 (MQ=255)
                                               gagaTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAACGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTa                                                                                                                                                                         >  1:842849/1‑239 (MQ=255)
                                                  aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCCCCTCCTCCGCCACTTTGGCCATTGCAGCAGGAAGGaaaa                                                                                                                                                                      >  1:425636/1‑239 (MQ=255)
                                                       cTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCa                                                                                                                                                                 >  1:1781863/1‑239 (MQ=255)
                                                       cTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAAACGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCCCCTCCTTGGCCCCCTTGGCCCATGCCGCAGGAAGTAAAACAtac                                                                                                                                                                 >  2:562359/1‑237 (MQ=255)
                                                          cGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAAACGCCAGGAAGGAAAGCTTAAGCGGGTGGTTTGTTGCTTTATAAGAACCCGCCCCTTTGGCCACTTTGGCTCTTGCAGGAGGGAAGAAAAAGACaaaa                                                                                                                                                              >  2:22162/1‑239 (MQ=255)
                                                                  acaaAAAGATCCAGCCCCAGAGAGAAGAAAGTGCCCCAACCCTGTTTTCCAGTGCCGAAGTGTTTGGTCCCTGTCCCAGTGTGTCCTGTAATAAAGAAGACGAATCCGTTTGAGTTAAAACCACCGGCGGTAACCGCTAGATAAAAAGTCTTAAGCGGATGTTTCGTTGTTTTATATACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTc                                                                                                                                                      <  2:1428089/239‑1 (MQ=255)
                                                                   cacaAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCt                                                                                                                                                                                                                            <  1:236217/170‑1 (MQ=255)
                                                                      aaaGAATCCGGCCGAGGGCAGAACAAAAGGCCGCCAACAAGTTGGACCAGCGCGAGGGGCTGTGGGCCGGTGCCGGTGGTGTGTGTGAAAAAGACGGAGAGTGGGGTGGAGATGAAAACCCGGGGGGTAATCGCAAGAAAGAAAGTCTAAAGCGGTTGTTTGGTTGTTTTATAGACACCCGCCTCTGCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGc                                                                                                                                                   <  2:1554164/239‑1 (MQ=255)
                                                                          aaTCACCCCGAGAGAGAAGAAAGTGCCGCGACCCATTTTGCCCAGGCCGAAGGGCTTTGTGCCTGGGCAAGTGGTTGCTGTAAAATAGAAGACAAATGCAATTGAAAGGAAACCCCCGGCGGTAATCCCAAGATAGAAAGTCTTAAGGGGATGTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCttt                                                                                                                                                <  2:4783/237‑1 (MQ=255)
                                                                            tCAGCCCCAGATAGAATCAAATGCCGCCACCCAGTTTAGCCATGCAGAAGGGCTTTGGCCCTGTGCGAGTGTGTGTGGTAAAATAGAAGAGAAATGCGATGAAGATGAAACCACCGGTGGAAATCCCAAGAAAAAAAGCCTTAGGTGGATTTTTGGTTGTTTTAGAGCCGCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTcaca                                                                                                                                            <  2:1591431/239‑1 (MQ=255)
                                                                             cAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGCTGAAAACACCGGCGGGAAACGCCAGGTAGAAAAGCTTAAGCGGAGGTTTCGTTGCTTTATATacac                                                                                                                                                                                                                  >  2:236217/1‑170 (MQ=255)
                                                                                ccccAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTa                                                                                                                                         >  1:1591431/1‑238 (MQ=255)
                                                                                ccccAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGGAAAAGATCAAAAGGGTTGTCAGCTTTCACACTa                                                                                                                                         >  1:637926/1‑238 (MQ=255)
                                                                                   acaagagcagacAAAAGCCACAAACCAGTTTTGCCCTGCGAGAGGGTTTTGGGCGGGTGCAGGGGGTGGCGGGGAGAAAGAAGAAGAATGCAATGGAGAAGACACCCCGGGAGTAAACCGCAAAAAAAAAAGTCTAAAGCGGAGTTTTGTTTGTTTTATAGACCACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAAct                                                                                                                                      <  2:2129506/228‑1 (MQ=255)
                                                                                    agagagAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCt                                                                                                                                                                                                                            <  1:1802662/153‑1 (MQ=255)
                                                                                            aGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTTCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTa                                                                                                                                         >  1:1438812/1‑227 (MQ=255)
                                                                                               aaaTGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGGTCTAAGCGGATGGTTCGTTGCCTTATAGacac                                                                                                                                                                                                                  >  2:1802662/1‑153 (MQ=255)
                                                                                                                     aaggacGAGGGACATTGTGCTGGGCCCAGTGGTTGTGGTAATAAAAAAGAAGAATGCGATTGAGAAAAACACCCCGGCGGTAACCGCAAGATAAAATGTCTTAGGGGTTTGTTTCGTTGTTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                                    <  2:425636/234‑1 (MQ=255)
                                                                                                                         ccGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGGAAAAGATCAAAAGGGTTGGCAGCTTTCACACTAACTCTCTCTTTATTAAg                                                                                                                        >  1:1685524/1‑215 (MQ=255)
                                                                                                                         ccGAAGGGCATTGTGCCTGTGCCAGCGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCCGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCCTTGCAGCACGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACCCTCTCTTTATTAAGCCGGCGCCGCGATACTAACaaaac                                                                                                >  2:1254309/1‑237 (MQ=255)
                                                                                                                                 cATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAACATTATAGa                                                                                        >  1:1372162/1‑239 (MQ=255)
                                                                                                                                           gtgccAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTttatt                                                                                                                           >  1:171247/1‑194 (MQ=255)
                                                                                                                                           gtgccAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTg                                                                              >  1:2104065/1‑239 (MQ=255)
                                                                                                                                           gtgccAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAATCTTAAGCGGATGTTTCGTTGCTTTATAAACACCCGCCTCTTCGGCCACTTTTGCCATTTCAGCAGGAAGTAAAAGATCAAAAAGGTTGTCAGCTTTCACACTAACTCTCTCTTTTTTTAGTCGGCGCCCAGAAACTAACAAAACATTAACGATGGGAAATTg                                                                              >  2:347815/1‑239 (MQ=255)
                                                                                                                                              ccAGTGGTTGCTGTAATATAAAAGACGAAGGCGTTTGAGATAAAAACCCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCTTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGat                                                                            <  1:22162/238‑1 (MQ=255)
                                                                                                                                                aGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTTGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCa                                                                                                                                                     >  1:100461/1‑163 (MQ=255)
                                                                                                                                                aGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAAATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGAtata                                                                         >  1:1941365/1‑239 (MQ=255)
                                                                                                                                                  tGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGGTTCGTTGCTTTATAGACACCCCCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAATTCGGCGACGAGATACTAACAAAGCATtata                                                                                          >  1:1831921/1‑220 (MQ=255)
                                                                                                                                                           tgATATAGAAGACGAATGCAATGGAGATAAAACCCCCGGGGGTAATCGCAAGAAAGAAAGTTTTAAGGGGGTGTTTCTTTTTTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgcagca                                                                                                                                                                                 <  1:654274/125‑1 (MQ=255)
                                                                                                                                                             atatAGAAGACGAAGGCGATTGATAAGAAACCCCCGGCGGTAACCGCGAAATAGAAATTTTTAAGGGGTGTTTGCGTTGCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGc                                                            <  1:1337100/239‑1 (MQ=255)
                                                                                                                                                             atatAGAAGACAAATGCGAGTGAGAAAAAACCCCCGGGAGAAAGCGACAGATAAAAAGTCTAAGGCGGTTGTTTCTGGGCTTAAAAGACCACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                                <  2:1473696/187‑1 (MQ=255)
                                                                                                                                                                     gaCGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAAGCTTAAGCGGGTGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                                      >  2:654274/1‑125 (MQ=255)
                                                                                                                                                                       cGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCCGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTaa                                                                                                      >  1:1473696/1‑187 (MQ=255)
                                                                                                                                                                         aaggggATTGAAATAAAACCCCGGGCGTTAATCGACAGATAAAAAGTTTAAAGGGGTTGTTTCTTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                           <  1:347815/223‑1 (MQ=255)
                                                                                                                                                                              gATTGAGAGAAAAACCCCGGCGGTAATCGCCAGAGAAAAAGTCTTAAGCGGTTTTTTCTTTGTTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGt                                           <  1:137487/239‑1 (MQ=255)
                                                                                                                                                                              atttGAGATAAACCCCCCGGCGTATACCGCCAGTAAAAAAGCTTAAGGGGGTTGTTGCTTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGt                                           <  1:615952/237‑1 (MQ=255)
                                                                                                                                                                                       gaaaccccgggcggTAATCGCAAGATAGAAAGTCTTAAGGGGATGTTTTGTTGTTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCAt                                  <  1:562359/230‑1 (MQ=255)
                                                                                                                                                                                          aaCACCGGCGGTAATCGCCAGATAGAAAGTTTTAAGCGGATTTTTCGTTGTTTTATTGACACCCGCCTCTTTTGGCACTTTTGCCATTTTAGCAGGAAGTaaaa                                                                                                                                                                      >  1:32514/1‑104 (MQ=255)
                                                                                                                                                                                               cggcggTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAg                                                 >  1:192786/1‑216 (MQ=255)
                                                                                                                                                                                                                       tttaaggCGGATGTTTAGTTGTTTTAGAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAg                                                                                                 <  2:1006162/138‑1 (MQ=255)
                                                                                                                                                                                                                         ttAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTgcagca                                                                                                                                                                                 >  1:1394116/1‑63 (MQ=255)
                                                                                                                                                                                                                            aGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAat        >  1:2041938/1‑228 (MQ=255)
                                                                                                                                                                                                                              cGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAatat      >  1:1472771/1‑228 (MQ=255)
                                                                                                                                                                                                                                 aTGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGAt                                                                                       >  1:1006162/1‑144 (MQ=255)
                                                                                                                                                                                                                                       tCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                                >  2:1072591/1‑114 (MQ=255)
                                                                                                                                                                                                                                       tCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                                <  1:1072591/114‑1 (MQ=255)
                                                                                                                                                                                                                                             gCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGa  >  1:1365783/1‑219 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                           
ATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGA  >  NC_000913/954153‑954605

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