Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 594,105 A→C N222K (AAT→AAG)  cusS ← copper‑sensing histidine kinase in two‑component regulatory system with CusR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913594,1050AC99.5% 695.4 / ‑5.7 194N222K (AAT→AAG) cusScopper‑sensing histidine kinase in two‑component regulatory system with CusR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (94/98);  total (95/99)
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.81e-01

Only 100 of 194 total aligned reads displayed.

ACGCTCGATCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGC  >  NC_000913/594009‑594200
                                                                                                |                                                                                               
aCGCTCGATCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTc                                                                                              <  2:2419100/100‑1 (MQ=255)
      gATCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGAtt                                                                                        >  1:3101550/1‑100 (MQ=255)
       aTCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGAttt                                                                                       <  2:3859025/100‑1 (MQ=255)
        tCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTg                                                                                      <  2:970831/100‑1 (MQ=255)
        tCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTg                                                                                      <  2:1759053/100‑1 (MQ=255)
          atatGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGAc                                                                                    <  2:3081688/100‑1 (MQ=255)
          atatGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGAc                                                                                    >  1:433522/1‑100 (MQ=255)
          atatGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGAc                                                                                    <  1:2192164/100‑1 (MQ=255)
           tatGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACg                                                                                   <  1:552025/100‑1 (MQ=255)
            atGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACgg                                                                                  <  1:1394494/100‑1 (MQ=255)
             tGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGc                                                                                 <  2:2710708/100‑1 (MQ=255)
               gTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTg                                                                               >  2:120252/1‑100 (MQ=255)
               gTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTg                                                                               <  2:1367654/100‑1 (MQ=255)
                ttGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGa                                                                              <  2:4351971/100‑1 (MQ=255)
                   aaCGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGc                                                                           >  1:3427686/1‑100 (MQ=255)
                     cGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCTAACGTCGAGCACTTTCAAGGGAATCTTCTGGATTAGACGGCTGAAGCTg                                                                         >  2:4670369/1‑100 (MQ=255)
                      gACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGc                                                                        >  1:3049985/1‑100 (MQ=255)
                       ccAGTACCAGCTGTTCCAGCTCAAGGGGCACGGTCAGCGGGTCGAGGCCACCGTCGACATCTTTCGAGGTAATCCTCTGGATTTGACGGCCGACGCTGCg                                                                       >  2:2669902/2‑100 (MQ=255)
                       aCAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCg                                                                       <  1:4526448/100‑1 (MQ=255)
                        cAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCgg                                                                      <  1:3075686/100‑1 (MQ=255)
                          gTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGAt                                                                    >  1:3837264/1‑100 (MQ=255)
                           tACCGGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGTCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACCGCTGACGCTGCGGATc                                                                   <  2:3644514/100‑1 (MQ=255)
                           tACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATc                                                                   <  1:2826796/100‑1 (MQ=255)
                           tACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATc                                                                   <  2:4844235/100‑1 (MQ=255)
                            aCCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCg                                                                  <  1:2100982/100‑1 (MQ=255)
                            aCCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCg                                                                  >  1:2505995/1‑100 (MQ=255)
                            aCCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCg                                                                  >  1:2292250/1‑100 (MQ=255)
                             ccAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCgg                                                                 <  2:785716/100‑1 (MQ=255)
                             ccAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCgg                                                                 >  2:3574484/1‑100 (MQ=255)
                              cAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGgc                                                                >  2:1434924/1‑100 (MQ=255)
                               aGCTGTTCCAGCTCCATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGgcg                                                               >  2:4119257/1‑100 (MQ=255)
                               aGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGgcg                                                               <  2:1095134/100‑1 (MQ=255)
                               aGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGgcg                                                               <  1:3580678/100‑1 (MQ=255)
                               aGCTGTTCCAGCTCAATGGACACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGgcg                                                               >  1:890730/1‑100 (MQ=255)
                                gCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGgcgc                                                              <  2:2019014/100‑1 (MQ=255)
                                 cTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGcgcg                                                             <  1:186058/100‑1 (MQ=255)
                                 cTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGcgcg                                                             <  1:1538810/100‑1 (MQ=255)
                                   gTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTg                                                           <  2:1653455/100‑1 (MQ=255)
                                   gTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTg                                                           >  2:606229/1‑100 (MQ=255)
                                   gTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTg                                                           >  2:2284086/1‑100 (MQ=255)
                                   gTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTg                                                           >  1:3250893/1‑100 (MQ=255)
                                     tCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGAc                                                         >  1:728821/1‑100 (MQ=255)
                                     tCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGAc                                                         <  2:4871710/100‑1 (MQ=255)
                                       cAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCt                                                       >  2:4472040/1‑100 (MQ=255)
                                         gCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCttt                                                     >  1:1001319/1‑100 (MQ=255)
                                         gCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCttt                                                     >  1:3953877/1‑100 (MQ=255)
                                         gCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCttt                                                     <  2:1471784/100‑1 (MQ=255)
                                          cTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTa                                                    <  2:3458203/100‑1 (MQ=255)
                                           tCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTAt                                                   <  2:4595471/100‑1 (MQ=255)
                                            cAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATg                                                  <  2:3873015/100‑1 (MQ=255)
                                              aTGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTa                                                >  2:3390469/1‑100 (MQ=255)
                                              aTGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTa                                                >  1:3616590/1‑100 (MQ=255)
                                               tGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTAc                                               >  1:5342775/1‑100 (MQ=255)
                                                 ggCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCg                                             >  1:5440239/1‑100 (MQ=255)
                                                 ggCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCg                                             >  1:5155603/1‑100 (MQ=255)
                                                   cACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGcc                                           >  2:1182000/1‑100 (MQ=255)
                                                    aCGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCa                                          >  1:850750/1‑100 (MQ=255)
                                                    aCGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCa                                          >  2:4847070/1‑100 (MQ=255)
                                                      ggTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAAc                                        <  1:2773888/100‑1 (MQ=255)
                                                      ggTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAAc                                        >  1:155185/1‑100 (MQ=255)
                                                      ggTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAAc                                        <  2:2501452/100‑1 (MQ=255)
                                                       gTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACa                                       >  1:4182486/1‑100 (MQ=255)
                                                        tCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAg                                      <  2:2138072/100‑1 (MQ=255)
                                                        tCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAg                                      <  2:5348002/100‑1 (MQ=255)
                                                        tCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAg                                      <  2:4777408/100‑1 (MQ=255)
                                                        tCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAg                                      <  1:2028968/100‑1 (MQ=255)
                                                          tGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTa                                    >  1:5566536/1‑100 (MQ=255)
                                                          tGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTa                                    >  2:5438503/1‑100 (MQ=255)
                                                           gCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTAc                                   <  2:900018/100‑1 (MQ=255)
                                                           gCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTAc                                   <  1:4466114/100‑1 (MQ=255)
                                                                tCGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATaa                              <  1:1765039/100‑1 (MQ=255)
                                                                 cGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATaaa                             >  2:3925441/1‑100 (MQ=255)
                                                                 cGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATaaa                             <  1:1049104/100‑1 (MQ=255)
                                                                 cGAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTCGATGTGACGGCTGACGCCGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATaaa                             >  2:5158216/1‑100 (MQ=255)
                                                                  gAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAg                            >  1:796583/1‑100 (MQ=255)
                                                                  gAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAg                            >  2:295233/1‑100 (MQ=255)
                                                                  gAGGCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAg                            >  1:2660062/1‑100 (MQ=255)
                                                                    ggCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGAc                          <  1:3770394/100‑1 (MQ=255)
                                                                     gCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACg                         <  2:1899822/100‑1 (MQ=255)
                                                                     gCGAACGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACg                         >  1:460735/1‑100 (MQ=255)
                                                                         aCGTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCa                     >  2:1992822/1‑100 (MQ=255)
                                                                           gTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGGTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGAGCAgg                   <  2:3534798/100‑1 (MQ=255)
                                                                           gTCGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAgg                   >  1:4018264/1‑100 (MQ=255)
                                                                             cGAGATCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGAt                 <  1:1423194/100‑1 (MQ=255)
                                                                                gaTCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCt              <  1:244454/100‑1 (MQ=255)
                                                                                 aTCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTg             >  1:5051545/1‑100 (MQ=255)
                                                                                 aTCTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTg             <  2:5084748/100‑1 (MQ=255)
                                                                                  tCTTTCGAGGTAATCTGCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGa            <  2:3602830/100‑1 (MQ=255)
                                                                                   cTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGAt           >  1:1381435/1‑100 (MQ=255)
                                                                                   cTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGAt           >  2:3281935/1‑100 (MQ=255)
                                                                                   cTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGAt           >  2:4911589/1‑100 (MQ=255)
                                                                                   cTTTCGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGAt           <  1:3422246/100‑1 (MQ=255)
                                                                                       cGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTAcc       >  1:1721848/1‑100 (MQ=255)
                                                                                       cGAGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTAcc       >  1:548931/1‑100 (MQ=255)
                                                                                       cGAGGTAAGCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTAcc       >  2:5270589/1‑100 (MQ=255)
                                                                                         aGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTAccag     >  1:4789112/1‑98 (MQ=255)
                                                                                         aGGTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGa     >  1:1175683/1‑100 (MQ=255)
                                                                                          ggTAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGAt    <  2:3781852/100‑1 (MQ=255)
                                                                                            tAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGc  >  2:4086594/1‑100 (MQ=255)
                                                                                            tAATCTTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGc  >  2:1858099/1‑100 (MQ=255)
                                                                                                |                                                                                               
ACGCTCGATCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGC  >  NC_000913/594009‑594200

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