Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,485,872 C→T R27H (CGC→CAC)  pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,485,8720CT100.0% 534.4 / NA 137R27H (CGC→CAC) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (75/62);  total (75/62)

Only 100 of 137 total aligned reads displayed.

AAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAAC  >  NC_000913/4485727‑4486017
                                                                                                                                                 |                                                                                                                                                 
taGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGt                                                                                                                                              >  4:250138/2‑151 (MQ=255)
aaGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGt                                                                                                                                              >  2:1300267/1‑151 (MQ=255)
    agaATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGtt                                                                                                                                          <  2:905205/151‑1 (MQ=255)
      aaTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCg                                                                                                                                        <  2:944345/151‑1 (MQ=255)
         tCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGa                                                                                                                                       >  3:173480/1‑149 (MQ=255)
           gCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGAc                                                                                                                                   >  4:49419/1‑151 (MQ=255)
            cTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACg                                                                                                                                  >  1:491353/1‑151 (MQ=255)
              cGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACg                                                                                                                                  <  1:393674/149‑1 (MQ=255)
              cGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGc                                                                                                                                 <  2:806000/150‑1 (MQ=255)
               ggAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGccc                                                                                                                               >  1:771338/1‑151 (MQ=255)
                gAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCa                                                                                                                              <  2:1108285/151‑1 (MQ=255)
                  aaGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCacg                                                                                                                            >  2:938423/1‑151 (MQ=255)
                   aGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCacg                                                                                                                            >  1:210822/1‑150 (MQ=255)
                      aCATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCacgacg                                                                                                                         >  2:487089/1‑150 (MQ=255)
                       cATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGAt                                                                                                                       <  2:868788/151‑1 (MQ=255)
                          tCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGc                                                                                                                     <  2:1715720/150‑1 (MQ=255)
                           cGGAACATGGTGCAGCAACAATGTCTGTCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGc                                                                                                                     <  1:1213710/149‑1 (MQ=255)
                           cGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCa                                                                                                                    <  3:5212/150‑1 (MQ=255)
                            ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAgg                                                                                                                  >  3:2944/1‑151 (MQ=255)
                             gAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                                                                                                 >  1:1413681/1‑151 (MQ=255)
                               aCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCAc                                                                                                               <  3:368217/151‑1 (MQ=255)
                               aCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCAc                                                                                                               <  4:475637/151‑1 (MQ=255)
                               aCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCAc                                                                                                               >  2:341861/1‑151 (MQ=255)
                                 aTGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTc                                                                                                             >  3:583249/1‑151 (MQ=255)
                                   ggTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCg                                                                                                           <  2:210854/151‑1 (MQ=255)
                                   ggTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCg                                                                                                           >  2:1345962/1‑151 (MQ=255)
                                      gcagcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTg                                                                                                        >  3:296055/1‑151 (MQ=255)
                                       cagcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGt                                                                                                       >  1:22338/1‑151 (MQ=255)
                                          caacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtc                                                                                                    >  2:768155/1‑151 (MQ=255)
                                          caacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtc                                                                                                    <  3:127597/151‑1 (MQ=255)
                                           aacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtct                                                                                                   >  2:1127129/1‑151 (MQ=255)
                                            acaaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTcc                                                                                                 <  1:1671725/151‑1 (MQ=255)
                                              aaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCg                                                                                                >  2:1727413/1‑151 (MQ=255)
                                                    tGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTAc                                                                                          >  1:100109/1‑151 (MQ=255)
                                                     gCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTAcc                                                                                         <  2:167094/151‑1 (MQ=255)
                                                        ccGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCt                                                                                      >  2:218397/1‑151 (MQ=255)
                                                           gTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCtttt                                                                                   <  4:149069/151‑1 (MQ=255)
                                                            ttttCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttttt                                                                                  >  4:38756/1‑151 (MQ=255)
                                                             tttCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGAAGGCACTCCGCTGTTTCTCCGGGCTACCGCttttt                                                                                  <  1:428712/150‑1 (MQ=255)
                                                                ccttccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac                                                                              >  1:44838/1‑151 (MQ=255)
                                                                  ttccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                                                                             >  1:411166/1‑150 (MQ=255)
                                                                   tccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCggg                                                                                              >  4:548099/1‑132 (MQ=255)
                                                                   tccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa                                                                            >  1:669643/1‑150 (MQ=255)
                                                                   tccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa                                                                            >  2:1388431/1‑150 (MQ=255)
                                                                    ccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa                                                                           <  2:493924/150‑1 (MQ=255)
                                                                     cAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc                                                                         <  1:167130/151‑1 (MQ=255)
                                                                       gTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc                                                                       <  1:1503821/151‑1 (MQ=255)
                                                                         tCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc                                                                      >  1:624560/1‑150 (MQ=255)
                                                                         tCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc                                                                      <  1:877105/150‑1 (MQ=255)
                                                                          cGCCCCGACGGAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGGTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGTGCTACCGCTTTTTACACTAAACTCCAt                                                                    >  1:927392/1‑151 (MQ=255)
                                                                            ccccGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGc                                                                  >  1:776051/1‑151 (MQ=255)
                                                                              ccGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCAc                                                                >  4:66255/1‑151 (MQ=255)
                                                                               cGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACt                                                               <  3:520184/151‑1 (MQ=255)
                                                                                gACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTa                                                              <  1:1083532/151‑1 (MQ=255)
                                                                                 aCGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTAc                                                             >  2:867107/1‑151 (MQ=255)
                                                                                 aCGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTAc                                                             <  2:154351/151‑1 (MQ=255)
                                                                                    tAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCt                                                          >  1:602197/1‑151 (MQ=255)
                                                                                       cAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCt                                                       <  4:204461/151‑1 (MQ=255)
                                                                                        aGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTg                                                      >  2:719011/1‑151 (MQ=255)
                                                                                         gggCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTg                                                      >  2:1463826/1‑150 (MQ=255)
                                                                                           gcgcTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAAt                                                   <  2:1528559/151‑1 (MQ=255)
                                                                                           gcgcTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAAt                                                   <  2:1324209/151‑1 (MQ=255)
                                                                                             gcTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCt                                                 <  2:1509029/151‑1 (MQ=255)
                                                                                             gcTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCt                                                 >  1:216610/1‑151 (MQ=255)
                                                                                               tGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGAAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCtt                                                >  4:510803/1‑150 (MQ=255)
                                                                                                gATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTa                                               >  2:158096/1‑150 (MQ=255)
                                                                                                 aTGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTaaa                                             <  1:863154/151‑1 (MQ=255)
                                                                                                   gTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTaa                                              >  1:706510/1‑148 (MQ=255)
                                                                                                    tACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGAc                                          >  3:471971/1‑151 (MQ=255)
                                                                                                     aCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACa                                         <  2:1003055/151‑1 (MQ=255)
                                                                                                       ccATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAAc                                       <  1:1282317/151‑1 (MQ=255)
                                                                                                        cATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAAcg                                      <  2:1269907/151‑1 (MQ=255)
                                                                                                        cATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAAcg                                      <  1:64196/151‑1 (MQ=255)
                                                                                                           cGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAAcggcg                                   >  2:808610/1‑151 (MQ=255)
                                                                                                            gCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAAcggcgg                                  >  1:1185446/1‑151 (MQ=255)
                                                                                                             cTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAAcggcgg                                  >  1:595672/1‑150 (MQ=255)
                                                                                                             cTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGGTCGAAGACACCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGATCCTGAATCTTAAAGACAAcggc                                    >  2:761958/1‑148 (MQ=255)
                                                                                                              tGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAAcggcgg                                  >  1:458816/1‑149 (MQ=255)
                                                                                                               gATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGt                                 <  1:1686686/149‑1 (MQ=255)
                                                                                                               gATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGt                                 <  1:938305/149‑1 (MQ=255)
                                                                                                               gATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTgg                               <  1:1451682/151‑1 (MQ=255)
                                                                                                                aTTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAAcggcgg                                  >  1:1030658/1‑147 (MQ=255)
                                                                                                                aTTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTgg                               >  1:65246/1‑150 (MQ=255)
                                                                                                                     aTCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAg                         <  1:1728560/151‑1 (MQ=255)
                                                                                                                      tCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTAc                           >  2:23892/1‑148 (MQ=255)
                                                                                                                       cGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGa                        <  3:508881/150‑1 (MQ=255)
                                                                                                                       cGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGAt                       <  2:1722338/151‑1 (MQ=255)
                                                                                                                        gAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGa                        >  3:216259/1‑149 (MQ=255)
                                                                                                                        gAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATa                      <  2:1689537/151‑1 (MQ=255)
                                                                                                                             gTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAAtt                 >  1:469873/1‑151 (MQ=255)
                                                                                                                             gTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAAt                  <  3:457361/150‑1 (MQ=255)
                                                                                                                              ttCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAAtt                 >  2:867369/1‑150 (MQ=255)
                                                                                                                                cTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGc               >  1:184562/1‑150 (MQ=255)
                                                                                                                                  gCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCaa             <  2:491430/150‑1 (MQ=255)
                                                                                                                                    aaTCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCt          <  1:75216/151‑1 (MQ=255)
                                                                                                                                       cGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTc       <  2:1061783/151‑1 (MQ=255)
                                                                                                                                       cGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTc       <  2:33547/151‑1 (MQ=255)
                                                                                                                                            aaaGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAAATGCAAACTTTTGTaa   >  1:612351/1‑150 (MQ=255)
                                                                                                                                             aaGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAAc  >  1:491456/1‑150 (MQ=255)
                                                                                                                                             aaGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAAc  >  2:1264473/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAAC  >  NC_000913/4485727‑4486017

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