Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,665,7020CA27.7% 196.5 / 31.0 83A37S (GCG→TCG) gadXacid resistance regulon transcriptional activator; autoactivator
Reads supporting (aligned to +/- strand):  ref base C (55/5);  new base A (0/23);  total (55/28)
Fisher's exact test for biased strand distribution p-value = 5.36e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAA  >  NC_000913/3665476‑3665912
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                    
aGCATGTCTGCATAGCGAAGTTTTTTCTTCTGCGTTACTTGCGCATCTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCGCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCa                                                                                                                                                                                                                     <  2:1930527/228‑1 (MQ=255)
  cATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaa                                                                                                                                                                                                                   >  1:942919/1‑228 (MQ=255)
          cATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGACCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACAAGGTc                                                                                                                                                                                                        >  1:1440665/1‑229 (MQ=255)
             aagagagtttgttATTCGGCATTGGTTGCAAAACTTCAGAAAAAAAAGTAATACCTCTACAACAGTGAGGAACGATAAAAGAAGAACTTTTTTAATGTGCAGGAAGATTTACTCCCCCCTCCTGGAGCATCGGGAAAAATAATGAAAGCGTGTCCCTTGACACAAAAAAAAAATTTCCAACCCATTTATCAACCAAATTTTGGCTTGCACCCCAAAAACCAAGGTCAcc                                                                                                                                                                                                     <  2:2576/218‑1 (MQ=255)
               cGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGc                                                                                                                                                                                                   >  1:1768404/1‑229 (MQ=255)
               cGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTAACCGc                                                                                                                                                                                                   >  1:421862/1‑229 (MQ=255)
                              tGCGATTGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAACAGGTCACCGCCATTAAAATAGCGat                                                                                                                                                                                    >  1:348745/1‑229 (MQ=255)
                                cGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCCTTAAAATAGCGata                                                                                                                                                                                   >  1:1677586/1‑228 (MQ=255)
                                cGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGata                                                                                                                                                                                   >  1:492091/1‑228 (MQ=255)
                                cGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGata                                                                                                                                                                                   >  1:1716427/1‑228 (MQ=255)
                                  atatTTGCGTAACTCGCTGAAAACCAGAAATGTCTCTTCTAAGGTTAAGAACGAATATAGAAAATATTTTTTAATGTGTAGAAAGATTTATGTCCACCCCCTTGATCATCGGAAGAATAAAGCAAAGGGTTCCCTTGCACAAAAACCAAAAATTCTCAACACATTTATTTATTTAATTTTGTCTTCCACCCAAAAAACCCGGGTCCCCGCCATTAAAATAGCGATATTc                                                                                                                                                                                <  2:1185709/229‑1 (MQ=255)
                                           cAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGACACCGCCATTAAAATAGCGATATTCACCATTAAc                                                                                                                                                                       >  1:834329/1‑229 (MQ=255)
                                             aCTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTAATCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCa                                                                                                                                                                     >  1:843340/1‑229 (MQ=255)
                                                 tataaaaaaCAGAAATATCTTTAATAAAGTGAAGAAAAAAAAAGAAAAAATCTTTTTTATGGTAATGAATTTTTATGCCCCCCTCCTTGGGCACCGGAAAAAAAAATAAAAGCTTTCCCTTTAACAAGAATAAAAAATTCCAAACACACTTATCTACTCAATTTTTGCTTCCACCAGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGt                                                                                                                                                                 <  2:1438196/225‑1 (MQ=255)
                                                 cagaaaaaaCAGAAATTCCCATACTACATTGATCAACTAGTAAGAAAAAATCTTTTTTATGTACATAAATTTTTACTCCCTCCCCCTAGAGCATGGGAAAAAAAAAAGAAGCCGTCCCCCTTGCCACGAAAAAAAAATTTCCCACCCCCTTATCTCCTCAATTTTGCTTTCCCCCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATgg                                                                                                                                                                  <  2:421109/224‑1 (MQ=255)
                                                 agtaaaaaaCAGAAATATCTCACCGAAAGTGATGAACGATTAAGAAAGAAATTTTTTATGTTGAAGGAAGATTGAATGCCTCCTCATGGAGCATCGGAAAAAATAATGAAAGCGTTTCCCTTGACACGAATACAAAATTTCCACCACACTTTTCTACTTGATTTTGTTTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGt                                                                                                                                                                 <  2:216872/225‑1 (MQ=255)
                                                                           caaggatgaATGAGTTATGAAAAAACTTTTTTTTGGGGAATAAATTTTTATTCCCTCTCCTTGGGAAGCGGGAAAAAATAAGAAAGGCGGCTCCTTGGCCACAAAAAAAAATTTTAAACCCCCCTTATCTATTAAAATTTGTTTTCCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATa                                                                                                                                         <  2:1544970/224‑1 (MQ=255)
                                                                            aGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTTGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGc                                                                                                                                      >  2:800321/1‑229 (MQ=255)
                                                                              tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCa                                                                                                                                                                                                                     >  1:1832442/1‑150 (MQ=255)
                                                                              tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCa                                                                                                                                                                                                                     <  2:1829327/150‑1 (MQ=255)
                                                                              tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCAcc                                                                                                                                                                                                     >  1:1276682/1‑164 (MQ=255)
                                                                              tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATAtt                                                                                                                                                                                 >  1:519160/1‑184 (MQ=255)
                                                                              tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCa                                                                                                                                     >  2:872463/1‑228 (MQ=255)
                                                                              tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCa                                                                                                                                     >  1:527357/1‑228 (MQ=255)
                                                                              tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCa                                                                                                                                     >  1:54682/1‑228 (MQ=255)
                                                                              tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCa                                                                                                                                     >  1:115530/1‑228 (MQ=255)
                                                                              tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTAATCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCa                                                                                                                                                                            >  1:1319447/1‑189 (MQ=255)
                                                                              tgatgaACGAGTAAGGAAAAAACTTTTTAATGTGAATGAAGATTTAATGCCCCCTCCTTGAGCATGGGGAAAAAAAATGAAGGCGTGCCCCTTGACAAAAAAAAAAAATTTCCAACACCCTTTCCTCCTCAATTTTGGCTTCCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATAtt                                                                                                                                                                                 <  2:518237/184‑1 (MQ=255)
                                                                                atgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACACAATTTTCAACACGCTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAt                                                                                                                                   >  1:1767646/1‑228 (MQ=255)
                                                                                atgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATa                                                                                                                                         >  1:1547637/1‑222 (MQ=255)
                                                                                atgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAt                                                                                                                                   >  1:784033/1‑228 (MQ=255)
                                                                                atgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAt                                                                                                                                   >  1:689692/1‑228 (MQ=255)
                                                                                atgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAt                                                                                                                                   >  1:1132430/1‑228 (MQ=255)
                                                                                   aaCGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGat                                                                                                                                                                                    >  1:1479393/1‑176 (MQ=255)
                                                                                   aaCGAATAATGAAGAACTTTTTTTAGGGGAATGAAATTTTAATCCCTCCTCCTGAGGCAGCGGAAGAAAAAATAAAAGCGTGCCCCTTAGCAAAAAAAAAAAATTTCCAACACTCTTATCTACTCGAATTTGGTTTCCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGa                                                                                                                               <  2:1895777/229‑1 (MQ=255)
                                                                                     cGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATAtt                                                                                                                                                                                 >  1:682557/1‑177 (MQ=255)
                                                                                      gagaAATGAAAAATTTTTTTAAGGTACAAAAAATTTTAATCCCTCCTCCTTGGGCATGGGAAAAAAAAAAGAAAGCGTCTCCCTTAACAGAAAAACAAAATTACAAACACTTTTATCTATTGGATTTTGGTTTGCTTCCAAAAAACCAAGGTATCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa                                                                                                                            <  2:1130459/229‑1 (MQ=255)
                                                                                      gaatAAAGAAAATACTTTATGAGGTGCATGAATATTTAACCCCTCCTACATGATCATCGGAAAAAATAAAGAAGGCTGGTCCCTTGACAAGAAAAAAAAATTTTCAACACTCTTATCACCAAAATTTTGGTTTCCACCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa                                                                                                                            <  2:1765400/226‑1 (MQ=255)
                                                                                        gTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAAtt                                                                                                                          >  1:1823456/1‑229 (MQ=255)
                                                                                           aaaaaaaattttttttGTGTGAAAGAATATTTAATCCCCCCTCCTGGAGAGGGGGGAAAAATAATGAAGGCGTTCCCTTTAACAAAAAAAAAAATTTTAAACCCCCTTTATCTCCTAAATTTGGGTTTCCACCCAAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTcc                                                                                                                        <  2:1746354/218‑1 (MQ=255)
                                                                                                gaaACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATAtt                                                                                                                                                                                 >  1:29953/1‑166 (MQ=255)
                                                                                                 acttttttttATGGATAGGAAGTTTAATTCCCCCCTCCTGAAGAGTGGGGAAAAAAAAAGAACGCGTGTCCCTTAACAAGAAAAAAAATTTTCAAACCCACTTATTTATACGTTTTTGGTTTCCATCCAAAAAAACGAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa                                                                                                                            <  2:688483/214‑1 (MQ=255)
                                                                                                  ctcTTTTTTATGTAAATAAATTTTGTATCCCCCCTCCGTGACCAGGGGAAAAAAAAAGAAAAGCGTCTCCTTTGACAGAAATAAAAAATTTCAAACCCCTTTATCAATGAGATTTTTGTTTGCCCCCAAAAAAACCGGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                                                                                                 <  2:1764642/226‑1 (MQ=255)
                                                                                                              gacaaTGAATATTTATCGCCCCCTCCTGAAGCATGGGGAAAAATAATAAAGGCGTTTCCCTAAAAACAAAAAAAAAATTTAAAACACACTTATCGCCTGAATTTTGGTTACCACCCAAAAAACCAAGGAAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                                                                                                     <  2:1722271/225‑1 (MQ=255)
                                                                                                                   tGAATATTAAATCCCTCCTCCTTGACAATGGGAAAAATAAATGAAAGCTTTCCCCTTGACACAAAAAAAAAATTTCAAACCCATTTTTTTAATCGATTTGGGCTTCCACCAAAAAAAACCGTCTCTCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACt                                                                                                <  2:491223/229‑1 (MQ=255)
                                                                                                                                cctcctcctTGAGAATCGGAAAAAATAATGAACGCGTGCCCCTTGACAAGAATAAAAAATTTTCAACCCTTTTATCTACTCGATTTTGGTTTACATCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCa                                                                                                         <  2:782664/207‑1 (MQ=255)
                                                                                                                                    ctcctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCAcc                                                                                                                                                                                                     >  1:208165/1‑111 (MQ=255)
                                                                                                                                    ctcctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAg                                                                                                                                                                                        >  1:1969858/1‑124 (MQ=255)
                                                                                                                                           aGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATAtt                                                                            >  1:1012853/1‑225 (MQ=255)
                                                                                                                                            acaacaggaaaaaaaaaTAAAAGTGTGCCCCTACACACAAAAAAAAATTTCCAAACACATTTATCAACTCAATTTGGCTTTCCACCCAAAAAAACGAGTCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTg                                                                        <  2:373080/218‑1 (MQ=255)
                                                                                                                                               tagggAAAAATAATGAAGGCGTGTCCCTTGACACAAATAAAAAATTTTCAACCCATTTATCTACTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                                                                                                   <  2:343444/180‑1 (MQ=255)
                                                                                                                                               tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                                                                                                   >  1:344066/1‑182 (MQ=255)
                                                                                                                                                cGGGCGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTa                                                                           >  1:1501746/1‑221 (MQ=255)
                                                                                                                                                cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAcc                                                                                                                                                                                                              <  1:1572340/91‑1 (MQ=255)
                                                                                                                                                cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAcc                                                                                                                                                                                                              >  2:1569643/1‑91 (MQ=255)
                                                                                                                                                cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                                                                                                                                            >  1:1886589/1‑140 (MQ=255)
                                                                                                                                                cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATga                                                                                                                   >  1:755049/1‑180 (MQ=255)
                                                                                                                                                cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATAtt                                                                                         >  1:403075/1‑207 (MQ=255)
                                                                                                                                                cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTg                                                                        >  1:373754/1‑224 (MQ=255)
                                                                                                                                                 gggaaaaaaTAAAGAAGGCGTGTCCTTGACAAAAAAAAAAAAATTTCCAACCCCATTTCCTCCACAATTTTGGTTGCCCTCCAAAAAAACCAGGTAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATAtt                                                                                         <  2:402350/201‑1 (MQ=255)
                                                                                                                                                 gggAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATAtt                                                                                                                                                                                 >  1:1291507/1‑118 (MQ=255)
                                                                                                                                                 agaaaaaaaTAATGAAAACGCGTCTATAACAACGAAAAAAAAATTTCAACCACACTTACCATTACGATGTTGACTTGCATCGAAAAAAACCAGGACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTa                                                                           <  2:1499153/215‑1 (MQ=255)
                                                                                                                                                                  aGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCAt                                                                                                                                                                           >  1:641591/1‑109 (MQ=255)
                                                                                                                                                                        tCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCtt                                             >  1:799840/1‑229 (MQ=255)
                                                                                                                                                                                      aaTACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACt                                                                                                >  1:306929/1‑166 (MQ=255)
                                                                                                                                                                                                                     cTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCAcgccg                                                    >  1:61707/1‑182 (MQ=255)
                                                                                                                                                                                                                      tCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGc                                                                                                                                                                                                   >  1:1151148/1‑38 (MQ=255)
                                                                                                                                                                                                                      tCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATAtt                                                                                                                                                                                 >  1:508363/1‑56 (MQ=255)
                                                                                                                                                                                                                      tCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                                                                                                                                                             >  2:1963044/1‑60 (MQ=255)
                                                                                                                                                                                                                      tCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGgag    >  2:1642100/1‑229 (MQ=255)
                                                                                                                                                                                                                      tCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGgag    >  2:1292254/1‑229 (MQ=255)
                                                                                                                                                                                                                      tCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTAATTGATAAATTAATTTAAACGAAACGGgag    >  2:1071465/1‑229 (MQ=255)
                                                                                                                                                                                                                        gattttGTCTTGCACCAAAAAAAACAAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATa                                                                                                                                         <  2:1161796/91‑1 (MQ=255)
                                                                                                                                                                                                                        gattttGGCTTGCACCAGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                                                                                                     <  2:552402/127‑1 (MQ=255)
                                                                                                                                                                                                                        gAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATAtt                                                                                                                                                                                 >  1:1804953/1‑54 (MQ=255)
                                                                                                                                                                                                                        gAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATa                                                                                                                                         >  1:1163824/1‑94 (MQ=255)
                                                                                                                                                                                                                        gAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                                                                                                     >  1:1321630/1‑130 (MQ=255)
                                                                                                                                                                                                                        gAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                                                                                                     >  1:553375/1‑130 (MQ=255)
                                                                                                                                                                                                                         gattttgcatcgCACCCACAAAAACCAGACTCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGaa  <  2:1820355/218‑1 (MQ=255)
                                                                                                                                                                                                                               gCTTGCCCCCGAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTccc                                                                                                                       <  2:1422647/105‑1 (MQ=255)
                                                                                                                                                                                                                               gCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTccc                                                                                                                       >  1:1425091/1‑105 (MQ=255)
                                                                                                                                                                                                                                     cATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGt                                                                                                                >  2:1428717/1‑107 (MQ=255)
                                                                                                                                                                                                                                     cATCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGt                                                                                                                <  1:1431171/107‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                    
AGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAA  >  NC_000913/3665476‑3665912

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