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,8500AG100.0% 776.3 / NA 241L27P (CTG→CCG) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (121/120);  total (121/120)

Only 100 of 251 total aligned reads displayed.

TTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCT  >  NC_000913/446760‑446939
                                                                                          |                                                                                         
ttAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCggg                                                                                         >  2:234097/1‑93 (MQ=255)
  aaaCACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGaa                                                                                       >  1:3253311/1‑93 (MQ=255)
   aaCACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAAt                                                                                      >  2:5725401/1‑93 (MQ=255)
   aaCACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAAt                                                                                      <  2:4799660/93‑1 (MQ=255)
    acacacAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATc                                                                                     >  1:4126235/1‑93 (MQ=255)
     cacacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGCTAATCAATGCTGCCCTTTGAGGCCAGCGGGAATcc                                                                                    >  2:3081789/1‑93 (MQ=255)
       cacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATcccc                                                                                  <  1:3995395/93‑1 (MQ=255)
       cacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATcccc                                                                                  >  2:6298654/1‑93 (MQ=255)
       cacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATcccc                                                                                  <  1:6017162/93‑1 (MQ=255)
         caACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCa                                                                                >  1:3191774/1‑93 (MQ=255)
           aCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAAt                                                                              >  2:788539/1‑93 (MQ=255)
            cccGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATc                                                                             >  1:2467079/1‑93 (MQ=255)
             ccGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCa                                                                            >  1:3531002/1‑93 (MQ=255)
              cGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAc                                                                           <  1:4564228/93‑1 (MQ=255)
               gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAcc                                                                          >  1:2441145/1‑93 (MQ=255)
               gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAcc                                                                          <  1:2337913/93‑1 (MQ=255)
               gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAcc                                                                          >  2:5335657/1‑93 (MQ=255)
                aCGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCAAATCACCg                                                                         >  1:1949599/1‑93 (MQ=255)
                 cGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCga                                                                        <  1:2216849/93‑1 (MQ=255)
                   ccACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCgaga                                                                      >  1:500960/1‑93 (MQ=255)
                   ccACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCgaga                                                                      <  2:751981/93‑1 (MQ=255)
                     acaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    <  2:5342661/93‑1 (MQ=255)
                     acaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    <  1:2729180/93‑1 (MQ=255)
                     acaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    <  1:5198476/93‑1 (MQ=255)
                     acaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    >  1:1161863/1‑93 (MQ=255)
                      caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATACCCCAATCACCGAGATCa                                                                   >  1:3350274/1‑93 (MQ=255)
                        aCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAgg                                                                 >  2:2175918/1‑93 (MQ=255)
                         ccAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGt                                                                >  1:5569893/1‑93 (MQ=255)
                         ccAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGt                                                                <  1:5481037/93‑1 (MQ=255)
                         ccAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGt                                                                <  1:5408732/93‑1 (MQ=255)
                          cAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGtt                                                               <  2:1543344/93‑1 (MQ=255)
                           aGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTg                                                              >  2:6054742/1‑93 (MQ=255)
                           aGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTg                                                              >  1:3635388/1‑93 (MQ=255)
                           aGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTg                                                              <  2:970282/93‑1 (MQ=255)
                           aGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTg                                                              <  1:518822/93‑1 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGc                                                             >  2:5514751/1‑93 (MQ=255)
                              gACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCa                                                           >  1:209632/1‑93 (MQ=255)
                               acacCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCaa                                                          >  1:1084080/1‑93 (MQ=255)
                                 acCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAg                                                        <  2:5856600/93‑1 (MQ=255)
                                 acCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAg                                                        >  1:825792/1‑93 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                       <  1:1504528/93‑1 (MQ=255)
                                    ccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatg                                                     <  1:5656475/93‑1 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    >  1:182772/1‑93 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    <  1:430935/93‑1 (MQ=255)
                                      agccAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatgat                                                   <  2:1000848/91‑1 (MQ=255)
                                       aCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGatgatg                                                  <  1:4219177/93‑1 (MQ=255)
                                        ccAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTAGCAAACATGATGc                                                 >  1:2675016/1‑93 (MQ=255)
                                         cAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                <  1:6163730/93‑1 (MQ=255)
                                         cAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                <  2:5200240/93‑1 (MQ=255)
                                            cGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCCAAGATGATGCCTgg                                             >  1:5326414/1‑93 (MQ=255)
                                              tgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                           >  1:1036258/1‑93 (MQ=255)
                                               gtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGCTCAGGTTGCCAAAGATGATGCCTGGttt                                          >  1:393850/1‑93 (MQ=255)
                                                 aGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCg                                        >  2:1826649/1‑93 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                       <  1:1893142/93‑1 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                       <  2:3745647/93‑1 (MQ=255)
                                                   aTAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGtt                                      >  2:214797/1‑93 (MQ=255)
                                                    tAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa                                     <  2:4557532/93‑1 (MQ=255)
                                                     aaaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTAc                                    <  2:3576244/93‑1 (MQ=255)
                                                      aaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACt                                   >  2:264316/1‑93 (MQ=255)
                                                       aCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACtt                                  >  1:3423055/1‑93 (MQ=255)
                                                       aCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACtt                                  <  1:2704979/93‑1 (MQ=255)
                                                         aGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                >  2:3598452/1‑93 (MQ=255)
                                                           gggATAAACAATGCTGCCCGTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAg                              <  1:5857293/93‑1 (MQ=255)
                                                            ggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAgg                             <  1:1276578/93‑1 (MQ=255)
                                                            ggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAgg                             <  2:4408358/93‑1 (MQ=255)
                                                             gATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                            <  2:3052614/93‑1 (MQ=255)
                                                              aTAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                           <  1:1856286/93‑1 (MQ=255)
                                                                aaTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAtt                         <  1:4764673/93‑1 (MQ=255)
                                                                aaTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAtt                         <  2:1114618/93‑1 (MQ=255)
                                                                 aTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTg                        >  1:3020814/1‑93 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGTATCCCGCAATCACCGAGATCAGGGTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTCTTGc                       >  2:5968904/1‑93 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                       >  1:4982327/1‑93 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                       <  2:3387098/93‑1 (MQ=255)
                                                                   cAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCt                      <  2:5452554/93‑1 (MQ=255)
                                                                    aaTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCtt                     >  2:2524345/1‑93 (MQ=255)
                                                                    aaTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCtt                     >  2:1405488/1‑93 (MQ=255)
                                                                     aTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                    <  2:1269422/93‑1 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                   >  2:3304710/1‑93 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                   >  1:235396/1‑93 (MQ=255)
                                                                        ctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc                 >  2:2765591/1‑93 (MQ=255)
                                                                        ctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc                 <  1:1406877/93‑1 (MQ=255)
                                                                        ctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc                 <  1:2561034/93‑1 (MQ=255)
                                                                        ctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc                 >  1:2914752/1‑93 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca                >  2:4219814/1‑93 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca                <  2:3971183/93‑1 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat               >  1:6144778/1‑93 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat               >  1:2807860/1‑93 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat               >  2:1441023/1‑93 (MQ=255)
                                                                           cccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTCCTTGCAGATATTGCTTAAAcatc              >  1:407749/1‑93 (MQ=255)
                                                                            ccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatca             >  1:5149769/1‑93 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat            <  1:5521393/93‑1 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat            <  1:4739567/93‑1 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat            <  2:4319474/93‑1 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           <  2:2670263/93‑1 (MQ=255)
                                                                               ttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATaa          <  1:2314468/93‑1 (MQ=255)
                                                                                tGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc         <  2:6155717/93‑1 (MQ=255)
                                                                                tGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc         >  2:2861666/1‑93 (MQ=255)
                                                                                  aGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccg       <  1:45406/93‑1 (MQ=255)
                                                                                   ggCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgc      >  2:4484798/1‑93 (MQ=255)
                                                                                       agcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGct  >  2:5032027/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: