Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,312,289 A→G Y235H (TAT→CAT)  alsR ← DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,2890AG98.0% 181.6 / NA 50Y235H (TAT→CAT) alsRDNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (28/21);  minor base T (1/0);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATT  >  NC_000913/4312041‑4312535
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
gtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATc                                                                                                                                                                                                                                                     <  1:180834‑M2/235‑1 (MQ=255)
gtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATc                                                                                                                                                                                                                                                     <  1:301308‑M2/235‑1 (MQ=255)
  cagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTCTTTCAGCGCGCCTGGAGAAAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCgg                                                                                                                                                                                                                                                   >  1:548359‑M2/15‑251 (MQ=255)
    gttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCACGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTCAGCAAATAATATGCTCGGcc                                                                                                                                                                                                                                                 >  2:132120‑M2/13‑251 (MQ=255)
         gacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGCAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACAGCGTTTCCGGGGCTGGTGCGCAAATAATATGATCGGCCAGttt                                                                                                                                                                                                                                            >  1:366982‑M2/8‑251 (MQ=255)
                         gacagTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCTCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATg                                                                                                                                                                                                                             <  1:8363/247‑1 (MQ=255)
                            agTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGCGAAGAAAACAACAATTGCGCCGGCTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGAGAGCAAATAATATGATCGGCCAGTTTCGCTACCGGTGAATGGtt                                                                                                                                                                                                                          >  2:561815/2‑250 (MQ=255)
                               aaCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGc                                                                                                                                                                                                                        <  1:78319/250‑1 (MQ=255)
                                     gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGt                                                                                                                                                                                                                 <  1:159248/251‑1 (MQ=255)
                                      cgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGCGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCACGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCCGTTTCGCTATCGGTGAATGGTAGCTATGGGtt                                                                                                                                                                                                                >  2:82071/1‑251 (MQ=255)
                                       gATGGCGGAATTACCTTTTTCTTATTTCAGTGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTta                                                                                                                                                                                                               <  2:322053/251‑1 (MQ=255)
                                             ggAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACaa                                                                                                                                                                                                          >  1:263229/1‑250 (MQ=255)
                                                       ttttCTTATTTCAGCGCGCCTTGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGTGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTg                                                                                                                                                                                                <  1:323889/250‑1 (MQ=255)
                                                              aTTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTc                                                                                                                                                                                                                                           <  1:388291/200‑1 (MQ=255)
                                                              aTTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGtt                                                                                                                                                                                         <  1:602985/250‑1 (MQ=255)
                                                              aTTTCAGCGCGCCTGGCGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGCCAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTc                                                                                                                                                                                                                                           >  2:388291/1‑200 (MQ=255)
                                                                ttCAGCGCGCCTGGTGAGAAGAAATCAACAATTGAGCCGGTTTTTTCCATATTCATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGCAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCtt                                                                                                                                                                                      >  2:33643/1‑251 (MQ=255)
                                                                  cAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGTTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCtttt                                                                                                                                                                                    <  2:552765/251‑1 (MQ=255)
                                                                  cAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGGTCAAGGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATAGTCTTGCAGAGGCATTACGACCTAATATCGGCGTTTCCGGGGCTGGTGAGCAAATAATAGGATCGGCCAGTTGCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCtttt                                                                                                                                                                                    <  2:199760/251‑1 (MQ=255)
                                                                    gcgcgcCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCtttttt                                                                                                                                                                                  <  1:453317/251‑1 (MQ=255)
                                                                               gagaAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTa                                                                                                                                                                        <  1:378777/250‑1 (MQ=255)
                                                                               gagaAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTCGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTa                                                                                                                                                                        <  2:597834/250‑1 (MQ=255)
                                                                                         caacaaTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCtttt                                                                                                                                                             >  1:587208/1‑251 (MQ=255)
                                                                                              aTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGCCAGAGACAAAAAAAGCGACCAGCAAGGCTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATc                                                                                                                                                        >  1:555407/1‑251 (MQ=255)
                                                                                                      ggTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTc                                                                                                                                                <  1:71423/251‑1 (MQ=255)
                                                                                                             tccATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTAAGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCggg                                                                                                                                         >  2:529501/3‑250 (MQ=255)
                                                                                                                                   cAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCacaac                                                                                                                    >  2:34038/1‑250 (MQ=255)
                                                                                                                                     aTGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTTAGCAAATAATATGATCGGCCAGTTTCGTTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTGTTCCCCGTTCTTTTTTGCCAGTTCTCCGTCCGCTTTTGCATCACTGTTTCGCCCTGAATGGGTCGCTATCTGCACAACAtc                                                                                                                 <  2:144389/251‑1 (MQ=255)
                                                                                                                                       gTTGAGCTGGGCAACAGAGACAAACACAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATTATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCTATCTTTTTTGCCAGTTCTACGGCCTCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCCCAACAtctc                                                                                                               >  2:276748/1‑251 (MQ=255)
                                                                                                                                                ggcgtaagagACAAAAAAAGCGTCCAGAAAGGTTGGTTGTAATGTTCTTGCCGAGGAATTACGACGTAATAACGGCGTTTACGGGGCTGGTGAGCAAATAATATGATCGGCCCGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCAGTTCTTTTTTGCCAGTTCTACGGACGCTTTTACATCTCTGGTTAGCCCGGAATGGGTCACTAACAGCACAAAATATCCTTCATGt                                                                                                       >  2:510699/7‑250 (MQ=255)
                                                                                                                                                       gagaCAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACg                                                                                                >  2:187283/1‑250 (MQ=255)
                                                                                                                                                             aaaaaaGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGCAGCgg                                                                                          >  2:517384/1‑250 (MQ=255)
                                                                                                                                                                     gTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCCTTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTCCGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCCCTACCAGCACAACATCTCCTTACTGTAAAAACGAAGCGGAcatcatc                                                                                  >  2:335558/1‑250 (MQ=255)
                                                                                                                                                                           caaggttggttGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGAGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGAtgt                                                                           >  1:510823/9‑251 (MQ=255)
                                                                                                                                                                             aggttggttGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGAtgtg                                                                          >  2:615548/7‑250 (MQ=255)
                                                                                                                                                                                                 gCCGAGGCATTTCGATCTAATAACGTCGTTTCCGGGGCTGGTGAGCAAATAATATGATCTGTCATTTTCGCTATCTGTGAATGGTAGCTTTGTGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTCCGGCCGCTTTTTCATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGAc                                                      <  2:631880/250‑1 (MQ=255)
                                                                                                                                                                                                 gCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTaac                                                                                                    >  2:109798/1‑204 (MQ=255)
                                                                                                                                                                                                 gCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTaac                                                                                                    <  1:109798/204‑1 (MQ=255)
                                                                                                                                                                                                 gCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATATTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCCTGATGTAAGCATCAGGATAGGCCTGAcc                                                     >  2:107081/1‑250 (MQ=255)
                                                                                                                                                                                                         aTTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGc                                             >  2:225766/1‑251 (MQ=255)
                                                                                                                                                                                                            aCGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATAGGATCGGCCAGTTTCGCTATCGGTGCATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCCCAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCCTCAGGATAGGCCTGACCGCGTACGCCa                                           >  2:558784/1‑250 (MQ=255)
                                                                                                                                                                                                                  taataaCGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATgggc                                     >  1:69442/1‑247 (MQ=255)
                                                                                                                                                                                                                          gCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTtg                            <  1:531811/251‑1 (MQ=255)
                                                                                                                                                                                                                           cGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTtgt                           >  1:200049/1‑251 (MQ=255)
                                                                                                                                                                                                                             tttCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGc                         >  1:250278/1‑251 (MQ=255)
                                                                                                                                                                                                                                ccGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTg                       <  1:202600/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gCTGGGGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATTGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATc                 <  1:43359/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGGTGAGCAAATAATATGATCGGCCCGTTTCGCTATCGGTGAATGGTAGCTATGAGTTATACAAATCATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGAATCAGGATAGGCCTGACAGCGTACGCCAATTGGCAAGAACTTGTGCTGTACATc                 >  2:406706/1‑250 (MQ=255)
                                                                                                                                                                                                                                         gTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGTATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGc              >  1:210807/1‑250 (MQ=255)
                                                                                                                                                                                                                                             gCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGcac            >  1:167336/1‑248 (MQ=255)
                                                                                                                                                                                                                                                   aatatgatCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCAtt  >  2:177508/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    atatgatCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCAGAGATAGCAtt  <  1:82071/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
ATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATT  >  NC_000913/4312041‑4312535

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