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.4% 520.5 / ‑4.8 167L27P (CTG→CCG) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (87/79);  total (88/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 172 total aligned reads displayed.

TTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGC  >  NC_000913/446760‑446938
                                                                                          |                                                                                        
 tAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGa                                                                                       >  2:2112111/1‑93 (MQ=255)
   aaCACACAACCCGACGCCACAACCAGTGACACCCCACCCAGCGTGTAGATAAACAGGGAATAATCAATGCTGCCCTTTGAGGCCAGCGGGAAt                                                                                     <  1:3365290/93‑1 (MQ=255)
     cacacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATcc                                                                                   >  1:3894995/1‑93 (MQ=255)
      acacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATccc                                                                                  >  2:3590541/1‑93 (MQ=255)
      acacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATccc                                                                                  >  1:4484860/1‑93 (MQ=255)
        acaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATccccc                                                                                <  1:4110242/93‑1 (MQ=255)
         caACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCa                                                                               <  2:3225128/93‑1 (MQ=255)
          aaCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCaa                                                                              <  1:3399571/93‑1 (MQ=255)
           aCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAAt                                                                             <  2:2313866/93‑1 (MQ=255)
           aCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAAt                                                                             >  2:3875973/1‑93 (MQ=255)
             ccGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCa                                                                           >  2:4636914/1‑93 (MQ=255)
              cGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAc                                                                          <  2:2816986/93‑1 (MQ=255)
              cGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAc                                                                          <  1:4407561/93‑1 (MQ=255)
              cGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAc                                                                          <  1:2358938/93‑1 (MQ=255)
                 cGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCga                                                                       <  2:4639433/93‑1 (MQ=255)
                 cGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCga                                                                       >  1:1617195/1‑93 (MQ=255)
                  gCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCgag                                                                      >  2:1380731/1‑93 (MQ=255)
                  gCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCACCAATCACCgag                                                                      >  2:1617143/1‑93 (MQ=255)
                   ccACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCgaga                                                                     >  2:592994/1‑93 (MQ=255)
                       aaCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAg                                                                 <  2:4359709/93‑1 (MQ=255)
                       aaCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCAAATCACCGAGATCAg                                                                 >  1:3075827/1‑93 (MQ=255)
                         ccAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGt                                                               <  2:4156128/93‑1 (MQ=255)
                          cAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGtt                                                              >  1:2570119/1‑93 (MQ=255)
                          cAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGtt                                                              <  1:977365/93‑1 (MQ=255)
                           aGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTg                                                             >  2:828623/1‑93 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGc                                                            >  2:1054574/1‑93 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGc                                                            >  2:4319300/1‑93 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGc                                                            >  1:736870/1‑93 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGTGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                      <  2:1624144/93‑1 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                      <  1:285104/93‑1 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                      <  2:1802960/93‑1 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                      <  1:24783/93‑1 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                      <  2:3397728/93‑1 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                   <  1:96635/93‑1 (MQ=255)
                                      aaCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatgat                                                  <  2:240533/93‑1 (MQ=255)
                                      aaCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatgat                                                  <  1:3695996/93‑1 (MQ=255)
                                      aaCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatgat                                                  >  2:2348986/1‑93 (MQ=255)
                                        ccAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGc                                                >  1:4291241/1‑93 (MQ=255)
                                          aGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCt                                              >  1:2127020/1‑93 (MQ=255)
                                           gCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTg                                             >  1:1030281/1‑93 (MQ=255)
                                            cGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTgg                                            >  2:3520068/1‑93 (MQ=255)
                                             gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGt                                           <  1:422982/93‑1 (MQ=255)
                                             gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGt                                           >  2:736111/1‑93 (MQ=255)
                                              tgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                          >  1:1146912/1‑93 (MQ=255)
                                                tAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                        >  2:2346135/1‑93 (MQ=255)
                                                 aGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCg                                       >  1:2979722/1‑93 (MQ=255)
                                                 aGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCg                                       >  2:1777363/1‑93 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                      <  2:3802445/93‑1 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                      <  2:2834965/93‑1 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                      <  2:654932/93‑1 (MQ=255)
                                                  gATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                      <  2:665985/93‑1 (MQ=255)
                                                   aTAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGtt                                     <  2:4192861/93‑1 (MQ=255)
                                                   aTAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGtt                                     >  2:1045671/1‑93 (MQ=255)
                                                    tAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa                                    >  1:41901/1‑93 (MQ=255)
                                                      aaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACt                                  >  2:3776463/1‑93 (MQ=255)
                                                       aCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACtt                                 >  1:3953681/1‑93 (MQ=255)
                                                       aCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACtt                                 <  1:3178707/93‑1 (MQ=255)
                                                        cAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTg                                <  1:658043/93‑1 (MQ=255)
                                                         aGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                               <  2:4404425/93‑1 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCGAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                              <  2:1370044/93‑1 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                              <  1:4614618/93‑1 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                              >  2:211288/1‑93 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                              >  1:1863741/1‑93 (MQ=255)
                                                           gggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAg                             <  2:1446002/93‑1 (MQ=255)
                                                            ggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAgg                            <  1:1382918/93‑1 (MQ=255)
                                                            ggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAgg                            >  1:143050/1‑93 (MQ=255)
                                                             gATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                           >  2:2655907/1‑93 (MQ=255)
                                                              aTAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                          <  2:2513055/93‑1 (MQ=255)
                                                              aTAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                          <  1:1532148/93‑1 (MQ=255)
                                                               tAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAt                         <  1:525695/93‑1 (MQ=255)
                                                                aaTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAtt                        >  2:1920426/1‑93 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                      >  1:3227834/1‑93 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAACAGGAATACCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                      <  1:3473900/93‑1 (MQ=255)
                                                                   cAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCt                     <  2:2010998/93‑1 (MQ=255)
                                                                   cAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCt                     >  2:3520402/1‑93 (MQ=255)
                                                                    aaTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCtt                    <  2:1173171/93‑1 (MQ=255)
                                                                    aaTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCtt                    >  1:4554506/1‑93 (MQ=255)
                                                                     aTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                   >  1:2907623/1‑93 (MQ=255)
                                                                     aTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                   >  1:1868870/1‑93 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                  <  1:3784894/93‑1 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                  >  2:1806033/1‑93 (MQ=255)
                                                                       gctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaaa                 <  1:2794354/93‑1 (MQ=255)
                                                                        ctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc                <  2:405515/93‑1 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca               >  2:3680440/1‑93 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca               <  1:2020121/93‑1 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat              <  1:1764380/93‑1 (MQ=255)
                                                                            ccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatca            >  1:1933195/1‑93 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat           >  1:1984944/1‑93 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat           >  1:1475469/1‑93 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa          <  2:158246/93‑1 (MQ=255)
                                                                               ttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATaa         <  2:3252196/93‑1 (MQ=255)
                                                                                 gAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAcc       >  1:1187521/1‑93 (MQ=255)
                                                                                 gAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAcc       >  2:2685117/1‑93 (MQ=255)
                                                                                   ggCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgc     <  1:2081837/93‑1 (MQ=255)
                                                                                   ggCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgc     <  1:3660938/93‑1 (MQ=255)
                                                                                   ggCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAcagc     <  2:494833/93‑1 (MQ=255)
                                                                                    gCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgcc    >  1:4008087/1‑93 (MQ=255)
                                                                                      cagcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACcgccgc  >  2:3979416/1‑93 (MQ=255)
                                                                                      cagcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACcgccgc  <  1:281561/93‑1 (MQ=255)
                                                                                      cagcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACcgccgc  <  1:3578685/93‑1 (MQ=255)
                                                                                          |                                                                                        
TTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGC  >  NC_000913/446760‑446938

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