Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 446,850 A→G L27P (CTG→CCG)  cyoE ← protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913446,8500AG99.0% 997.0 / ‑1.4 311L27P (CTG→CCG) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (160/148);  total (163/148)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

Only 100 of 325 total aligned reads displayed.

TTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCT  >  NC_000913/446760‑446939
                                                                                          |                                                                                         
 tAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGa                                                                                        <  1:1216858/93‑1 (MQ=255)
 tAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGa                                                                                        <  1:463587/93‑1 (MQ=255)
   aaCACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCACGCGGGaac                                                                                      >  1:4860226/1‑92 (MQ=255)
       cacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATcccc                                                                                  >  2:8778901/1‑93 (MQ=255)
          aaCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCaa                                                                               <  2:1187048/93‑1 (MQ=255)
          aaCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCaa                                                                               <  1:8783525/93‑1 (MQ=255)
           aCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAAt                                                                              >  2:6015363/1‑93 (MQ=255)
            cccGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATc                                                                             >  2:7679095/1‑93 (MQ=255)
              cGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAc                                                                           >  1:857698/1‑93 (MQ=255)
              cGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAc                                                                           >  1:9780589/1‑93 (MQ=255)
               gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAcc                                                                          >  2:3190530/1‑93 (MQ=255)
                aCGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCg                                                                         >  1:10022631/1‑93 (MQ=255)
                aCGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCg                                                                         >  1:8456095/1‑93 (MQ=255)
                aCGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCg                                                                         >  1:8015764/1‑93 (MQ=255)
                aCGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCg                                                                         >  2:6798835/1‑93 (MQ=255)
                 cGCCACAACCAGAGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTCTGAGGCCAGCGGGAATCCCCCAATCACCga                                                                        >  2:9455468/1‑93 (MQ=255)
                   ccACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCgaga                                                                      <  1:3537355/93‑1 (MQ=255)
                   ccACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCgaga                                                                      >  1:6603237/1‑93 (MQ=255)
                    cacaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGAt                                                                     <  2:4534557/93‑1 (MQ=255)
                     acaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    >  1:6178763/1‑93 (MQ=255)
                     acaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    <  1:5242697/93‑1 (MQ=255)
                       aaCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAg                                                                  <  2:3835358/93‑1 (MQ=255)
                        aCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAgg                                                                 >  2:4861487/1‑93 (MQ=255)
                         ccAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGt                                                                <  2:6148843/93‑1 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGc                                                             >  2:6607418/1‑93 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGc                                                             <  2:5650421/93‑1 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGc                                                             >  1:2991430/1‑93 (MQ=255)
                             tGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGcc                                                            >  2:5821768/1‑93 (MQ=255)
                              gACACCCCAACCAGCGTGTAGATAAACAGCGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCa                                                           >  1:751604/1‑93 (MQ=255)
                                cacCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCaaa                                                         >  1:6489093/1‑93 (MQ=255)
                                 cccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAg                                                        <  2:1859332/92‑1 (MQ=255)
                                 acCCCAACGAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAg                                                        <  1:766506/93‑1 (MQ=255)
                                 acCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAg                                                        >  1:5698784/1‑93 (MQ=255)
                                 acCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAg                                                        <  1:7846341/93‑1 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                       <  1:1649555/93‑1 (MQ=255)
                                   cccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgat                                                      <  1:4928891/93‑1 (MQ=255)
                                   cccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgat                                                      <  2:8511613/93‑1 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    <  1:3220028/93‑1 (MQ=255)
                                      aaCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatgat                                                   <  1:3706411/93‑1 (MQ=255)
                                       aCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGatgatg                                                  <  2:618726/93‑1 (MQ=255)
                                        ccAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGc                                                 >  1:2305857/1‑93 (MQ=255)
                                        ccAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGc                                                 >  1:1038307/1‑93 (MQ=255)
                                        ccAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGc                                                 <  2:7510876/93‑1 (MQ=255)
                                         cAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                <  1:5460135/93‑1 (MQ=255)
                                         cAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                <  1:9656975/93‑1 (MQ=255)
                                         cAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                <  2:6097443/93‑1 (MQ=255)
                                           gCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTg                                              >  2:6393728/1‑93 (MQ=255)
                                           gCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTg                                              <  2:8910886/93‑1 (MQ=255)
                                           gCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGGTGCCAAAGATGATGCCTg                                              >  1:2696834/1‑93 (MQ=255)
                                             gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGt                                            >  2:5053610/1‑93 (MQ=255)
                                             gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGt                                            <  1:9977200/93‑1 (MQ=255)
                                              tgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                           >  1:2860923/1‑93 (MQ=255)
                                              tgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                           >  2:5310038/1‑93 (MQ=255)
                                                tAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                         >  1:6415351/1‑93 (MQ=255)
                                                tAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                         >  2:2791227/1‑93 (MQ=255)
                                                tAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                         >  2:7076321/1‑93 (MQ=255)
                                                tAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                         >  2:1657753/1‑93 (MQ=255)
                                                tAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                         >  1:9130449/1‑93 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACTGAGATCAGGTTGCAAAAGATGATGCCTGGTTTCGt                                       >  1:8294781/1‑93 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                       <  1:8116159/93‑1 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                       <  1:4229133/93‑1 (MQ=255)
                                                    tAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa                                     >  2:8262368/1‑93 (MQ=255)
                                                    tAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa                                     >  2:3581771/1‑93 (MQ=255)
                                                       aCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACtt                                  >  1:6271612/1‑93 (MQ=255)
                                                       aCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACtt                                  <  2:3937127/93‑1 (MQ=255)
                                                         aGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                <  2:4312698/93‑1 (MQ=255)
                                                         aGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                >  1:9478752/1‑93 (MQ=255)
                                                         aGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                >  1:1848711/1‑93 (MQ=255)
                                                         aGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                <  2:3307165/93‑1 (MQ=255)
                                                         aGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                <  2:194619/93‑1 (MQ=255)
                                                          ggggAGATTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                               <  2:1468718/93‑1 (MQ=255)
                                                           gggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAg                              >  2:1719120/1‑93 (MQ=255)
                                                             gATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                            <  2:7343990/93‑1 (MQ=255)
                                                               tAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAt                          >  2:5970148/1‑93 (MQ=255)
                                                                aaTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    >  2:2539811/1‑66 (MQ=255)
                                                                aaTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    <  1:2539811/66‑1 (MQ=255)
                                                                 aTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTg                        >  1:9042956/1‑93 (MQ=255)
                                                                 aTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTg                        <  2:7878590/93‑1 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                       >  2:6239099/1‑93 (MQ=255)
                                                                     aTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                    <  1:4591303/93‑1 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                   >  2:6526526/1‑93 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                   >  1:4879348/1‑93 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                   >  1:3845820/1‑93 (MQ=255)
                                                                       gctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaaa                  <  1:7297823/93‑1 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca                >  1:6857750/1‑93 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca                <  2:9210562/93‑1 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca                >  2:9143625/1‑93 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat               >  1:5095372/1‑93 (MQ=255)
                                                                           cccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatc              <  1:1351873/93‑1 (MQ=255)
                                                                            ccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatca             >  1:891026/1‑93 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat            >  1:8875856/1‑93 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat            >  2:4191333/1‑93 (MQ=255)
                                                                              tttGGGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           <  2:6368616/93‑1 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           <  2:3507054/93‑1 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           <  2:9787632/93‑1 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           >  1:3959336/1‑93 (MQ=255)
                                                                                tGAGGCCAGCGTGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc         <  1:2257631/93‑1 (MQ=255)
                                                                                tGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc         <  1:4756942/93‑1 (MQ=255)
                                                                                   ggCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgc      <  1:5232929/93‑1 (MQ=255)
                                                                                       agcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGct  >  1:6941985/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCT  >  NC_000913/446760‑446939

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