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,8500AG98.8% 522.0 / ‑2.9 163L27P (CTG→CCG) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (78/82);  total (79/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

Only 100 of 172 total aligned reads displayed.

TTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCT  >  NC_000913/446760‑446939
                                                                                          |                                                                                         
       cacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATcccc                                                                                  <  1:4359834/93‑1 (MQ=255)
        acaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATccccc                                                                                 >  1:371297/1‑93 (MQ=255)
        acaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATccccc                                                                                 >  2:3249241/1‑93 (MQ=255)
          aaCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCaa                                                                               <  2:1069971/93‑1 (MQ=255)
           aCCCGACGCCACCACCAGTGACACCCCAACCAGCGTGTAGATACACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCCAt                                                                              <  2:4220630/93‑1 (MQ=255)
           aCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAAt                                                                              >  1:884614/1‑93 (MQ=255)
            ccccACGCCATAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCATTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATc                                                                             <  2:7594/93‑1 (MQ=255)
              cGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCAc                                                                           >  2:3726392/1‑93 (MQ=255)
                     ccaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    <  1:2603414/92‑1 (MQ=255)
                     acaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    <  1:4546584/93‑1 (MQ=255)
                     acaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    >  2:1158025/1‑93 (MQ=255)
                     acaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    >  1:1939861/1‑93 (MQ=255)
                     acaACCAGTGACACCCCAAACAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATc                                                                    >  1:4597926/1‑93 (MQ=255)
                      ccaCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCa                                                                   <  1:3865041/91‑1 (MQ=255)
                      caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCCCCGAGATCa                                                                   <  1:1857529/93‑1 (MQ=255)
                         ccAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGt                                                                >  2:2251606/1‑93 (MQ=255)
                           aGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTg                                                              >  2:1681666/1‑93 (MQ=255)
                           aGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTg                                                              >  2:3682650/1‑93 (MQ=255)
                            gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGc                                                             >  2:2141586/1‑93 (MQ=255)
                             tGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGcc                                                            <  1:4121637/93‑1 (MQ=255)
                              gACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCa                                                           >  1:303164/1‑93 (MQ=255)
                                 acCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCaaaa                                                        >  1:3221671/1‑92 (MQ=255)
                                 acCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAg                                                        >  1:1191939/1‑93 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                       <  1:353123/93‑1 (MQ=255)
                                  ccccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                                                       <  1:3226864/93‑1 (MQ=255)
                                   cccAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgat                                                      <  1:2547754/93‑1 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    <  2:531983/93‑1 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    >  2:1345146/1‑93 (MQ=255)
                                     cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatga                                                    >  2:644340/1‑93 (MQ=255)
                                      aaCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatgat                                                   <  1:513860/93‑1 (MQ=255)
                                      aaCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatgat                                                   >  2:76752/1‑93 (MQ=255)
                                       aCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGatgatg                                                  >  1:3591340/1‑93 (MQ=255)
                                        ccAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGc                                                 <  1:1111831/93‑1 (MQ=255)
                                         cAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                >  1:710792/1‑93 (MQ=255)
                                         cAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                >  1:1467266/1‑93 (MQ=255)
                                         cAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                >  2:961098/1‑93 (MQ=255)
                                         cAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                <  1:4345954/93‑1 (MQ=255)
                                          aGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCt                                               <  1:1982097/93‑1 (MQ=255)
                                           gCGTGTAGATAAACAGGGGATACTCAATGCTGCCCTTTGAGACCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTg                                              <  1:4212635/93‑1 (MQ=255)
                                           gCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTg                                              >  1:3698759/1‑93 (MQ=255)
                                            cGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTgg                                             <  1:4150466/93‑1 (MQ=255)
                                             gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGt                                            <  2:2654968/93‑1 (MQ=255)
                                              tgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                           >  1:3299843/1‑93 (MQ=255)
                                               gtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGttt                                          <  1:560633/93‑1 (MQ=255)
                                               gtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGttt                                          >  2:3264883/1‑93 (MQ=255)
                                                 aGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCg                                        <  1:4016034/93‑1 (MQ=255)
                                                   aTAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGtt                                      <  2:1215287/93‑1 (MQ=255)
                                                    tAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa                                     >  1:194003/1‑93 (MQ=255)
                                                     aaaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTAc                                    >  1:2270390/1‑93 (MQ=255)
                                                     aaaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTAc                                    <  2:1994072/93‑1 (MQ=255)
                                                      aaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACt                                   <  2:2770519/93‑1 (MQ=255)
                                                      aaCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACt                                   >  1:1159479/1‑93 (MQ=255)
                                                       aCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACtt                                  >  2:1451162/1‑93 (MQ=255)
                                                       aCAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACtt                                  >  1:3502120/1‑93 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                               >  1:3041348/1‑93 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                               >  1:3359033/1‑93 (MQ=255)
                                                          ggggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa                               >  1:4334011/1‑93 (MQ=255)
                                                           gggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAg                              <  2:3492347/93‑1 (MQ=255)
                                                            ggATAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAgg                             <  2:1871991/93‑1 (MQ=255)
                                                              atatTCAATGCGGTCGTTTGAGGCCAGCGGGATGCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAt                          <  1:2008929/93‑1 (MQ=255)
                                                              aTAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                           >  2:3687477/1‑93 (MQ=255)
                                                               tAATCAATGCTGCCCTTTGAGGCCAGCTGGAATCCCCCAATCCCCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAt                          <  2:3627595/93‑1 (MQ=255)
                                                               tAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAt                          >  1:1032549/1‑93 (MQ=255)
                                                               tAATCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAt                          >  1:273000/1‑93 (MQ=255)
                                                                aaTCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAtt                         >  1:2616276/1‑93 (MQ=255)
                                                                  tCAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                       >  2:524045/1‑93 (MQ=255)
                                                                   cAATGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCt                      <  2:166128/93‑1 (MQ=255)
                                                                    aaTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCtt                     >  2:2891479/1‑93 (MQ=255)
                                                                    aaTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCtt                     <  2:3650140/93‑1 (MQ=255)
                                                                    aaTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCtt                     >  1:106462/1‑93 (MQ=255)
                                                                    aaTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCtt                     >  1:2419124/1‑93 (MQ=255)
                                                                     aTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                    >  2:2102743/1‑93 (MQ=255)
                                                                     aTGCTGCCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                    >  2:447915/1‑93 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                   <  1:4653742/93‑1 (MQ=255)
                                                                      tgctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTAGTTACTTGCAGGTATTGCTTaa                   <  1:802698/93‑1 (MQ=255)
                                                                        ctgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc                 >  1:629147/1‑93 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca                <  1:2487558/93‑1 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca                >  1:2137278/1‑93 (MQ=255)
                                                                         tgcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAca                <  1:2857408/93‑1 (MQ=255)
                                                                          gcCCTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat               <  2:57557/93‑1 (MQ=255)
                                                                            ccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGCTGAGGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatca             >  2:1749331/1‑93 (MQ=255)
                                                                            ccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatca             >  2:978074/1‑93 (MQ=255)
                                                                            ccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatca             <  1:3340767/93‑1 (MQ=255)
                                                                            ccTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatca             <  1:4633855/93‑1 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat            <  2:4492201/93‑1 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat            <  1:2314344/93‑1 (MQ=255)
                                                                             cTTTGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat            >  2:4064753/1‑93 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           <  1:2572699/93‑1 (MQ=255)
                                                                              tttGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATa           >  2:3852232/1‑93 (MQ=255)
                                                                               ttGAGGCCAGCGCGGATCCCCCAATCCGCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTATTTGCAGGTATTGCTTAAACATCATaa          <  1:1576526/93‑1 (MQ=255)
                                                                                tGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc         <  1:827596/93‑1 (MQ=255)
                                                                                tGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc         <  1:950531/93‑1 (MQ=255)
                                                                                tGAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc         <  2:1787019/93‑1 (MQ=255)
                                                                                 gAGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAcc        >  2:2250155/1‑93 (MQ=255)
                                                                                  aGGCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccg       >  2:1008587/1‑93 (MQ=255)
                                                                                   ggCCAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgc      >  1:1759862/1‑93 (MQ=255)
                                                                                     ccAGCGGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgccg    <  2:3014703/93‑1 (MQ=255)
                                                                                      cagcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACcgccgc   >  2:4650101/1‑93 (MQ=255)
                                                                                      cagcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACcgccgc   <  2:2881253/93‑1 (MQ=255)
                                                                                       agcggGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGct  >  1:145518/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCT  >  NC_000913/446760‑446939

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