Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,926,991 G→A L5L (CTG→TTG)  asnC ← DNA‑binding transcriptional dual regulator AsnC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,926,9910GA94.0% 135.5 / ‑0.3 51L5L (CTG→TTG) asnCDNA‑binding transcriptional dual regulator AsnC
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (23/24);  total (24/27)
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.60e-01

GCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTC  >  NC_000913/3926750‑3927233
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
gcgGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGCTTGTCGATCAAATAATTTTc                                                                                                                                                                                                                                             >  1:252048/1‑251 (MQ=255)
  ggAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTcc                                                                                                                                                                                                                                            >  1:159335/1‑250 (MQ=255)
   gAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCAt                                                                                                                                                                                                                                          >  1:154337/1‑251 (MQ=255)
                          aatgaaacacccGATAAATCAGCCCACGCTATACCCGAGCTGCTTCGGTTTGACATCAATACGCGCCCCGGTAATGTCTCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTTTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTGTCTTATGCGGATtgat                                                                                                                                                                                                                   <  2:480125/241‑1 (MQ=255)
                           aaTATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATtgat                                                                                                                                                                                                                   >  1:480125/1‑250 (MQ=255)
                           aaTAAAATGCCGATAAAGCAGCCTACGTCAAAACCGGGCTGCTTCGGGCTGACATCAATCCGCGCCCCGGTAATGACCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATtgatg                                                                                                                                                                                                                  <  2:159335/251‑1 (MQ=255)
                                    ccGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCt                                                                                                                                                                                                         <  1:135801/251‑1 (MQ=255)
                                        aaaatcagcCTACGTCATAACCGAGCTACTTCGGGCTGCCCTCAACTCGCTCCCCTGTAATGACCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCTCGCATTTCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATTATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTAttt                                                                                                                                                                                                     <  2:457943/246‑1 (MQ=255)
                                               gtcctctcccgatctgggctgctTCGGGCTGACATCAATACGCGCCCCGGTAACTATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCctt                                                                                                                                                                                              <  1:458742/234‑1 (MQ=255)
                                                 ctaagTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTACTGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGAAGCCAAATTGATTCGCCAATTCGGCGTAACCGGTGCGCGCATTACCAATTAATGCTTCCAGGATGCCACGGTACAAATTGTCGATCAAATAATTTTACTTAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAACcttat                                                                                                                                                                                            >  1:412808/5‑249 (MQ=255)
                                                 cTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCcttct                                                                                                                                                                                            >  1:412814/1‑251 (MQ=255)
                                                       aaaaaCCGAGTTGCGTCGGGTTTACATAAATCCGCGCCCCGGTAATGACCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTa                                                                                                                                                                                      <  2:294166/248‑1 (MQ=255)
                                                            ccGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGa                                                                                                                                                                                  >  1:485083/1‑250 (MQ=255)
                                                              gAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAACGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGGGCGCGCATTGCCCATTAATGCTTCCAGGATGCACCGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGGTGATTCATTATATTTTAGCctt                                                                                                                                                                                              >  2:458742/1‑236 (MQ=255)
                                                                      tCGGGCTTACATAAACACGGGCCCCTGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTg                                                                                                                                                                        <  2:568121/250‑1 (MQ=255)
                                                                       cGGGATGACAGCAAGACGCGCCCCTGATATGACCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCATATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAg                                                                                                                                                                     <  2:603709/251‑1 (MQ=255)
                                                                        ggggatacaccAATACGCCCCCCGGTACTGCTCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAAACGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTTCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGTATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGt                                                                                                                                                                    <  2:336155/241‑1 (MQ=255)
                                                                                aTCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGACCGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCAGCAGTTCGGCATAAGCGGTGCGCGCATTGCCCATTCATGCTTCCAGGGTGCCACGGTCCAGATTATCGATCAAATAATTTTCCATAGGATTTTCTTATGCGACTTGATGGTTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCttt                                                                                                                                                            >  1:222042/1‑251 (MQ=255)
                                                                                    aTACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGCTTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTAt                                                                                                                                                        >  1:413575/1‑251 (MQ=255)
                                                                                               ggTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCCGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCATAAGTGAGTTAGGCTTTTTATTGAATGattat                                                                                                                                             >  2:515143/1‑251 (MQ=255)
                                                                                                   aTGATCACCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTGATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGCTTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAATAGTGAGTTAGGCTTTTTATTGAATGATTATTGc                                                                                                                                          >  2:21454/1‑250 (MQ=255)
                                                                                                          cccGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCAtgtgtg                                                                                                                                   <  1:489634/251‑1 (MQ=255)
                                                                                                           ccGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCCAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCa                                                                                                                                                                              >  1:284467/1‑208 (MQ=255)
                                                                                                             cccTGCTTCATTTTCTCTTCTCTAACGTTAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCTCGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTTTGTCg                                                                                                                                <  2:367106/250‑1 (MQ=255)
                                                                                                                 gCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAGAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGttt                                                                                                                            <  2:485083/251‑1 (MQ=255)
                                                                                                                 gCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGttt                                                                                                                            <  1:111656/251‑1 (MQ=255)
                                                                                                                 gCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGGTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGttt                                                                                                                            <  2:265075/251‑1 (MQ=255)
                                                                                                                 gCTTCATTTTCTCTACTCGAACGTGAATCGGCCCCGGACTGACGTCAAAGTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGtt                                                                                                                             <  1:141447/250‑1 (MQ=255)
                                                                                                                    tCATTTTCTATACTCGAAAGAGAATCGTCCCCGGCATGAAGCCAAGTTGTTTCGCAAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACAGACCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGCTGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTttg                                                                                                                         >  2:449399/1‑251 (MQ=255)
                                                                                                                         ggggacgattcaccGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCt                                                                                                                    >  1:87000/14‑251 (MQ=255)
                                                                                                                               ccacgaacggtaaTCGGCCCCGGACTGACTCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTTACCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTGTGTGTTGCTTAAta                                                                                                               <  2:640050/239‑2 (MQ=255)
                                                                                                                               aCTCGAACGTGAATCGTCCCCTGACTGACGCCATATTGTTTCGCCAGTTCGTCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGTCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCa                                                                                                              <  1:657879/251‑1 (MQ=255)
                                                                                                                                         gAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAg                                                                                                   <  1:133887/251‑1 (MQ=255)
                                                                                                                                          aaTCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCCTGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAgg                                                                                                  <  1:60511/251‑1 (MQ=255)
                                                                                                                                                ccccGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAg                                                                                            <  2:190000/251‑1 (MQ=255)
                                                                                                                                                     gACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGta                                                                                       <  1:532736/251‑1 (MQ=255)
                                                                                                                                                         gACGCCAAATTGTTTCTCCAGTGCGGCGTAAGCGGATCGTGCATTGCCCATTAATGCTTCCTGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATaa                                                                                    <  2:469847/250‑1 (MQ=255)
                                                                                                                                                              cAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATg                                                                              >  1:169237/1‑251 (MQ=255)
                                                                                                                                                               aaaTTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATg                                                                              >  1:310890/1‑250 (MQ=255)
                                                                                                                                                                          tcAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTATATTTTAGCCTTCTTTTTTAATGAATCAAAGTGGAGTTAGGCTTTTTATTGAATGATTATTGTATGTGTGTCGGTTTTTGTTGCTTAATCAGAAGAACACGGACGCAGGAGTATAAAACATAAAAAcctctaac                                                                  >  2:116641/2‑245 (MQ=255)
                                                                                                                                                                                  gCGTAAGCGGTTCTCGCATTGCCCATTTATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAATCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCaa                                                          <  2:471108/251‑1 (MQ=255)
                                                                                                                                                                                    gTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAaac                                                        >  1:177761/1‑251 (MQ=255)
                                                                                                                                                                                      aaGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGAGTGATGATTCATTCTATTTTAGCCTTCTTTGTTAATGAATCAAAAGTGATTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATACAAAATGAAAACCGCTCACATTGCCAaacaa                                                      <  2:63770/251‑1 (MQ=255)
                                                                                                                                                                                                           ttaatgtgtCCATGATGCCTCGGTCCAGATTGTCGATCTAATAATTTTCCATATGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTTTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTTCCAATCTACGTTAAATTAGCTTCGGGaa                                  <  2:392769/242‑1 (MQ=255)
                                                                                                                                                                                                                 cTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCAc                           >  1:608992/1‑251 (MQ=255)
                                                                                                                                                                                                                 cTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGAATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCAc                           >  1:391607/1‑251 (MQ=255)
                                                                                                                                                                                                                   tacAGGATGCCACGGTCCAGATTGTAGATCAAATAATTTTCCATAGGATTTTCTTATGCGGAGTGATGATTCATTCTAATTTAGCCTTCTGTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTAGTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGAAGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACt                          <  2:559993/248‑1 (MQ=255)
                                                                                                                                                                                                                    ccAGGATTCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTTCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACttt                        <  2:572019/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTAGTCAACTgaacgac      >  2:204080/1‑245 (MQ=255)
                                                                                                                                                                                                                                       aTTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGCGTATAAAAAATGAAAAACGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAAc     >  1:157528/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTc  >  2:200275/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
GCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTC  >  NC_000913/3926750‑3927233

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