Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,356,775 C→T A279T (GCC→ACC)  cadB ← transport of lysine/cadaverine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,7750CT99.4% 1205.8 / ‑4.7 358A279T (GCC→ACC) cadBtransport of lysine/cadaverine
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (186/170);  total (188/170)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 372 total aligned reads displayed.

CTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAG  >  NC_000913/4356685‑4356864
                                                                                          |                                                                                         
 tttCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAggtggt                                                                                        >  2:3105743/1‑93 (MQ=255)
 tttCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAggtggt                                                                                        <  1:4510497/93‑1 (MQ=255)
 tttCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAggtggt                                                                                        >  1:1950131/1‑93 (MQ=255)
 tttCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAggtggt                                                                                        <  2:6621530/93‑1 (MQ=255)
 tttCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTggg                                                                                        >  2:5884443/1‑92 (MQ=255)
    cGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAAAATCCTCCGGGGGCCAAGAGAAGTCAGCCACGCAAGGGTGGGGaa                                                                                     >  1:2511510/1‑93 (MQ=255)
      ggAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTTAATg                                                                                   <  2:5048282/93‑1 (MQ=255)
      ggAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATg                                                                                   >  1:2940659/1‑93 (MQ=255)
        aaGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCa                                                                                 <  2:4107057/93‑1 (MQ=255)
        aaGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCa                                                                                 <  1:7136556/93‑1 (MQ=255)
          gTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGa                                                                               <  2:1010097/93‑1 (MQ=255)
          gTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGa                                                                               <  1:4437199/93‑1 (MQ=255)
          gTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGa                                                                               >  1:2121199/1‑93 (MQ=255)
             aCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAAc                                                                            <  1:7366962/93‑1 (MQ=255)
              ccGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGGAGAAAca                                                                           >  2:5963053/1‑92 (MQ=255)
              ccGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAAcc                                                                           >  2:4640140/1‑93 (MQ=255)
               cgtcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCa                                                                          <  1:7050296/93‑1 (MQ=255)
               cgtcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCa                                                                          <  1:4929291/93‑1 (MQ=255)
                 tcgtCAGCGGCTCGTGCACCGGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAgc                                                                        <  1:4364770/93‑1 (MQ=255)
                 gcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAgc                                                                        <  2:6816144/92‑1 (MQ=255)
                   gtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAgcgg                                                                      >  2:6768406/1‑93 (MQ=255)
                   gtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAgcgg                                                                      <  1:1695489/93‑1 (MQ=255)
                    ttAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGTTGGTGAATGCAGAAACCAgcggc                                                                     >  2:4315439/1‑93 (MQ=255)
                    ttAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAgcggc                                                                     <  1:192746/93‑1 (MQ=255)
                    ttAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAgcggc                                                                     >  2:557195/1‑93 (MQ=255)
                    ttAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAgcggc                                                                     >  2:5429950/1‑93 (MQ=255)
                      aGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGgcgc                                                                   >  1:1288927/1‑93 (MQ=255)
                       gCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCa                                                                  >  2:2034872/1‑93 (MQ=255)
                         ggCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGc                                                                >  1:1586336/1‑93 (MQ=255)
                            aCGTTCACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGGAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCa                                                             <  1:6696087/93‑1 (MQ=255)
                            aCGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCa                                                             <  1:112877/93‑1 (MQ=255)
                            aCGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCa                                                             >  2:7830179/1‑93 (MQ=255)
                            aCGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCa                                                             <  2:6193828/93‑1 (MQ=255)
                                 cacCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTAc                                                        >  2:1415940/1‑93 (MQ=255)
                                   cctgcctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCg                                                      >  1:3623608/1‑93 (MQ=255)
                                    ctgcctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGa                                                     >  2:1694274/1‑93 (MQ=255)
                                     tgcctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAg                                                    >  2:3477314/1‑93 (MQ=255)
                                      acctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAgg                                                   <  1:765536/92‑1 (MQ=255)
                                        ctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGGGGAt                                                 >  2:5289208/1‑93 (MQ=255)
                                          ggCCTACCAACATCATCCAGGAGCCCCGAGAAGTCAGGCACGCAAAGGTGGAGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATac                                               >  2:5191610/1‑92 (MQ=255)
                                          ggCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAg                                               <  1:3262538/93‑1 (MQ=255)
                                          ggCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAg                                               >  1:5706242/1‑93 (MQ=255)
                                          ggCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAg                                               <  1:7586642/93‑1 (MQ=255)
                                           gCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAACGGTGGTGAATGCAGCAACCAGCGGCGGCGCCCAGTACCCGAGGACAGt                                              >  1:614669/1‑93 (MQ=255)
                                           gCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGt                                              <  1:2642382/93‑1 (MQ=255)
                                             cTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGtta                                            >  1:6516353/1‑92 (MQ=255)
                                             cTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTg                                            >  1:5114743/1‑93 (MQ=255)
                                             cTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTg                                            <  2:1627085/93‑1 (MQ=255)
                                              tACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGa                                           >  2:4242956/1‑93 (MQ=255)
                                              tACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGa                                           >  1:6812214/1‑93 (MQ=255)
                                              tACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGa                                           >  2:3632195/1‑93 (MQ=255)
                                              tACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGa                                           <  1:6338669/93‑1 (MQ=255)
                                              tACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGa                                           >  1:2820162/1‑93 (MQ=255)
                                               aCCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGaa                                          <  1:2013239/93‑1 (MQ=255)
                                               aCCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGaa                                          >  2:4512735/1‑93 (MQ=255)
                                               aCCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGaa                                          <  2:6138542/93‑1 (MQ=255)
                                                ccAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAg                                         >  2:836189/1‑93 (MQ=255)
                                                    catcatCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACt                                     <  1:7453700/93‑1 (MQ=255)
                                                      tcatcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGa                                   <  1:3549601/93‑1 (MQ=255)
                                                      tcatcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGa                                   >  2:7617630/1‑93 (MQ=255)
                                                       catcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGAt                                  <  2:632313/93‑1 (MQ=255)
                                                        atcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGAtt                                 <  1:5414130/93‑1 (MQ=255)
                                                         tcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTg                                >  1:3898522/1‑93 (MQ=255)
                                                            aGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCaa                             <  1:204815/93‑1 (MQ=255)
                                                             ggAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCaaa                            <  1:5992324/93‑1 (MQ=255)
                                                             ggAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCaaa                            <  2:7498708/93‑1 (MQ=255)
                                                             ggAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATACTTGAAGCACTGATTGCaac                            >  1:3803758/1‑92 (MQ=255)
                                                              gAGCCCAGAGAAGTCAGGCACGCAAAGGTGTTGACTGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAAc                           <  1:80654/93‑1 (MQ=255)
                                                                 cccAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGa                        <  2:3338512/93‑1 (MQ=255)
                                                                 cccAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGa                        >  2:7330966/1‑93 (MQ=255)
                                                                 cccAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGa                        >  2:2994999/1‑93 (MQ=255)
                                                                 cccAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGa                        >  1:827273/1‑93 (MQ=255)
                                                                 cccAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGa                        >  1:3355977/1‑93 (MQ=255)
                                                                  ccAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAg                       <  1:181880/93‑1 (MQ=255)
                                                                   cAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGc                      >  2:7664235/1‑93 (MQ=255)
                                                                   cAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGc                      <  1:3131225/93‑1 (MQ=255)
                                                                   cAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGc                      >  2:7122176/1‑93 (MQ=255)
                                                                   cAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGc                      >  2:2713345/1‑93 (MQ=255)
                                                                   cAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGc                      <  2:2159495/93‑1 (MQ=255)
                                                                    agagAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCa                     >  2:6498590/1‑93 (MQ=255)
                                                                    agagAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCa                     <  2:6912819/93‑1 (MQ=255)
                                                                       gaAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCg                  >  2:2415916/1‑93 (MQ=255)
                                                                         aGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGa                >  1:7900575/1‑93 (MQ=255)
                                                                         aGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGa                <  1:3450097/93‑1 (MQ=255)
                                                                          gTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGaa               <  1:5554445/93‑1 (MQ=255)
                                                                          gTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGaa               >  2:7198275/1‑93 (MQ=255)
                                                                          gTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGaa               >  1:7736185/1‑93 (MQ=255)
                                                                           tCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAg              >  2:3881217/1‑93 (MQ=255)
                                                                            cAGTCACGCAAAGGTGGTGAATGCAGAAACCAGCTGCGCAACCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAgc             >  2:4917707/1‑93 (MQ=255)
                                                                            cAGGCACGCAAAGGTGGTGAATGCAGAAGCCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCGCCGGAAgc             <  1:1471653/93‑1 (MQ=255)
                                                                             gggCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAgcc            >  1:1616310/2‑93 (MQ=255)
                                                                             aGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAgcc            >  1:4145398/1‑93 (MQ=255)
                                                                              ggCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAgccg           >  2:5748016/1‑93 (MQ=255)
                                                                                 aCGCAAAGGTGGTGAATGCAGAAGCCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCa        <  2:2913055/93‑1 (MQ=255)
                                                                                 aCGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCa        >  1:2324450/1‑93 (MQ=255)
                                                                                     aaaGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTAc    >  1:5956535/1‑93 (MQ=255)
                                                                                     aaaGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTAc    <  1:7216574/93‑1 (MQ=255)
                                                                                      aaGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACa   <  1:479237/93‑1 (MQ=255)
                                                                                       aGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTAccg  >  2:128758/1‑91 (MQ=255)
                                                                                       aGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACag  >  2:17628/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAG  >  NC_000913/4356685‑4356864

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