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

GCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTC  >  NC_000913/3926750‑3927233
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
gcgGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGCTTGTCGATCAAATAATTTTc                                                                                                                                                                                                                                             >  1:252048/1‑251 (MQ=255)
  ggAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTcc                                                                                                                                                                                                                                            >  1:159335/1‑250 (MQ=255)
   gAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCAt                                                                                                                                                                                                                                          >  1:154337/1‑251 (MQ=255)
                          aatgaaacacccGATAAATCAGCCCACGCTATACCCGAGCTGCTTCGGTTTGACATCAATACGCGCCCCGGTAATG‑TCTCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTTTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTGTCTTATGCGGATtgat                                                                                                                                                                                                                   <  2:480125/241‑1 (MQ=255)
                           aaTATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATtgat                                                                                                                                                                                                                   >  1:480125/1‑250 (MQ=255)
                           aaTAAAATGCCGATAAAGCAGCCTACGTCAAAACCGGGCTGCTTCGGGCTGACATCAATCCGCGCCCCGGTAATGACC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATtgatg                                                                                                                                                                                                                  <  2:159335/251‑1 (MQ=255)
                                    ccGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCt                                                                                                                                                                                                         <  1:135801/251‑1 (MQ=255)
                                        aaaatcagcCTACGTCATAACCGAGCTACTTCGGGCTGCCCTCAACTCGCTCCCCTGTAATGACC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCTCGCATTTCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATTATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTAttt                                                                                                                                                                                                     <  2:457943/246‑1 (MQ=255)
                                               gtcctctcccgatctgggctgctTCGGGCTGACATCAATACGCGCCCCGGTAACTATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCctt                                                                                                                                                                                              <  1:458742/234‑1 (MQ=255)
                                                 ctaagTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTACTGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGAAGCCAAATTGATTCGCCAATTCGGCGTAACCGGTGCGCGCATTACCAATTAATGCTTCCAGGATGCCACGGTACAAATTGTCGATCAAATAATTTTACTTAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAACcttat                                                                                                                                                                                            >  1:412808/5‑249 (MQ=255)
                                                 cTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCcttct                                                                                                                                                                                            >  1:412814/1‑251 (MQ=255)
                                                       aaaaaCCGAGTTGCGTCGGGTTTACATAAATCCGCGCCCCGGTAATGACC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTa                                                                                                                                                                                      <  2:294166/248‑1 (MQ=255)
                                                            ccGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGa                                                                                                                                                                                  >  1:485083/1‑250 (MQ=255)
                                                              gAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAACGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGGGCGCGCATTGCCCATTAATGCTTCCAGGATGCACCGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGGTGATTCATTATATTTTAGCctt                                                                                                                                                                                              >  2:458742/1‑236 (MQ=255)
                                                                      tCGGGCTTACATAAACACGGGCCCCTGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTg                                                                                                                                                                        <  2:568121/250‑1 (MQ=255)
                                                                       cGGGATGACAGCAAGACGCGCCCCTGATATGACC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCATATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑Ag                                                                                                                                                                     <  2:603709/251‑1 (MQ=255)
                                                                        ggggatacaccAATACGCCCCCCGGTACTGCTC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAAACGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTTCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGTATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGt                                                                                                                                                                    <  2:336155/241‑1 (MQ=255)
                                                                                aTCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGACCGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCAGCAGTTCGGCATAAGCGGTGCGCGCATTGCCCATTCATGCTTCCAGGGTGCCACGGTCCAGATTATCGATCAAATAATTTTCCATAGGATTTTCTTATGCGACTTGATGGTTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCttt                                                                                                                                                            >  1:222042/1‑251 (MQ=255)
                                                                                    aTACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGCTTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTAt                                                                                                                                                        >  1:413575/1‑251 (MQ=255)
                                                                                               ggTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCCGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCATAAGTG‑AGTTAGGCTTTTTATTGAATGattat                                                                                                                                             >  2:515143/1‑251 (MQ=255)
                                                                                                   aTGATC‑ACCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTGATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGCTTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAATAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGc                                                                                                                                          >  2:21454/1‑250 (MQ=255)
                                                                                                          cccGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCAtgtgtg                                                                                                                                   <  1:489634/251‑1 (MQ=255)
                                                                                                           ccGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCCAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCa                                                                                                                                                                              >  1:284467/1‑208 (MQ=255)
                                                                                                             cccTGCTTCATTTTCTCTTCTCTAACGTTAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCTCGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTTTGTCg                                                                                                                                <  2:367106/250‑1 (MQ=255)
                                                                                                                 gCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAGAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGttt                                                                                                                            <  2:485083/251‑1 (MQ=255)
                                                                                                                 gCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGttt                                                                                                                            <  1:111656/251‑1 (MQ=255)
                                                                                                                 gCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGGTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGttt                                                                                                                            <  2:265075/251‑1 (MQ=255)
                                                                                                                 gCTTCATTTTCTCTACTCGAACGTGAATCGGCCCCGGACTGACGTCAAAGTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGtt                                                                                                                             <  1:141447/250‑1 (MQ=255)
                                                                                                                    tCATTTTCTATACTCGAAAGAGAATCGTCCCCGGCATGAAGCCAAGTTGTTTCGCAAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACAGACCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGCTGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTttg                                                                                                                         >  2:449399/1‑251 (MQ=255)
                                                                                                                         ggggacgattcaccGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCt                                                                                                                    >  1:87000/14‑251 (MQ=255)
                                                                                                                               ccacgaacggtaaTCGGCCCCGGACTGACTCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTTACCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTGTGTGTTGCTTAAta                                                                                                               <  2:640050/239‑2 (MQ=255)
                                                                                                                               aCTCGAACGTGAATCGTCCCCTGACTGACGCCATATTGTTTCGCCAGTTCGTCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGTCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCa                                                                                                              <  1:657879/251‑1 (MQ=255)
                                                                                                                                         gAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑g                                                                                                   <  1:133887/251‑1 (MQ=255)
                                                                                                                                          aaTCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCCTGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑gg                                                                                                  <  1:60511/251‑1 (MQ=255)
                                                                                                                                                ccccGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAg                                                                                            <  2:190000/251‑1 (MQ=255)
                                                                                                                                                     gACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGta                                                                                       <  1:532736/251‑1 (MQ=255)
                                                                                                                                                         gACGCCAAATTGTTTCTCCAGTGCGGCGTAAGCGGATCGTGCATTGCCCATTAATGCTTCCTGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATaa                                                                                    <  2:469847/250‑1 (MQ=255)
                                                                                                                                                              cAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATg                                                                              >  1:169237/1‑251 (MQ=255)
                                                                                                                                                               aaaTTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATg                                                                              >  1:310890/1‑250 (MQ=255)
                                                                                                                                                                          tcAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTATATTTTAGCCTTCTTTTTTAATGAATCAAA‑GTGGAGTTAGGCTTTTTATTGAATGATTATTGTATGTGTGTCGGTTTTTGTTGCTTAATCAGAAG‑AACACGGACGCAGGAGTATAAAACATAAAAAcctctaac                                                                  >  2:116641/2‑245 (MQ=255)
                                                                                                                                                                                  gCGTAAGCGGTTCTCGCATTGCCCATTTATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAATCAACA‑GGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCaa                                                          <  2:471108/251‑1 (MQ=255)
                                                                                                                                                                                    gTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAaac                                                        >  1:177761/1‑251 (MQ=255)
                                                                                                                                                                                      aaGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGAGTGATGATTCATTCTATTTTAGCCTTCTTTGTTAATGAATCAAAAGTG‑ATTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATACAAAATGAAAACCGCTCACATTGCCAaacaa                                                      <  2:63770/251‑1 (MQ=255)
                                                                                                                                                                                                           ttaatgtgtCCATGATGCCTCGGTCCAGATTGTCGATCTAATAATTTTCCATATGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTTTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTTCCAATCTACGTTAAATTAGCTTCGGGaa                                  <  2:392769/242‑1 (MQ=255)
                                                                                                                                                                                                                 cTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCAc                           >  1:608992/1‑251 (MQ=255)
                                                                                                                                                                                                                 cTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGAATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCAc                           >  1:391607/1‑251 (MQ=255)
                                                                                                                                                                                                                   tacAGGATGCCACGGTCCAGATTGTAGATCAAATAATTTTCCATAGGATTTTCTTATGCGGAGTGATGATTCATTCTAATTTAGCCTTCTGTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTAGTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGAAGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACt                          <  2:559993/248‑1 (MQ=255)
                                                                                                                                                                                                                    ccAGGATTCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTTCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACttt                        <  2:572019/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTAGTCAACTgaacgac      >  2:204080/1‑245 (MQ=255)
                                                                                                                                                                                                                                       aTTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGCGTATAAAAAATGAAAAACGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAAc     >  1:157528/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTc  >  2:200275/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
GCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATC‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTG‑AGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACA‑GGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTC  >  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: ‑