Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,321,019 C→A 6.8% E199* (GAG→TAG)  yhbX ← putative EptAB family phosphoethanolamine transferase, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,321,0190CA 6.8% 198.3 / 5.5 59E199* (GAG→TAG) yhbXputative EptAB family phosphoethanolamine transferase, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base A (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

TCGGGGTAATTATGAATGTCATGACTCAAAACAGAATCAGCAGTTAAAGAAAGGGGAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGCACATAGATAAGGGATATACATCCCCAGCATTTTGA  >  NC_000913/3320733‑3321297
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
tCGGGGTAATTATGAATGTCATGACTCAAAACAGAATCAGCAGTTAAAGAAAGGGGAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCtt                                                                                                                                                                                                                                                                                      >  1:118846/1‑289 (MQ=255)
tCGGGGTAATTATGAATGTCATGACTCAAAACAGAATCAGCAGTTAAAGAAAGGGGAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCtt                                                                                                                                                                                                                                                                                      >  2:373431/1‑289 (MQ=255)
                   cATGACTCAAAACAGAATCAGCAGTTAAAGAAAGGGGAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATa                                                                                                                                                                                                                                                                   >  1:79017/1‑289 (MQ=255)
                                  aaTCAGCAGTTAAAGAAAGGGGAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGaaaaa                                                                                                                                                                                                                                                    <  1:373431/289‑1 (MQ=255)
                                   aTCAGCAGTTAAAGAAAGGGGAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGaaaaa                                                                                                                                                                                                                                                    <  1:520740/288‑1 (MQ=255)
                                   aTCAGCAGTTAAAGAAAGGGGAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGaaaaa                                                                                                                                                                                                                                                    <  1:453692/288‑1 (MQ=255)
                                   aTCAGCAGTTAAAGAAAGGGGAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGaaaa                                                                                                                                                                                                                                                     <  2:626745/287‑1 (MQ=255)
                                                   agGTGAACCGATAGCGCAGTGTAAATTGAGCCGCTATTTGCTTGATGAATCAGTTTGGTCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGTGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAGTCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGtat                                                                                                                                                                                                                                            <  2:75674/280‑1 (MQ=255)
                                                     gggAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTCATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGCAGCaaa                                                                                                                                                                                                                                  >  2:79890/1‑288 (MQ=255)
                                                            gacagcgaaagtgaaGATGCGCAGCTTATTACTTGATTAAGCAGTTTGTTCTGTTTTCTTTGTGCTTCCACTTGCGGTGTCTGAGAGCGTGTATATCCATACAAAGACATATTGTCTACACGTACAGACTCCCCCACACTCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAATCGGCGTATAGGTAGCAAATCGCGACg                                                                                                                                                                                                                          <  2:112775/276‑1 (MQ=255)
                                                            gACAACGCAGTGTAAGGTGAGCCGCTTATTGCGTGATTAAACATGTTGATCTGTTTTCTTTGTGGTTCAACTTGCGGTGTCGTATAGCGGGTATATCCATACAAAGACATATTGTAGACTCGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCATGTGTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAGATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACg                                                                                                                                                                                                                          <  2:592511/289‑1 (MQ=255)
                                                             aCAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCaaa                                                                                                                                                                                                                                  >  1:75674/1‑280 (MQ=255)
                                                                 cgcAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAg                                                                                                                                                                                                                      <  2:575359/288‑1 (MQ=255)
                                                                 cgcAGTGGAGGGTGCGCGGCTGATTGCTTGGTTAAACAGTTTGTCCCGTTTTCTTGGTGCTTCAACTTGCGGGGTCGTAGAGCGGGTATAGCCATGCAAAGACATATAGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATGAATACCTGTATGCCTGACAGATAATTGAAAATACGGCACCGTGGTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAg                                                                                                                                                                                                                      <  2:245205/288‑1 (MQ=255)
                                                                                gccgcTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACg                                                                                                                                                                                                                          >  1:343526/1‑269 (MQ=255)
                                                                                gccgcTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACg                                                                                                                                                                                                                          <  2:343526/269‑1 (MQ=255)
                                                                                gccgcTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAAc                                                                                                                                                                                                      >  2:7294/1‑289 (MQ=255)
                                                                                gccgcTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTATTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAAc                                                                                                                                                                                                      >  2:250302/1‑289 (MQ=255)
                                                                                                  aaaCAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCa                                                                                                                                                                                    <  2:79017/289‑1 (MQ=255)
                                                                                                         ttGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGtat                                                                                                                                                                                                                                            <  2:409352/226‑1 (MQ=255)
                                                                                                         ttGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGtat                                                                                                                                                                                                                                            >  1:409352/1‑226 (MQ=255)
                                                                                                          tGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGaaaaa                                                                                                                                                                                                                                                    >  1:445088/1‑217 (MQ=255)
                                                                                                          tGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGaaaaa                                                                                                                                                                                                                                                    <  2:445088/217‑1 (MQ=255)
                                                                                                                     cTTTTTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGAc                                                                                                                                                                 <  2:465490/289‑1 (MQ=255)
                                                                                                                      tttGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCATTACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACa                                                                                                                                                                >  2:141547/1‑289 (MQ=255)
                                                                                                                         gtgCTTCAACGTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGc                                                                                                                                                             <  2:161941/289‑1 (MQ=255)
                                                                                                                           gCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGaaaaa                                                                                                                                                                                                                                                    <  1:287155/200‑1 (MQ=255)
                                                                                                                           gCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGaaaaa                                                                                                                                                                                                                                                    >  2:287155/1‑200 (MQ=255)
                                                                                                                                                   cGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACaa                                                                                                                                   >  1:428949/1‑289 (MQ=255)
                                                                                                                                                     tgtATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATc                                                                                                                                 >  2:353863/1‑289 (MQ=255)
                                                                                                                                                     tgtATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTACATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATc                                                                                                                                 >  2:658888/1‑289 (MQ=255)
                                                                                                                                                        atatCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTATTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAg                                                                                                                               <  1:250302/288‑1 (MQ=255)
                                                                                                                                                                   aaGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACg                                                                                                                                                                                                                          >  2:352308/1‑186 (MQ=255)
                                                                                                                                                                   aaGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACg                                                                                                                                                                                                                          <  1:352308/186‑1 (MQ=255)
                                                                                                                                                                                gACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACttttttt                                                                                                      >  2:539694/1‑289 (MQ=255)
                                                                                                                                                                                       aCAGACTCCCCCACAATGAACACGTAGGGATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCGGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCGAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCATTTTTTTTGTCGGCa                                                                                               <  2:308975/289‑1 (MQ=255)
                                                                                                                                                                                               ccccACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAg                                                                                                                                                                                                                      >  1:439695/1‑162 (MQ=255)
                                                                                                                                                                                               ccccACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAg                                                                                                                                                                                                                      <  2:439695/162‑1 (MQ=255)
                                                                                                                                                                                                               tAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGATCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTAATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCATTAATATTCCTGTACCTTTTTTTGTCGGCAAGGAAAAATCATATTttattatta                                                                       >  2:671543/1‑289 (MQ=255)
                                                                                                                                                                                                                 ggTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACt                                                                     >  2:600610/1‑289 (MQ=255)
                                                                                                                                                                                                                   tATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTg                                                                    >  1:296723/1‑288 (MQ=255)
                                                                                                                                                                                                                     tCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCa                                                                 >  1:493948/1‑289 (MQ=255)
                                                                                                                                                                                                                            cTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGaaa                                                                                         >  2:436945/1‑258 (MQ=255)
                                                                                                                                                                                                                            cTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGaaa                                                                                         <  1:436945/258‑1 (MQ=255)
                                                                                                                                                                                                                               tATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAaagaag                                                       >  2:260678/1‑289 (MQ=255)
                                                                                                                                                                                                                                      gactgaTAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATaaa                                                >  1:367765/1‑289 (MQ=255)
                                                                                                                                                                                                                                       actgaTAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTATTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGa                                                                                           >  2:503125/1‑245 (MQ=255)
                                                                                                                                                                                                                                       actgaTAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTATTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGa                                                                                           <  1:503125/245‑1 (MQ=255)
                                                                                                                                                                                                                                           aTAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGc                                           >  2:438758/1‑289 (MQ=255)
                                                                                                                                                                                                                                                   aaaTACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGcac                                   <  1:260678/289‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaTACGGGACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGcaca                                  <  2:24629/289‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaTACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGcaca                                  <  2:118846/289‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaTACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGcaca                                  <  1:600610/289‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGCACATAGATAAg                          <  1:438758/289‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGCACATAGATAAg                          <  1:141547/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                tgtTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGCACATAGATAAGGGat                      >  1:194704/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                       aaTTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAgg                                                                                            <  2:62858/212‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aaTTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAgg                                                                                            >  1:62858/1‑212 (MQ=255)
                                                                                                                                                                                                                                                                                    cTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGCACATAGATAAGGGATATACATCCCCAGCATTTTGa  <  1:7294/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
TCGGGGTAATTATGAATGTCATGACTCAAAACAGAATCAGCAGTTAAAGAAAGGGGAACCGACAGCGCAGTGTAAGGTGCGCCGCTTATTGCTTGATTAAACAGTTTGATCTGTTTTCTTTGTGCTTCAACTTGCGGTGTCGTAGAGCGTGTATATCCATACAAAGACATATTGTCGACACGTACAGACTCCCCCACAATCAACACGTAGGTATCAATACCTGTATCCCTGACTGATAATTGAAAATACGGCACCGTGTTTGCAATTGAGAGTAATCTTTGATGCTCTTTCGCTGCTAAAGCAAAATAGTTGAGATTGAAAAACGGCGTATAGGTAGCAAATCGCGACGCTAGTATATATGGACTGAACGCTTTTTTATTTTTTGCATCTTTATAAGCAAATTGACAAGCGGAAAATAAACTGCCCGAAATGACAATCAGCAATAATATTCCTGTCACTTTTTTTGTCGGCAAGGAAACATCATATTTTATTATTACTGCCAAAAAAAGAAGGGATAAAAAAGCAAAGGCACATAGATAAGGGATATACATCCCCAGCATTTTGA  >  NC_000913/3320733‑3321297

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