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.5% 241.0 / ‑3.7 68coding (518/3594 nt)evgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (33/34);  total (33/35)
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.79e-01

CGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATT  >  NC_000913/2484756‑2485018
                                                                                                                                        |                                                                                                                                
cgccgcCTCTTAATAAGGACATTGCTGCAACCAAACCATTGATAATTACCTTTCTGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccc                                                                                                                              >  2:798238/1‑139 (MQ=255)
  ccgccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccccc                                                                                                                            <  2:628840/139‑1 (MQ=255)
  ccgccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccccc                                                                                                                            >  1:1217293/1‑139 (MQ=255)
       tctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                       <  1:1317501/139‑5 (MQ=255)
              aatGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTa                                                                                                                 <  2:234787/139‑1 (MQ=255)
               atGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTaa                                                                                                                <  2:52920/139‑1 (MQ=255)
                tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAAt                                                                                                               >  1:1070209/1‑139 (MQ=255)
                   cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtca                                                                                                            >  1:718548/1‑139 (MQ=255)
                        ctgcAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                   <  1:1076743/112‑1 (MQ=255)
                        ctgcAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                   >  2:1076743/1‑112 (MQ=255)
                        ctgcAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCatcaat                                                                                                       >  1:1277371/1‑139 (MQ=255)
                             accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAc                                                                                                                                 >  2:208160/1‑108 (MQ=255)
                             accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAc                                                                                                                                 <  1:208160/108‑1 (MQ=255)
                               caaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTcc                                                                                                <  2:12167/139‑1 (MQ=255)
                                  accaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                             >  1:923661/1‑139 (MQ=255)
                                    caTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAg                                                                                           <  2:1277371/139‑1 (MQ=255)
                                     aTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                          <  2:718548/139‑1 (MQ=255)
                                      ttGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTaa                                                                                                                >  1:1196872/1‑116 (MQ=255)
                                      ttGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTaa                                                                                                                <  2:1196872/116‑1 (MQ=255)
                                       tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAGACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCCAAAg                                                                                           <  2:833152/136‑1 (MQ=255)
                                       tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTcc                                                                                                >  2:553457/1‑131 (MQ=255)
                                       tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTcc                                                                                                <  1:553457/131‑1 (MQ=255)
                                       tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAg                                                                                           >  1:833152/1‑136 (MQ=255)
                                        gATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaac                                                                                       >  2:1284514/1‑139 (MQ=255)
                                                ccTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                             >  2:45767/1‑125 (MQ=255)
                                                ccTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                             <  1:45767/125‑1 (MQ=255)
                                                   ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTcc                                                                                                <  2:1369989/119‑1 (MQ=255)
                                                   ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTcc                                                                                                >  1:1369989/1‑119 (MQ=255)
                                                   ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCtt                                                                            >  2:1024418/1‑139 (MQ=255)
                                                   ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCtt                                                                            >  2:411342/1‑139 (MQ=255)
                                                   ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCtt                                                                            >  2:950605/1‑139 (MQ=255)
                                                   ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCtt                                                                            >  1:177872/1‑139 (MQ=255)
                                                      cGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                         <  2:1217293/139‑1 (MQ=255)
                                                       ggCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                        <  1:572679/139‑1 (MQ=255)
                                                           cTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCGTTTCCAAAAGCAACAATTATCTCTTTTACAAAt                                                                    <  1:950605/139‑1 (MQ=255)
                                                           cTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAAt                                                                    <  1:1284514/139‑1 (MQ=255)
                                                                  ccaccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                             <  1:1024418/139‑1 (MQ=255)
                                                                  ccaccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                             <  2:923661/139‑1 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                   >  1:389993/1‑64 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                   <  2:389993/64‑1 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCa                                                                                               >  2:1323350/1‑99 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCa                                                                                               <  1:1323350/99‑1 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                         <  2:711528/121‑1 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                         >  1:711528/1‑121 (MQ=255)
                                                                              aCTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                             >  1:153622/1‑127 (MQ=255)
                                                                              aCTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                             <  2:153622/127‑1 (MQ=255)
                                                                              aCTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCAt                                                 <  1:411342/139‑1 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                   >  2:738576/1‑54 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                   <  1:738576/54‑1 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATc                                                >  1:733082/1‑136 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATc                                                <  2:733082/136‑1 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGt                                             >  2:17540/1‑139 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                           <  2:1070209/139‑1 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                       <  1:36167/62‑5 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                       >  2:36167/1‑58 (MQ=255)
                                                                                            cTTACCTCACCAAAATCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                          <  1:889133/84‑1 (MQ=255)
                                                                                            cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                          >  2:889133/1‑84 (MQ=255)
                                                                                            cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                         <  1:807489/101‑1 (MQ=255)
                                                                                            cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                         >  2:807489/1‑101 (MQ=255)
                                                                                             ttACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGc                                  <  2:657001/139‑1 (MQ=255)
                                                                                                ccTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                          >  1:643076/1‑80 (MQ=255)
                                                                                                ccTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                          <  2:643076/80‑1 (MQ=255)
                                                                                                         aaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGc                                  <  2:16602/127‑1 (MQ=255)
                                                                                                         aaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGc                                  >  1:16602/1‑127 (MQ=255)
                                                                                                               tAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaa                                                                                              <  1:729243/61‑1 (MQ=255)
                                                                                                               tAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaa                                                                                              >  2:729243/1‑61 (MQ=255)
                                                                                                               tAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTg                >  2:879832/1‑139 (MQ=255)
                                                                                                                      tGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                  <  1:21011/98‑1 (MQ=255)
                                                                                                                      tGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                  >  2:21011/1‑98 (MQ=255)
                                                                                                                      tGCTCGGGTAGCAAATTACCACCCCAGACGAGGTAATTCATCAATCATTTCCAAAATCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTaa         <  1:386531/139‑1 (MQ=255)
                                                                                                                       gCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGgtag            >  2:68426/1‑135 (MQ=255)
                                                                                                                       gCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGgtag            <  1:68426/135‑1 (MQ=255)
                                                                                                                       gCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAc        >  1:504654/1‑139 (MQ=255)
                                                                                                                        cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACTATGATTACTTTAtt                 >  2:814648/1‑129 (MQ=255)
                                                                                                                        cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACTATGATTACTTTAtt                 <  1:814648/129‑1 (MQ=255)
                                                                                                                        cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAca       >  1:278981/1‑139 (MQ=255)
                                                                                                                        cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAca       <  2:504654/139‑1 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAtt                                                      >  1:169966/1‑87 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAtt                                                      <  2:169966/87‑1 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAca       <  1:1142043/134‑1 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAca       >  2:1142043/1‑134 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCAtt  >  1:495946/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                
CGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATT  >  NC_000913/2484756‑2485018

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: