Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,335,654 T→A 9.8% intergenic (‑64/+40) eptA ← / ← adiC lipid A phosphoethanolamine transferase/arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,335,6540TA 9.8% 144.0 / 5.8 51intergenic (‑64/+40)eptA/adiClipid A phosphoethanolamine transferase/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (24/22);  new base A (3/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

GCCAGTGGTCGATTAAGCCATAAGAAAGAGCTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCACAGGTGTAAAGATATGGCACCAGT  >  NC_000913/4335355‑4335949
                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                       
gcCAGTGGTCGATTAAGCCATAAGAAAGAGCTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAACGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAAACTTAAGtt                                                                                                                                                                                                                                                                                                        >  2:135135/1‑301 (MQ=255)
      ggTCGATTAATCCATAAGAAATAGCTTAGTGTCAGTACAATATTAATCACGCTGAAAGCGACGACCGGCAGCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAg                                                                                                                                                                                                                                                                                                  <  1:396490/301‑1 (MQ=255)
           tttAAGCCATAAGAAAGAGCTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGAACTTAAGGTTgg                                                                                                                                                                                                                                                                                             <  2:618019/300‑1 (MQ=255)
             tttcagacgtgtgctcttccgatctCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAACATAAAAGCCAACCTTAAGAACTTAAGGTTgg                                                                                                                                                                                                                                                                                             <  2:538669/277‑1 (MQ=255)
                   aTAAGAAAGAGCTTAGTGTCAGGAAAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATg                                                                                                                                                                                                                                                                                     <  1:211974/301‑1 (MQ=255)
                         aatagCTTGGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAATAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCt                                                                                                                                                                                                                                                                               <  1:370579/298‑1 (MQ=255)
                            agCTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATATTTGCTGGa                                                                                                                                                                                                                                                                            >  2:508584/1‑301 (MQ=255)
                                    agttcagacgtttgctcttccgatctaGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAgggg                                                                                                                                                                                                                                                                    <  2:480152/275‑1 (MQ=255)
                                      cAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAACATAAAAGCCAACCTTAAGAACTTAAGGTTggaaaacggaagagcgtcgggaaggag                                                                                                                                                                                                                                                                    >  1:537996/1‑276 (MQ=255)
                                                 tAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCCAAGCAAA‑AT‑TAAGCCAAACTTAAGTTCTTAAGGTTGGCTTTAATGTTTACTGGATACAAGGGAAATAAACTTTgt                                                                                                                                                                                                                                                       >  1:206173/1‑300 (MQ=255)
                                                           gAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                                                                                                                                                                                                                             <  2:519458/301‑1 (MQ=255)
                                                              aGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaaatcggaagagcgtcgggtagggaa                                                                                                                                                                                                                                          >  1:479551/1‑277 (MQ=255)
                                                                  acgacCGGGATAGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAAAACCTGTTTAAAAAAGGAAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                                                                                                                                                                                                                      <  1:248658/301‑1 (MQ=255)
                                                                           aTCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCAAAGGGATACggg                                                                                                                                                                                                                             <  2:510890/301‑1 (MQ=255)
                                                                                   gAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAtt                                                                                                                                                                                                                     >  2:556651/1‑300 (MQ=255)
                                                                                               ggtATGCAGCTAATCCAGCGGCAGCGCCTGTAACACTTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTa                                                                                                                                                                                                         <  1:555937/299‑1 (MQ=255)
                                                                                                    gCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAACGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTCAGGGGAATTAATATTTGCTTATTGGTGAATATAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTct                                                                                                                                                                                                    >  1:26123/1‑300 (MQ=255)
                                                                                                     cAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAg                                                                                                                                                                                                   <  2:390552/301‑1 (MQ=255)
                                                                                                             ttcagacgtgtgctcttccgatctttttAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACa                                                                                                                                                                                           <  2:533066/276‑1 (MQ=255)
                                                                                                                        gcCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatg                                                                                                                                                                                <  2:26174/301‑1 (MQ=255)
                                                                                                                                     gTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCAAACaagaccggaaaagcgtcgtggaggg                                                                                                                                                                   >  1:532393/1‑277 (MQ=255)
                                                                                                                                        tAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTACGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaa                                                                                                                                                                >  1:461129/1‑301 (MQ=255)
                                                                                                                                                    tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATaa                                                                                                                                                    >  1:28111/1‑301 (MQ=255)
                                                                                                                                                     atTCAGGCAGATAGAGATAAAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTAGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTAATGACCAACAGGGAGACAAATGACCACATAAc                                                                                                                                                   >  1:154182/1‑301 (MQ=255)
                                                                                                                                                            ttgATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTc                                                                                                                                            <  2:38393/299‑1 (MQ=255)
                                                                                                                                                              tccgatcttataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAGGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGc                                                                                                                                          <  2:7511/293‑1 (MQ=255)
                                                                                                                                                              gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTcag                                                                                                                                          >  1:38341/1‑299 (MQ=255)
                                                                                                                                                                      ttctaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGAcc                                                                                                                                  <  2:28784/297‑1 (MQ=255)
                                                                                                                                                                      tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAGGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATAAGGGTGACAAATGACCACAAAACCCCTTTCGcagatcgga                                                                                                                                  >  1:7511/1‑293 (MQ=255)
                                                                                                                                                                          aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACAAAACCTCTTTCGCTCCGGAccagaa                                                                                                                              >  1:28733/1‑297 (MQ=255)
                                                                                                                                                                                        aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGc                                                                                                                <  2:461693/301‑1 (MQ=255)
                                                                                                                                                                                             cAGGCGAGTAAATTCAAATAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTag                                                                                                           <  2:154333/301‑1 (MQ=255)
                                                                                                                                                                                                        aTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCaaa                                                                                                >  2:316642/1‑300 (MQ=255)
                                                                                                                                                                                                         ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGAAGAAGAGGAAGGCAATg                                                                                               >  2:149820/1‑301 (MQ=255)
                                                                                                                                                                                                          tCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCACTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGAAGAGGAAGGCAAAgg                                                                                              >  1:39847/1‑301 (MQ=255)
                                                                                                                                                                                                             aaGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATCCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTCCTTAAGGTTGGCTTTAACGTTTGCTGGATTAAGGGGAATTCATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCACATACATGGCGGTGATGACCATGAGGGTGACAAATGACCACATAACCTCATTCGCTCCGGACCCCAACACGGCCCACAAGCAGTAGAGGAAGGCAATGgaag                                                                                           >  1:18250/1‑301 (MQ=255)
                                                                                                                                                                                                                               ggCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACAACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAAATGCAAGATATGACg                                                                         >  2:285449/1‑301 (MQ=255)
                                                                                                                                                                                                                                   cgtCAACAGCTCGGTTTTTCCATCGAACAGAGTGCGAATATGATCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTTGCTTTTATGTGTGCTGGATTAAGGTGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGAGAATCAGGGTGACAAATGACCACATAACCTCTTACGCTCCGGACCACACCACGGCCCGGATGCAGTAGAGGAAGGCAAGGGGAGTAACTGCCAGATATggcgggcg                                                                     <  1:95258/299‑1 (MQ=255)
                                                                                                                                                                                                                                        aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCtt                                                                >  2:611220/1‑301 (MQ=255)
                                                                                                                                                                                                                                                    ttCCATCGAACAAAGTGCGCATATGCTCTCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGAc                                                    <  2:141360/301‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTCTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACAGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGAc                                                    >  2:349357/1‑301 (MQ=255)
                                                                                                                                                                                                                                                     tCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGAcc                                                   >  2:532838/1‑301 (MQ=255)
                                                                                                                                                                                                                                                       cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTCGAGGAAGGCAATGGTAGAAACTGCCAGAAATGCCGGGCGTGATTTACCAAAGTAACCgg                                                 >  1:93486/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                   gctcttccgattttgCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAAGCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCACAg                       <  1:533601/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         a‑aT‑TAAGTGTACCTTAATTTCTTAAGGTTGGCTTTTATGGTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTGTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCACAGGTGTAAAg               <  2:18301/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         a‑aT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCACGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCACAGGTGTaaaa               >  2:161907/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                aGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCACAgagatcggaagaga          >  2:534274/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                                 gCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGACAAGGAGTAAAGAAATgg         >  1:630166/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                  ccAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCACAGGAGTAAAGATATgga        >  2:93840/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                   cAACCTTAAGAACTTAAGGTTGGC‑TTAAT‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCCGGGTGACAACTGACCACCTAACCTCTTTCGCTCCGGACCCCACCCCGGCCCAGATGCCGTAGAGGAAGGCAATGGTAGTAACTGCCAGACATACCGGGCGTGCTTTACAAAAGTAACCACGTCCGAGCAACAATAACGCAGCACAACTGTAAAAAAATgacccc     >  1:512077/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                                                       tttAAGTTCTTAAGGTTGGTTTTTATGTTTGCTGGATTTAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTGTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTGTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCACAGGTGTAAAGATATGGCACCAg   <  1:532165/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                        ttAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCACAGGTGTAAAGATATGGCACCAgt  <  2:93610/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                       
GCCAGTGGTCGATTAAGCCATAAGAAAGAGCTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCACAGGTGTAAAGATATGGCACCAGT  >  NC_000913/4335355‑4335949

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑