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.6% 877.0 / ‑5.0 276L27P (CTG→CCG) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (147/128);  total (147/129)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

Only 100 of 289 total aligned reads displayed.

TTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCT  >  NC_000913/446760‑446939
                                                                                          |                                                                                         
 tAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAAGAGGGGATAATCAATGCTGCCCTTTGAGGCCACCGGGa                                                                                        >  1:3971561/1‑93 (MQ=255)
     cacacaACCCGACGCCACGACCAGTGACATCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATcc                                                                                    <  1:4877738/93‑1 (MQ=255)
      acacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATccc                                                                                   >  2:5527172/1‑93 (MQ=255)
       cacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGGGTAGATAAACAGGGGATAATCAATACTGCCCTTTGAGGCTAGCGGGAATcccc                                                                                  >  2:5333130/1‑93 (MQ=255)
         caACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCa                                                                                <  1:6136189/93‑1 (MQ=255)
          aaCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTGGAGGCCAGCGGGAATCCCCCaa                                                                               <  1:1031593/93‑1 (MQ=255)
            cccGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATc                                                                             <  1:766752/93‑1 (MQ=255)
             ccGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCa                                                                            >  2:7815518/1‑93 (MQ=255)
              cGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAAAAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAc                                                                           <  1:230662/93‑1 (MQ=255)
                 cGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCga                                                                        <  2:7348311/93‑1 (MQ=255)
                 cGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCga                                                                        >  2:1309612/1‑93 (MQ=255)
                  gCCACAACCAGAGACACCCCAACCAGAGTGTAGATACACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGAGGGAATCCCACCATAACCgag                                                                       <  2:7001374/93‑1 (MQ=255)
                    cacaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGAt                                                                     >  2:7897963/1‑93 (MQ=255)
                    cacaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGAt                                                                     <  2:1645188/93‑1 (MQ=255)
                     ccaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    <  1:1508162/92‑1 (MQ=255)
                      caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGTCCAGCGGGAATCCCCCAATCACCGAGATCa                                                                   >  2:7565439/1‑93 (MQ=255)
                      caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCa                                                                   >  2:2764178/1‑93 (MQ=255)
                      caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCa                                                                   <  2:2910451/93‑1 (MQ=255)
                      caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCa                                                                   >  1:5943139/1‑93 (MQ=255)
                       aaCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAg                                                                  >  1:2271570/1‑93 (MQ=255)
                          cAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGtt                                                               <  1:6443548/93‑1 (MQ=255)
                           aGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGttc                                                              >  2:4952669/1‑92 (MQ=255)
                           aGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTg                                                              >  2:282335/1‑93 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGc                                                             >  1:2043806/1‑93 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGc                                                             <  2:7690331/93‑1 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                       <  2:123373/93‑1 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    <  2:3603478/93‑1 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    >  2:3026348/1‑93 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    <  2:6087046/93‑1 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    <  2:629326/93‑1 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    <  2:1428156/93‑1 (MQ=255)
                                      aaCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCCATCACCGAGATCAGGTTGCCAAAgaggat                                                   >  1:7298655/1‑93 (MQ=255)
                                      aaCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatgat                                                   <  1:6228109/93‑1 (MQ=255)
                                      aaCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAACCCACGAGATCAGGTTGCCAAAgatgct                                                   >  1:5273225/1‑91 (MQ=255)
                                       aCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGatgatg                                                  >  1:3788389/1‑93 (MQ=255)
                                       aCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGatgatg                                                  <  1:990173/93‑1 (MQ=255)
                                       aCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGatgatg                                                  <  1:186150/93‑1 (MQ=255)
                                        ccAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGc                                                 <  1:3244702/93‑1 (MQ=255)
                                        ccAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGc                                                 <  2:949401/93‑1 (MQ=255)
                                          aGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCt                                               <  2:4339511/93‑1 (MQ=255)
                                           gCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTg                                              <  2:4004816/93‑1 (MQ=255)
                                            cGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGACAAAGATGATGCCTgg                                             >  1:4724196/1‑93 (MQ=255)
                                             gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGt                                            >  1:7849886/1‑93 (MQ=255)
                                             gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGt                                            >  1:4253847/1‑93 (MQ=255)
                                              tgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                           >  2:3960637/1‑93 (MQ=255)
                                              tgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                           >  2:6462339/1‑93 (MQ=255)
                                              tgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                           <  1:2527633/93‑1 (MQ=255)
                                               gtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGttt                                          <  1:4205868/93‑1 (MQ=255)
                                               gtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGttt                                          >  1:1431737/1‑93 (MQ=255)
                                                tAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                         <  2:4308604/93‑1 (MQ=255)
                                                 aGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCg                                        >  2:697171/1‑93 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                       <  2:7213815/93‑1 (MQ=255)
                                                    tAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa                                     >  2:2309083/1‑93 (MQ=255)
                                                     aaaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCGCCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTAc                                    <  1:5838279/93‑1 (MQ=255)
                                                     aaaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTAc                                    <  2:1577655/93‑1 (MQ=255)
                                                      aaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACt                                   >  1:6843806/1‑93 (MQ=255)
                                                         aGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                >  1:514195/1‑93 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                               <  2:21292/93‑1 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                               >  1:5557803/1‑93 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                               >  1:2149467/1‑93 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                               <  1:3120869/93‑1 (MQ=255)
                                                           gggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAg                              <  1:1657560/93‑1 (MQ=255)
                                                             gATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                            >  2:3206735/1‑93 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTTGTTTCGTTACTTGCAGGTATTGc                       >  1:3065339/1‑93 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                       <  2:6035233/93‑1 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                       >  1:1546328/1‑93 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                       <  2:3791797/93‑1 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                       >  2:2616648/1‑93 (MQ=255)
                                                                    aaTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCtt                     >  2:5822388/1‑93 (MQ=255)
                                                                     aTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGTTTa                    >  1:5512876/1‑93 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                   <  1:3642729/93‑1 (MQ=255)
                                                                       gctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaaa                  >  2:3273000/1‑93 (MQ=255)
                                                                       gctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaaa                  <  1:7048412/93‑1 (MQ=255)
                                                                        ctgcCCTTTGAGGCCAGCGGGAATCCCCCAATTAACGAGCTCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc                 >  2:630182/1‑93 (MQ=255)
                                                                        ctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTGTTGCTTAAAc                 <  1:2227356/93‑1 (MQ=255)
                                                                        ctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc                 <  1:2790020/93‑1 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca                >  1:2394897/1‑93 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcct               >  1:6075957/1‑91 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat               >  1:4641743/1‑93 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat               >  1:6710424/1‑93 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat               <  2:335361/93‑1 (MQ=255)
                                                                           cccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCCGGTAGTGCTTAACcatc              >  1:1248614/1‑93 (MQ=255)
                                                                           cccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatc              >  1:7085916/1‑93 (MQ=255)
                                                                           cccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatc              <  2:6585867/93‑1 (MQ=255)
                                                                            ccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatca             <  1:7454661/93‑1 (MQ=255)
                                                                            ccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatca             <  2:3297776/93‑1 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat            <  1:1300953/93‑1 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           >  2:594422/1‑93 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           >  2:5722272/1‑93 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           >  2:2174770/1‑93 (MQ=255)
                                                                               ttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATaa          <  1:3036122/93‑1 (MQ=255)
                                                                                tGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc         <  1:657989/93‑1 (MQ=255)
                                                                                tGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc         <  2:752055/93‑1 (MQ=255)
                                                                                  aGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccg       <  1:4482777/93‑1 (MQ=255)
                                                                                   ggCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgc      <  1:1190975/93‑1 (MQ=255)
                                                                                   ggCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgc      >  1:590649/1‑93 (MQ=255)
                                                                                   ggCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgc      >  2:5043179/1‑93 (MQ=255)
                                                                                     ccAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgccg    >  1:5384467/1‑93 (MQ=255)
                                                                                       agcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGct  >  2:6800379/1‑93 (MQ=255)
                                                                                       agcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGct  >  2:1835125/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCT  >  NC_000913/446760‑446939

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