Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,188,434 A→G 46.1% V448A (GTA→GCA)  phoQ ← sensory histidine kinase in two‑component regulatory system with PhoP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,188,4340AG46.1% 17.5 / 169.4 128V448A (GTA→GCA) phoQsensory histidine kinase in two‑component regulatory system with PhoP
Reads supporting (aligned to +/- strand):  ref base A (31/38);  new base G (26/33);  total (57/71)
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.90e-01

Only 100 of 129 total aligned reads displayed.

ATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCA  >  NC_000913/1188337‑1188532
                                                                                                 |                                                                                                  
aTGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTAcc                                                                                                  >  1:422170/1‑100 (MQ=255)
 tGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTATCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCg                                                                                                 >  1:623063/1‑100 (MQ=255)
  gCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTGCGCGGGCTACCGc                                                                                                >  2:2463582/1‑100 (MQ=255)
   cTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGcc                                                                                               >  1:3272765/1‑100 (MQ=255)
    tGGCGACCAAACATCGCCTCCATCCGCGCACCGACGAGCATGTTCTCCCCGGCGACGATTTTACCCTCATATTTCTCGGTGATTTCGCACGCTACCGCCa                                                                                              <  2:2013860/100‑1 (MQ=255)
      gCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGc                                                                                            <  2:3644147/100‑1 (MQ=255)
      gCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGc                                                                                            >  2:2610149/1‑100 (MQ=255)
       cGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGcc                                                                                           >  2:2085636/1‑100 (MQ=255)
       cGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGcc                                                                                           <  1:2346183/100‑1 (MQ=255)
       cGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGcc                                                                                           <  1:780430/100‑1 (MQ=255)
        gACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGccc                                                                                          <  1:3527609/100‑1 (MQ=255)
        gACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGccc                                                                                          <  2:2224389/100‑1 (MQ=255)
         aCCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCt                                                                                         <  1:1887850/100‑1 (MQ=255)
          ccAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTa                                                                                        >  2:3252741/1‑100 (MQ=255)
          ccAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACCATTTTACCCTCCTATTGCCCGGTGATTTCGCGGGCTGCCGCCAGCCCt                                                                                         >  2:1176075/1‑99 (MQ=255)
           cAAAATTCACCTCCATCCGTGCACCGCCCAGCATGCTCTCTCTGACGCCGGTTTTACCGTCATATTACCCGGTGATTTCGCGAGATACCGCCAGccc                                                                                          >  2:2518260/1‑97 (MQ=255)
            aaaaaTCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTaca                                                                                      <  2:3632714/100‑1 (MQ=255)
            aaaaaTCACCTCCATCCGCACACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTaca                                                                                      <  1:833847/100‑1 (MQ=255)
            aaaaaTCACCTCCATCCGCACACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTaca                                                                                      <  1:833997/100‑1 (MQ=255)
                  cACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGc                                                                                <  1:2416477/100‑1 (MQ=255)
                   aCCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGcc                                                                               >  1:453068/1‑100 (MQ=255)
                      tCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAg                                                                            >  1:1677900/1‑100 (MQ=255)
                        cATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGGTTTCTCGGGCTGCCGCCAGCCCTACACCTTGCCCTggg                                                                          >  1:2665750/1‑100 (MQ=255)
                        cATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAggg                                                                          >  2:1252430/1‑100 (MQ=255)
                         aTCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGc                                                                         <  2:475329/100‑1 (MQ=255)
                          tCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCg                                                                        <  2:5468/100‑1 (MQ=255)
                          tCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCg                                                                        <  2:3330765/100‑1 (MQ=255)
                          tCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCg                                                                        <  2:1546900/100‑1 (MQ=255)
                           ccGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGt                                                                       >  2:3598461/1‑100 (MQ=255)
                           ccGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGt                                                                       >  1:2318020/1‑100 (MQ=255)
                            cgcgCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTa                                                                      >  2:4010821/1‑100 (MQ=255)
                            cgcgCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTa                                                                      <  2:1265089/100‑1 (MQ=255)
                            cgcgCACCGCCCAGCATGCTCTCTCCGGCAACGATTTTACCCTCATATTGCTCGGTGCTTTCGCGGGCTACCGCCAGCCCTAAACCTTGCCCCGGGCGt                                                                       >  2:2909998/1‑99 (MQ=255)
                            cgcgCACCGCCCAGCATGCGCTCTCCGGCGACGGTTTTACCCTCATTTTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCg                                                                        <  2:3307179/98‑1 (MQ=255)
                              cgcACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTaaa                                                                    <  1:3076272/100‑1 (MQ=255)
                                  ccGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTAt                                                                <  2:3010170/100‑1 (MQ=255)
                                  ccGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTAt                                                                >  1:1541534/1‑100 (MQ=255)
                                   cGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATc                                                               >  1:614991/1‑100 (MQ=255)
                                   cGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATc                                                               <  1:3961456/100‑1 (MQ=255)
                                   cGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATc                                                               <  1:670299/100‑1 (MQ=255)
                                    gCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCa                                                              >  1:2329210/1‑100 (MQ=255)
                                     cctaacATGCTCTCTCCGGCGACGATTTTACCATCATATTGCTCGGTTATTTCGCGGGCTACCGCCAGCCCTCCACCTTGCCCAGGGCGTAAAGTCTCa                                                              >  2:2417728/6‑99 (MQ=255)
                                     cccAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCaa                                                             <  1:154354/100‑1 (MQ=255)
                                       cAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAAcc                                                           <  2:2919725/100‑1 (MQ=255)
                                         gCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCg                                                         >  1:2036609/1‑100 (MQ=255)
                                             gCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGtt                                                       >  2:663719/1‑98 (MQ=255)
                                             gCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTg                                                      <  1:2463596/99‑1 (MQ=255)
                                              ctctctCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGAc                                                    <  1:3267130/100‑1 (MQ=255)
                                              ctctctCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGAc                                                    >  1:2625241/1‑100 (MQ=255)
                                               tctctcCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGACGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGAcc                                                   <  2:2861736/100‑1 (MQ=255)
                                                ctctcCGGCGACGATTTTACCCTCATATTACTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCa                                                  <  1:1340754/100‑1 (MQ=255)
                                                 nctcCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCAc                                                 >  1:2057357/2‑100 (MQ=255)
                                                    ccGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGt                                              >  1:2125901/1‑100 (MQ=255)
                                                      ggCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCg                                            >  1:1353537/1‑100 (MQ=255)
                                                      ggCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCg                                            >  2:850567/1‑100 (MQ=255)
                                                      ggCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCg                                            <  1:1085665/100‑1 (MQ=255)
                                                       gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGa                                           <  2:1850637/100‑1 (MQ=255)
                                                           cgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAAt                                       <  1:1905372/100‑1 (MQ=255)
                                                            gaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATg                                      <  2:812987/100‑1 (MQ=255)
                                                            gaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATg                                      <  2:3362748/100‑1 (MQ=255)
                                                              ttttACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGAc                                    <  1:344362/100‑1 (MQ=255)
                                                               tttACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGAcc                                   <  2:3871628/100‑1 (MQ=255)
                                                               tttACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCATTGACCACGGTCGAAAATGAcc                                   <  2:1021472/100‑1 (MQ=255)
                                                                 tACCCTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACct                                  >  2:1606882/1‑99 (MQ=255)
                                                                   cccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCCGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCtctc                               >  2:2875256/1‑100 (MQ=255)
                                                                     cTCATTTTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCGTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgc                             <  2:4096830/100‑1 (MQ=255)
                                                                     cTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgc                             >  1:1223983/1‑100 (MQ=255)
                                                                     cTCATATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgc                             >  1:3195470/1‑100 (MQ=255)
                                                                     cTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgc                             >  1:242172/1‑100 (MQ=255)
                                                                      tCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgct                            >  1:147892/1‑100 (MQ=255)
                                                                       cATATTGCTCGGTGATTTCGCGGGGTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAACGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgctt                           <  2:227115/100‑1 (MQ=255)
                                                                       cATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgctt                           <  1:1252438/100‑1 (MQ=255)
                                                                        atatTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCAGGGTCGAAAATGACCTCTCgcttg                          <  1:3500691/100‑1 (MQ=255)
                                                                         tatTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgc                         >  1:3843921/1‑100 (MQ=255)
                                                                         tatTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgc                         >  2:2946030/1‑100 (MQ=255)
                                                                         tatTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgc                         >  2:2760283/1‑100 (MQ=255)
                                                                         tatTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgc                         >  2:2993220/1‑100 (MQ=255)
                                                                          atTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgct                        <  1:1366476/100‑1 (MQ=255)
                                                                          atTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgct                        >  2:294281/1‑100 (MQ=255)
                                                                           ttGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgctt                       <  1:3283897/100‑1 (MQ=255)
                                                                           ttGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgctt                       <  1:2610165/100‑1 (MQ=255)
                                                                            tGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTa                      >  1:3998680/1‑100 (MQ=255)
                                                                             gCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTaa                     >  2:4071662/1‑100 (MQ=255)
                                                                             gCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTaa                     >  1:3072926/1‑100 (MQ=255)
                                                                              cTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTaa                     >  2:1006077/1‑99 (MQ=255)
                                                                              cTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAAt                    >  1:718822/1‑100 (MQ=255)
                                                                               tCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATg                   <  1:2085649/100‑1 (MQ=255)
                                                                               tCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATg                   <  1:3228283/100‑1 (MQ=255)
                                                                                cGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATgg                  >  1:1927846/1‑100 (MQ=255)
                                                                                cGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATgg                  <  2:614985/100‑1 (MQ=255)
                                                                                cGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATgg                  <  1:1233709/100‑1 (MQ=255)
                                                                                cGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATgg                  >  1:955324/1‑100 (MQ=255)
                                                                                  gTGATTTCGCGGGCTACCACCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACAACGGTCGAAAATGACATCTCGCTTGCTTAATGGaa                >  2:1996980/1‑100 (MQ=255)
                                                                                    gATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATa              >  2:1333804/1‑100 (MQ=255)
                                                                                        tCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCgg          <  1:2256380/100‑1 (MQ=255)
                                                                                          gcgGGCTGCCGCCAGCCCTGCACCTGGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAGGATGACCGCTCGCTTGCTTAATTGAATACCgggg        <  2:231318/100‑1 (MQ=255)
                                                                                          gcgGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCgggg        <  1:540187/100‑1 (MQ=255)
                                                                                           cgGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGc       <  1:3770358/100‑1 (MQ=255)
                                                                                             ggCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACGCGTTGACCACGGCCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCca     >  2:1031522/1‑100 (MQ=255)
                                                                                                tACCTCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCcatca  <  2:3309409/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
ATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCA  >  NC_000913/1188337‑1188532

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