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.C97.1% 234.5 / ‑0.5 68coding (518/3594 nt)evgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (33/33);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

GCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAG  >  NC_000913/2484760‑2485023
                                                                                                                                    |                                                                                                                                    
gccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccccccc                                                                                                                               <  2:24239/139‑2 (MQ=255)
   tctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                            <  2:567404/139‑5 (MQ=255)
   tctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                            >  1:1517/1‑135 (MQ=255)
    ctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccagac                                                                                                                           <  2:1517/139‑6 (MQ=255)
    ctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccagac                                                                                                                           >  1:1057678/1‑134 (MQ=255)
     ttAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACg                                                                                                                           <  2:993098/139‑1 (MQ=255)
          aatGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTa                                                                                                                      <  1:646214/139‑1 (MQ=255)
            tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACg                                                                                                                           <  1:1009796/132‑1 (MQ=255)
            tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACg                                                                                                                           >  2:1009796/1‑132 (MQ=255)
              aCATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtc                                                                                                                  >  1:906314/1‑139 (MQ=255)
               cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtca                                                                                                                 >  2:190227/1‑139 (MQ=255)
               cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtca                                                                                                                 >  2:272485/1‑139 (MQ=255)
               cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtca                                                                                                                 >  2:1149862/1‑139 (MQ=255)
                  tgctgcAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtcatca                                                                                                              >  1:1186610/1‑139 (MQ=255)
                    ctgcAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAAACGACGTAATTCatcaat                                                                                                            >  2:460362/1‑139 (MQ=255)
                       cAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTa                                                                                                                      <  2:319168/126‑1 (MQ=255)
                       cAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTa                                                                                                                      >  1:319168/1‑126 (MQ=255)
                              accaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                                  >  1:842849/1‑139 (MQ=255)
                              accaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                                  >  1:1075244/1‑139 (MQ=255)
                                caTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAg                                                                                                <  2:1075244/139‑1 (MQ=255)
                                 aTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                               <  1:272485/139‑1 (MQ=255)
                                   tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                             >  2:876131/1‑139 (MQ=255)
                                                tCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCttt                                                                                >  2:989514/1‑139 (MQ=255)
                                                tCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                            >  1:31366/1‑90 (MQ=255)
                                                tCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                            <  2:31366/94‑5 (MQ=255)
                                                tCCGGCGCTGGTAACCACCCTTCACGACACAATGAGACCGCTTACCTCATCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCttt                                                                                >  1:789525/1‑139 (MQ=255)
                                                   ggCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                             <  2:516654/139‑1 (MQ=255)
                                                   ggCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                             <  1:1149862/139‑1 (MQ=255)
                                                       cTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAAt                                                                         >  1:330841/1‑139 (MQ=255)
                                                             accaccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                   <  1:386880/139‑1 (MQ=255)
                                                             accaccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                   <  2:330841/139‑1 (MQ=255)
                                                                    ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAt                                                            <  2:1186610/139‑1 (MQ=255)
                                                                    ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTa                                                          >  1:442434/1‑139 (MQ=255)
                                                                         gACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                       <  2:842849/139‑1 (MQ=255)
                                                                          aCTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                            <  2:296230/68‑5 (MQ=255)
                                                                          aCTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                            >  1:296230/1‑64 (MQ=255)
                                                                            tCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATcc                                                    >  1:947013/1‑139 (MQ=255)
                                                                                   gACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTaa                                                                                                                     >  2:783990/1‑67 (MQ=255)
                                                                                   gACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGCTaa                                                                                                                     <  1:783990/67‑1 (MQ=255)
                                                                                     ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCt                                           <  1:190227/139‑1 (MQ=255)
                                                                                        cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAtt                                                                                        >  1:893820/1‑91 (MQ=255)
                                                                                        cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAtt                                                                                        <  2:893820/91‑1 (MQ=255)
                                                                                        cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTTAGCTggg                                        >  1:629067/1‑139 (MQ=255)
                                                                                        cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTggg                                        <  1:876131/139‑1 (MQ=255)
                                                                                         ttACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGcac                                     <  2:442434/139‑1 (MQ=255)
                                                                                             cTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACaa                                   >  1:862302/1‑139 (MQ=255)
                                                                                                aaccaaaCCAGTTAATATTGCTCGGGTAACAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCCACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGa                                <  1:230798/135‑1 (MQ=255)
                                                                                                aCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGa                                >  2:1168770/1‑139 (MQ=255)
                                                                                                 ccAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                               <  2:862302/139‑1 (MQ=255)
                                                                                                 ccAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                               <  1:989514/139‑1 (MQ=255)
                                                                                                 ccAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                               <  1:1168770/139‑1 (MQ=255)
                                                                                                  cAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCa                                                                                                    <  1:235682/69‑1 (MQ=255)
                                                                                                  cAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCa                                                                                                    >  2:235682/1‑69 (MQ=255)
                                                                                                        cAGGAAATATTACTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTa                        <  2:789525/139‑1 (MQ=255)
                                                                                                          gTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                      >  1:190076/1‑139 (MQ=255)
                                                                                                          gTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                      <  2:1057678/139‑1 (MQ=255)
                                                                                                                atTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTGGCATCCGTCTCAGCTGGGCACAATGAt                               >  2:951031/1‑124 (MQ=255)
                                                                                                                atTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTGGCATCCGTCTCAGCTGGGCACAATGAt                               <  1:951031/124‑1 (MQ=255)
                                                                                                                atTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                               >  2:664574/1‑124 (MQ=255)
                                                                                                                atTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                               <  1:664574/124‑1 (MQ=255)
                                                                                                                  tGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                       >  1:1144343/1‑98 (MQ=255)
                                                                                                                  tGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                       <  2:1144343/98‑1 (MQ=255)
                                                                                                                  tGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACt                           <  1:617199/126‑1 (MQ=255)
                                                                                                                  tGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACt                           >  2:617199/1‑126 (MQ=255)
                                                                                                                  tGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTaa              >  2:1158025/1‑139 (MQ=255)
                                                                                                                    cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                                  <  1:265011/85‑1 (MQ=255)
                                                                                                                    cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                                  >  2:265011/1‑85 (MQ=255)
                                                                                                                     tCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATATATATCAGGCATTAGCa                                                       >  2:750589/1‑95 (MQ=255)
                                                                                                                     tCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAACAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                       <  1:750589/95‑1 (MQ=255)
                                                                                                                       gggTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTa                        >  1:1191858/1‑124 (MQ=255)
                                                                                                                       gggTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTa                        <  2:1191858/124‑1 (MQ=255)
                                                                                                                         gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAg                                             <  2:1036219/101‑1 (MQ=255)
                                                                                                                         gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAg                                             >  1:1036219/1‑101 (MQ=255)
                                                                                                                         gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                      >  1:2500/1‑124 (MQ=255)
                                                                                                                         gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                      <  2:2500/124‑1 (MQ=255)
                                                                                                                         gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCAtt       >  1:914270/1‑139 (MQ=255)
                                                                                                                              aaaTTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACcag  <  2:190076/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAG  >  NC_000913/2484760‑2485023

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