Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,468,115 A→G 6.4% R87R (AGA→AGG)  gtrA → CPS‑53 (KpLE1) prophage; bactoprenol‑linked glucose translocase/flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,468,1150AG 6.4% 185.5 / 3.1 64R87R (AGA→AGG) gtrACPS‑53 (KpLE1) prophage; bactoprenol‑linked glucose translocase/flippase
Reads supporting (aligned to +/- strand):  ref base A (34/24);  new base G (2/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CTTCGCAGGTTTCGTTGTGGCTGTGAGCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATA  >  NC_000913/2467968‑2468261
                                                                                                                                                   |                                                                                                                                                  
cTTCGCATGTTTCGTTGTGGCTGTGAGCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGAt                                                                                                                                                   <  6:186926/149‑1 (MQ=255)
 ttCGCAGGTTTCGTTGTGGCTGTGAGCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATg                                                                                                                                                  >  8:141682/1‑149 (MQ=255)
      aGGTTTCGTTGTGGCTGTGAGCTTTAGCTTCTTCGCGAATGCAAAAATCACATTCAAGGCATCGACTACAACGATGCGCAACATGCTATATGTTGGGTTCATGGGGACACTGCGTGCTACTCTTGGATGGCCTGCTGATACATGCGCAc                                                                                                                                             >  5:309392/1‑149 (MQ=255)
          ttCGTTGTGGCTGTGAGCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTcc                                                                                                                                         >  1:379437/1‑149 (MQ=255)
                gtgGCTGTGAGCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCga                                                                                                                                    <  8:126836/148‑1 (MQ=255)
                       tgAGCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATaa                                                                                                                              <  1:110643/147‑1 (MQ=255)
                       tgAGCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAAct                                                                                                                            >  3:244309/1‑149 (MQ=255)
                        gAGCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAActc                                                                                                                           <  4:67689/149‑1 (MQ=255)
                          gCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCtt                                                                                                                         <  8:174431/149‑1 (MQ=255)
                                   tcttcGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTtc                                                                                                                >  5:272704/1‑149 (MQ=255)
                                     ttcGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTtctc                                                                                                              >  6:117307/1‑149 (MQ=255)
                                        gcgAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGc                                                                                                           >  7:121955/1‑149 (MQ=255)
                                                aaaaTTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCa                                                                                                                     >  4:286759/1‑131 (MQ=255)
                                                aaaaTTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCa                                                                                                                     <  3:286759/131‑1 (MQ=255)
                                                aaaaTTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGcc                                                                                                   <  2:379437/149‑1 (MQ=255)
                                                aaaaTTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGcc                                                                                                   >  7:384732/1‑149 (MQ=255)
                                                  aaTTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTg                                                                                                 <  4:244309/149‑1 (MQ=255)
                                                    ttCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGgt                                                                                               <  6:272704/149‑1 (MQ=255)
                                                    ttCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGgt                                                                                               >  7:187505/1‑149 (MQ=255)
                                                         aTTCAAGGCATCGACTACAACGATGAGCTACATTCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGTCTGCTGATAGCTGCGCACTTCCCACGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCg                                                                                          >  4:263139/1‑149 (MQ=255)
                                                          ttCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCgg                                                                                         >  8:232656/1‑149 (MQ=255)
                                                             aaGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGttt                                                                                      >  3:234281/1‑149 (MQ=255)
                                                              aGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTc                                                                                     <  3:167332/149‑1 (MQ=255)
                                                               ggCATCGACTACAACGCTGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGGTGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGcc                                                                                                   >  1:380315/1‑134 (MQ=255)
                                                               ggCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGGTGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGcc                                                                                                   <  2:380315/134‑1 (MQ=255)
                                                               ggCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGGTGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGcc                                                                                                   <  2:374934/134‑1 (MQ=255)
                                                               ggCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGGTGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGcc                                                                                                   >  1:374934/1‑134 (MQ=255)
                                                               ggCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCg                                                                                    >  8:86351/1‑149 (MQ=255)
                                                                  aTCGACTACAACGATGCGCTATATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTc                                                                                  >  8:74251/1‑148 (MQ=255)
                                                                       cTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTc                                                                             >  4:256773/1‑148 (MQ=255)
                                                                          cAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCaaa                                                                          <  7:141682/148‑1 (MQ=255)
                                                                          cAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCaaa                                                                          <  7:86351/148‑1 (MQ=255)
                                                                           aaCGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGt                                                                        >  7:276238/1‑149 (MQ=255)
                                                                           aaCGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGt                                                                        <  8:121955/149‑1 (MQ=255)
                                                                                       cATGCTATATGTTGGGTTCATGGGGACAGTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTa                                                            <  1:38373/149‑1 (MQ=255)
                                                                                        aTGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGt                                                                                   <  6:61361/125‑1 (MQ=255)
                                                                                        aTGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGt                                                                                   >  5:61361/1‑125 (MQ=255)
                                                                                        aTGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAg                                                           <  7:232656/149‑1 (MQ=255)
                                                                                        aTGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAg                                                           >  2:227568/1‑149 (MQ=255)
                                                                                           cTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAggg                                                         <  1:227568/148‑1 (MQ=255)
                                                                                                  ttGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGaaa                                                 <  8:384732/149‑1 (MQ=255)
                                                                                                   tGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAAt                                                >  2:320986/1‑149 (MQ=255)
                                                                                                    gggTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTGAGGGATGCGAAAt                                                >  1:302854/1‑148 (MQ=255)
                                                                                                      ggtcctGGGGCCACTGCGGGCTACTGCTCGATGGGCTGCTGCCAGCGGCCCCCTTCCCCCGAGGCTAACCCTTGCCCCCTTCTCCGCCACCAGCCTGGTGGGCGCTCTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGaa                                             <  2:269372/144‑1 (MQ=255)
                                                                                                             gggACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATAtct                                      >  5:33044/1‑149 (MQ=255)
                                                                                                              ggACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATAtctc                                     <  3:256773/149‑1 (MQ=255)
                                                                                                               gACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATctct                                    >  7:32462/1‑149 (MQ=255)
                                                                                                                acagTGAGTGTTACTGTTGGATGGGCTGCTTATAGATGCGCACTTCCCCCTATTATAACTCTTGTCAGCTTCTCCTCCATCAGCCTGGTGTGCGGTTTCGTCTATTCATAGTGCATTGTGTTTAGGGATGCGAAATGAAGTTATCTGtt                                   <  6:309392/149‑1 (MQ=255)
                                                                                                                acacTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCtt                                   >  7:360878/1‑149 (MQ=255)
                                                                                                                 cacTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGTGATGCGAAATGAAGATATCTCTTg                                  <  4:234281/149‑1 (MQ=255)
                                                                                                                  acTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGt                                 <  8:187505/149‑1 (MQ=255)
                                                                                                                  acTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGCGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTg                                  >  4:145453/1‑148 (MQ=255)
                                                                                                                    tGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAg                               >  5:111018/1‑149 (MQ=255)
                                                                                                                          cTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGtgtg                                                                                            >  4:238499/1‑82 (MQ=255)
                                                                                                                           tACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGt                        >  8:336047/1‑149 (MQ=255)
                                                                                                                                ttGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTc                    >  7:313367/1‑148 (MQ=255)
                                                                                                                                  ggATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAAt                 >  2:352326/1‑149 (MQ=255)
                                                                                                                                          tgctgATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATgaagaag          >  1:81950/1‑148 (MQ=255)
                                                                                                                                              gATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATgaagaa           <  3:166300/143‑1 (MQ=255)
                                                                                                                                              gATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATgaagaa           >  4:166300/1‑143 (MQ=255)
                                                                                                                                              gATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGc      >  1:354424/1‑148 (MQ=255)
                                                                                                                                              gATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGGTTCGTCGATTCAAAGTTCATTGTCTTTAGGGCGGCGAAAGGCAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGc      >  1:359861/1‑148 (MQ=255)
                                                                                                                                               aTAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGa    <  1:320986/149‑1 (MQ=255)
                                                                                                                                                 aGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATa  >  4:51091/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTTCGCAGGTTTCGTTGTGGCTGTGAGCTTTAGCTTCTTCGCGAATGCAAAATTCACATTCAAGGCATCGACTACAACGATGCGCTACATGCTATATGTTGGGTTCATGGGGACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATA  >  NC_000913/2467968‑2468261

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