Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,463,952 G→A S253F (TCC→TTC)  treC ← trehalose‑6‑phosphate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,463,9520GA94.9% 166.1 / ‑1.8 59S253F (TCC→TTC) treCtrehalose‑6‑phosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (25/31);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

ATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTACCGCTGGATTCTCCCA  >  NC_000913/4463705‑4464193
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
atcaCCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACa                                                                                                                                                                                                                                                <  2:82458/251‑1 (MQ=255)
  ccccAAAGCGAGAAACAATGCGCGGCTGATCTTGGTTACTCCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCAGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACAtt                                                                                                                                                                                                                                              <  1:296513/249‑1 (MQ=255)
  caCCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCCGGTGGAGAACAtt                                                                                                                                                                                                                                              >  1:54831/1‑251 (MQ=255)
     cAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGGAAATTAAAGGTCATCTACAATTCACTGCCTGTCAGAGCCCCGTATCGCTGGCAAAGCTAAAGGCTGGTGGAGGACATTTcc                                                                                                                                                                                                                                           >  1:124849/1‑250 (MQ=255)
         gCGAGAAACAATGCGCGGCTGATCGTGGTTACCACAGAACAAGGCATTCCATGCTAAGTCGTGAATTCCTTGTTGCCAGGGGCAGAAAAATGTTTCCAACGCCACAAAGTCAGGTTTAGCTAGCGTCCATTTTTCCCCTCCGGGATAAGCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCTTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTa                                                                                                                                                                                                                                       <  2:181247/251‑1 (MQ=255)
                        cgGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGGTTTCCACGCCACAAAGTCCGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCCGGTGATGAAAATTAAAGGCCATCGCCAATTACCTGCCTGTAAAAGCCGCGTATCCCTGGCAATGCTCAAGGTTGGGGGGGAACATTTACCCTCCGGTCATTACCCcgc                                                                                                                                                                                                                        >  2:176399/1‑251 (MQ=255)
                          gCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGGCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGt                                                                                                                                                                                                                      <  1:437328/251‑1 (MQ=255)
                          gCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAAAGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTAACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCAATGCCTGTCAGAGCCGCGTATCGCTTGCAATGCTCAAGGCTGGTGGAGAAAATTTCAACTACGGGAATTAACCCACGt                                                                                                                                                                                                                      >  1:323225/1‑251 (MQ=255)
                                   ggTTACACCAGAACAAGGCTTTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAac                                                                                                                                                                                                             <  2:382228/251‑1 (MQ=255)
                                           cAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATcgcg                                                                                                                                                                                                     <  2:443376/251‑1 (MQ=255)
                                                     cATTCCATGCTACGTTGTGCATTTCGTGTTGCCAGTGGCGTAACAATGTTGTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTGCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTAGCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCTCGGTTCATCTCg                                                                                                                                                                                           <  2:216330/251‑1 (MQ=255)
                                                        tCCATGCTACGTTGTGCATTCCTGGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTTGCGTAAACACATCGCGGTTCATCTCGTg                                                                                                                                                                                         <  2:378785/250‑1 (MQ=255)
                                                         ccATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCa                                                                                                                                                                                       <  2:206554/251‑1 (MQ=255)
                                                           aTGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCaaa                                                                                                                                                                                     <  1:145145/251‑1 (MQ=255)
                                                               tACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACACCACAAAATCAGGGTTAGCCAGCCTCCATTTTTCACCACCCGGATAATCGACCTTCAGGTGATGAAAATTAAAGGGCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGGGGAGAACATTTCACCTTCGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACt                                                                                                                                                                                 >  1:149289/1‑251 (MQ=255)
                                                                 cGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCg                                                                                                                                                                               <  1:264761/251‑1 (MQ=255)
                                                                       gCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCt                                                                                                                                                                         <  1:483732/251‑1 (MQ=255)
                                                                                   gCCAGGGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAAATAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTc                                                                                                                                                              <  1:176399/250‑1 (MQ=255)
                                                                                       gTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGgtgt                                                                                                                                                         <  1:29841/251‑1 (MQ=255)
                                                                                                 aaTGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACg                                                                                                                                               <  1:132652/251‑1 (MQ=255)
                                                                                                  aTGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGc                                                                                                                                              >  1:135314/1‑251 (MQ=255)
                                                                                                       ttCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTACAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTc                                                                                                                                         >  1:484468/1‑251 (MQ=255)
                                                                                                              ccACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTAGAGAAAATTTCACCTACGGTCATTACCCCGCGTGGCGTAAACACATCGCGGATAATCTCGTACAAAACCTCGCGTGCTCGTGGCCCGTCGGTGTAGAAGCCACGCCCGTCGCCGTcc                                                                                                                                  >  1:437834/1‑251 (MQ=255)
                                                                                                               cacaAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAACCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCa                                                                                                                                 <  1:186163/251‑1 (MQ=255)
                                                                                                                  aaaGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTaa                                                                                                                                                                                                               <  2:409788/170‑1 (MQ=255)
                                                                                                                  aaaGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTaa                                                                                                                                                                                                               >  1:409788/1‑170 (MQ=255)
                                                                                                                             tAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTGTAGGTGATTAAAATTAAAGGTCATCTACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCGCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGaa                                                                                                                    <  2:448284/250‑1 (MQ=255)
                                                                                                                                ccAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAAcgc                                                                                                                >  1:219547/1‑251 (MQ=255)
                                                                                                                                  aGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGCCCAGGTCTTCAGGGAAACGCgg                                                                                                              >  1:354350/1‑251 (MQ=255)
                                                                                                                                                            aTCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACAt                                                                                    >  1:161633/1‑251 (MQ=255)
                                                                                                                                                                         tgatgaAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACg                                                                                                                                               <  2:457597/179‑1 (MQ=255)
                                                                                                                                                                         tgatgaAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACg                                                                                                                                               >  1:457597/1‑179 (MQ=255)
                                                                                                                                                                              aaaaTTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCg                                                                   >  1:394363/1‑250 (MQ=255)
                                                                                                                                                                                  ttAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGAcccc                                                              >  2:114260/1‑251 (MQ=255)
                                                                                                                                                                                       ggTCATCGACAATTAACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGACACCCGTCGAccccccgg                                                          >  2:421992/1‑246 (MQ=255)
                                                                                                                                                                                         tCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCg                                                       >  1:283653/1‑251 (MQ=255)
                                                                                                                                                                                              ggtaaTTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCAGCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGACTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCa                                                  <  2:480346/248‑1 (MQ=255)
                                                                                                                                                                                              gACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGcccc                                                  >  1:273780/1‑250 (MQ=255)
                                                                                                                                                                                                 aaTTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGa                                                >  1:350710/1‑250 (MQ=255)
                                                                                                                                                                                                 aaTTCACTGCCCGGCAGAGCCGCGTTTCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACAGCCAGGCGCAACCCGTCGACCCCACGATCGGTCCAGaa                                               <  2:54831/251‑1 (MQ=255)
                                                                                                                                                                                                     cACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTca                                           <  1:215284/251‑1 (MQ=255)
                                                                                                                                                                                                          ccTGTCAGAGCCGTGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGCCGACCCCACGATCGGCCCAGAACTCACAGAc                                      <  1:106241/251‑1 (MQ=255)
                                                                                                                                                                                                           cTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACt                                     >  1:236384/1‑251 (MQ=255)
                                                                                                                                                                                                             gTCATAGCCGCGTATCGCTGGCCATGCTCAAGGGTGGTGTAGAACATGTCACCTACGGTCATTAGCCCGCGTGGCGTATACCCATCGCGGGTCATCTCGTGCAAAAACGCGTGTGCTCGTGGCCCGGCGGTGTAGAAGCGACGCCCGTCGGCGGCCAGGTTTTCAGGGAAGCGCGGGTCTTTGGAGATCAGATTCACCACATCCGTGCGCAACCCGTCGTCCCCACGATCGGCCCAGAACTCACAGACt                                     <  2:288137/249‑1 (MQ=255)
                                                                                                                                                                                                              tCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCATGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACACACAGACtttt                                  >  1:197670/1‑251 (MQ=255)
                                                                                                                                                                                                               cAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACtttt                                  <  1:114260/250‑1 (MQ=255)
                                                                                                                                                                                                                agagCCGCGTATCGCTGGCAATGCTCAAGGCTGGAGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGCAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCCCGATCGGCCCGGAACTCACAGACtttttt                                >  1:316378/1‑251 (MQ=255)
                                                                                                                                                                                                                 gagCCGCGTATAGCTGGCAATGCTCAAGGCTGGTGGCGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCCCATCCCGGCGCAACCCGTCGACCCCCCGCTCGGCCCCGAACTCACAGACTTTTTTc                               >  2:291647/1‑251 (MQ=255)
                                                                                                                                                                                                                  agCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCa                              <  1:141566/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGCCTTTTTTCAGc                            >  2:229313/1‑251 (MQ=255)
                                                                                                                                                                                                                        gaatcccTGGCAATTCTCAAGGCTGGTAGAGAACATGTCACCTCCGGTCATTAGCCCGCGTGTCGTAAACACACCGCTGTTCATCTCGTGCAAAACCTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCTGATTCACCACATCCAGGCGCATCCCGTCGACCCCACTATCTGCCCAGATCTTACAGACTTTTTGCAGctcgg                        <  2:324917/245‑3 (MQ=255)
                                                                                                                                                                                                                         tATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGCTCAGATTCACCACATCCAGGCGAAACCCGTCGACCCCACGATCGGCCCAGAACTCCCAGACTTTTTTCAGCTCTgc                       >  2:449221/1‑251 (MQ=255)
                                                                                                                                                                                                                          aTCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTgcg                      <  2:149289/251‑1 (MQ=255)
                                                                                                                                                                                                                              cTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGt                   >  1:48008/1‑250 (MQ=255)
                                                                                                                                                                                                                                   aaTGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCCCGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTACCGCt             >  1:405945/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aaTGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCATCTCTGCGCGTTCCGCt             <  2:165039/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTACCGCTGGAtt        <  1:449221/250‑1 (MQ=255)
                                                                                                                                                                                                                                                ggtggAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTACCGCTGGATTCTCCCa  <  1:456369/249‑1 (MQ=255)
                                                                                                                                                                                                                                                ggtggAGAACATTTCACCTACGGGCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTACCGCTGGATTCTCCCa  <  2:280830/249‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
ATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTACCGCTGGATTCTCCCA  >  NC_000913/4463705‑4464193

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