Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,484,897 (C)6→7 coding (524/3594 nt) evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,484,8911.C98.7% 261.4 / ‑3.7 76coding (518/3594 nt)evgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (36/39);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑C‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACAT  >  NC_000913/2484756‑2485014
                                                                                                                                         |                                                                                                                            
cgccgcCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c‑cc                                                                                                                          >  1:125719/1‑139 (MQ=255)
    gccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAAACACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c‑cccccc                                                                                                                      <  1:500419/139‑2 (MQ=255)
              aatGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTa                                                                                                             <  1:719635/139‑1 (MQ=255)
                tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTa                                                                                                                               <  2:958720/120‑1 (MQ=255)
                tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTa                                                                                                                               >  1:958720/1‑120 (MQ=255)
                tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAAt                                                                                                           >  2:164992/1‑139 (MQ=255)
                   cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATtca                                                                                                        >  2:1591758/1‑139 (MQ=255)
                   cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCACTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATtca                                                                                                        >  2:1173241/1‑139 (MQ=255)
                     ttGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c                                                                                                                             <  2:582462/116‑1 (MQ=255)
                     ttGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c                                                                                                                             >  1:582462/1‑116 (MQ=255)
                        ctgcAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCatcaat                                                                                                   >  2:990471/1‑139 (MQ=255)
                             accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCAttt                                                                                              >  1:110218/1‑139 (MQ=255)
                             accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c‑ccccccaga                                                                                                                   <  1:230026/117‑5 (MQ=255)
                             accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c‑ccccccaga                                                                                                                   >  2:230026/1‑113 (MQ=255)
                               caaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTcc                                                                                            <  1:1593193/139‑1 (MQ=255)
                                  accaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                         >  1:762017/1‑139 (MQ=255)
                                     aTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                      <  2:125719/139‑1 (MQ=255)
                                       tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                    >  1:371905/1‑139 (MQ=255)
                                       tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTAACTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                    <  1:950885/139‑1 (MQ=255)
                                        gATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaac                                                                                   >  2:277371/1‑139 (MQ=255)
                                             tcaCCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTa                                                                              >  1:225841/3‑139 (MQ=255)
                                              tACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt                                                                             >  2:1374980/1‑139 (MQ=255)
                                              tACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c                                                                                                                             <  2:1398695/91‑1 (MQ=255)
                                              tACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c                                                                                                                             >  1:1398695/1‑91 (MQ=255)
                                               aCCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAtc                                                                            <  1:1552264/139‑1 (MQ=255)
                                                ccTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                      >  2:684650/1‑128 (MQ=255)
                                                ccTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                      <  1:684650/128‑1 (MQ=255)
                                                   ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c‑ccccccagac                                                                                                                  >  1:1402937/1‑91 (MQ=255)
                                                   ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c‑ccccccagac                                                                                                                  <  2:1402937/96‑6 (MQ=255)
                                                    tCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                         <  2:736824/121‑1 (MQ=255)
                                                    tCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                         >  1:736824/1‑121 (MQ=255)
                                                      cGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                     <  2:762017/139‑1 (MQ=255)
                                                      cGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                     <  2:110218/139‑1 (MQ=255)
                                                              gTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                     >  2:127507/1‑131 (MQ=255)
                                                              gTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                     <  1:127507/131‑1 (MQ=255)
                                                              gTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTta                                                             <  1:1591758/139‑1 (MQ=255)
                                                                     cccTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATtca                                                                                                        >  2:1550086/1‑89 (MQ=255)
                                                                     cccTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATtca                                                                                                        <  1:1550086/89‑1 (MQ=255)
                                                                      ccTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTa                                                                                                             >  2:1581358/1‑83 (MQ=255)
                                                                      ccTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTa                                                                                                             <  1:1581358/83‑1 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACg                                                                                                                  >  1:975971/1‑76 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACg                                                                                                                  <  2:975971/76‑1 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAt                                                   >  2:578206/1‑139 (MQ=255)
                                                                        ttCACGACTCAATGCGACAGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTACCAAATTA‑CGCTCCCAGACGAGGTAATTCATCAATAATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGACAt                                                   >  1:218277/1‑139 (MQ=255)
                                                                          cACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCatcaatca                                                                                                 <  1:262572/91‑1 (MQ=255)
                                                                          cACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCatcaatca                                                                                                 >  2:262572/1‑91 (MQ=255)
                                                                             gACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              <  1:164992/139‑1 (MQ=255)
                                                                             gACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              <  2:292130/139‑1 (MQ=255)
                                                                             gACCCAATGCTACCGCCTACCGCGCCAACACCGGTAAATA‑TGGATAGGATAGCAAATTACC‑CCCCCAGACGAAGTAACTCAACAAGCATTCCCAAATGGATCTATTATCTCTGTGACAAATTTATATCAGGGATTAGCa                                              <  1:349968/139‑1 (MQ=255)
                                                                               cTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATc                                            <  2:225841/139‑1 (MQ=255)
                                                                                tCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATcc                                           >  2:550854/1‑139 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                     <  1:1370992/111‑1 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                     >  2:1370992/1‑111 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              >  1:1548852/1‑134 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              <  2:1548852/134‑1 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                       >  2:503204/1‑139 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c‑ccccccaga                                                                                                                   <  1:1358788/62‑5 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑c‑ccccccaga                                                                                                                   >  2:1358788/1‑58 (MQ=255)
                                                                                        aCCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGc                                   <  1:578206/139‑1 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAAt                                                                >  2:774916/1‑109 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAAt                                                                <  1:774916/109‑1 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCt                                  <  2:371905/139‑1 (MQ=255)
                                                                                                     tcAAAGCCAGTAAATATTTGCT‑GGGTAGCAAATGACC‑CCCCCAGACGAGGTAATTTATCAAACGTTTCCAACAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCGCAATGAt                      <  1:12151/138‑1 (MQ=255)
                                                                                                            cAGTAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTa               <  1:1374980/139‑1 (MQ=255)
                                                                                                             aGTAAATA‑TTGCTCGTGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAGAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAt              <  2:747579/139‑1 (MQ=255)
                                                                                                               tAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                         <  2:1432745/94‑1 (MQ=255)
                                                                                                               tAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                         >  1:1432745/1‑94 (MQ=255)
                                                                                                               tAAATA‑TTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGCTAATTCATCCATCATTTCCAAAAGCAGCAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTg            <  1:277371/139‑1 (MQ=255)
                                                                                                                    a‑tTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                    >  2:1323664/1‑62 (MQ=255)
                                                                                                                    a‑tTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                    <  1:1323664/62‑1 (MQ=255)
                                                                                                                    a‑tTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                       >  2:1159434/1‑107 (MQ=255)
                                                                                                                    a‑tTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                       <  1:1159434/107‑1 (MQ=255)
                                                                                                                    a‑tTGCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGgtagt       <  1:990471/139‑1 (MQ=255)
                                                                                                                        gCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                      >  2:208185/1‑57 (MQ=255)
                                                                                                                        gCTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                      <  1:208185/57‑1 (MQ=255)
                                                                                                                         cTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGca                                                                                     <  1:1202701/57‑1 (MQ=255)
                                                                                                                         cTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGca                                                                                     >  2:1202701/1‑57 (MQ=255)
                                                                                                                         cTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                    >  1:515284/1‑74 (MQ=255)
                                                                                                                         cTCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                    <  2:515284/74‑1 (MQ=255)
                                                                                                                          tCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTa                                                 >  1:1602774/1‑92 (MQ=255)
                                                                                                                          tCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTa                                                 <  2:1602774/92‑1 (MQ=255)
                                                                                                                          tCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              >  1:1185026/1‑95 (MQ=255)
                                                                                                                          tCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              <  2:1185026/95‑1 (MQ=255)
                                                                                                                          tCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAcat  >  1:619206/1‑139 (MQ=255)
                                                                                                                          tCGGGTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAcat  >  1:543422/1‑139 (MQ=255)
                                                                                                                              gTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              <  1:71497/91‑1 (MQ=255)
                                                                                                                              gTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              >  2:71497/1‑91 (MQ=255)
                                                                                                                              gTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                      <  2:730720/115‑1 (MQ=255)
                                                                                                                              gTAGCAAATTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                      >  1:730720/1‑115 (MQ=255)
                                                                                                                                   aaaTTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              >  2:881845/1‑86 (MQ=255)
                                                                                                                                   aaaTTACC‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                              <  1:881845/86‑1 (MQ=255)
                                                                                                                                         |                                                                                                                            
CGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATA‑TTGCTCGGGTAGCAAATTA‑C‑CCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACAT  >  NC_000913/2484756‑2485014

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