Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,333,152 A→G L16P (CTA→CCA)  yjdB ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,333,1520AG100.0% 851.1 / NA 263L16P (CTA→CCA) yjdBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (126/137);  total (126/137)

Only 100 of 288 total aligned reads displayed.

ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241
                                                                                          |                                                                                         
acCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggt                                                                                         <  1:3794996/93‑1 (MQ=255)
acCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggt                                                                                         >  2:3526787/1‑93 (MQ=255)
acCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggt                                                                                         <  2:3548712/93‑1 (MQ=255)
acCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggt                                                                                         >  1:2792689/1‑93 (MQ=255)
 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        >  2:1233848/1‑93 (MQ=255)
 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        >  1:5235479/1‑93 (MQ=255)
 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        <  1:4398859/93‑1 (MQ=255)
 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        >  1:205626/1‑93 (MQ=255)
  cTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggtag                                                                                       <  2:1628148/93‑1 (MQ=255)
  cTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggtag                                                                                       <  1:3567895/93‑1 (MQ=255)
   tGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAgg                                                                                      >  1:4826727/1‑93 (MQ=255)
      ttAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGgcgc                                                                                   >  2:1482215/1‑93 (MQ=255)
         aaaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTc                                                                                <  1:4341290/93‑1 (MQ=255)
         aaaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTc                                                                                <  1:3124868/93‑1 (MQ=255)
            aaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAAc                                                                             <  1:6257850/93‑1 (MQ=255)
            aaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAAc                                                                             <  2:6221241/93‑1 (MQ=255)
             aGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACa                                                                            <  2:4453166/93‑1 (MQ=255)
             aGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACa                                                                            <  2:6101641/93‑1 (MQ=255)
              ggCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACAt                                                                           >  2:4638996/1‑93 (MQ=255)
              ggCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACAt                                                                           >  2:2625873/1‑93 (MQ=255)
               gCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATc                                                                          >  2:344242/1‑93 (MQ=255)
                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      >  2:3939430/1‑93 (MQ=255)
                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      <  1:1968600/93‑1 (MQ=255)
                    atTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgt                                                                     <  2:475392/93‑1 (MQ=255)
                    atTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgt                                                                     >  2:5348196/1‑93 (MQ=255)
                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    >  1:769307/1‑93 (MQ=255)
                      tCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtgt                                                                   <  1:236855/93‑1 (MQ=255)
                      tCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtgt                                                                   <  2:4081336/93‑1 (MQ=255)
                      tCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtgt                                                                   <  1:4278733/93‑1 (MQ=255)
                          gCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTcc                                                               <  2:2640439/93‑1 (MQ=255)
                          gCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTcc                                                               <  2:5576665/93‑1 (MQ=255)
                           cAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCa                                                              >  2:2030173/1‑93 (MQ=255)
                            aGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCAt                                                             >  1:5644017/1‑93 (MQ=255)
                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            <  2:5260950/93‑1 (MQ=255)
                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            <  2:5899239/93‑1 (MQ=255)
                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGAGCTTCAACATCACGGTGATTCCATc                                                            >  1:951158/1‑93 (MQ=255)
                               tAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGa                                                          >  1:1477379/1‑93 (MQ=255)
                                agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGTGCTTCAACATCACGGTGTTTCCATCGaa                                                         <  1:3878300/93‑1 (MQ=255)
                                agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGaa                                                         <  2:489132/93‑1 (MQ=255)
                                 gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAAc                                                        >  1:5877309/1‑93 (MQ=255)
                                 gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAAc                                                        <  2:4350685/93‑1 (MQ=255)
                                   gaTATAAAAAGCGGCCCACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGTGGTTTCCATCGaaaca                                                      >  1:1711164/1‑90 (MQ=255)
                                   gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACaa                                                      >  2:3236065/1‑93 (MQ=255)
                                    atatAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACaaa                                                     <  1:5015984/93‑1 (MQ=255)
                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    >  1:1646720/1‑93 (MQ=255)
                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    <  1:5061593/93‑1 (MQ=255)
                                      ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGt                                                   >  1:2424052/1‑93 (MQ=255)
                                       taAAAATCGTCAAACAATAGCCAGGCGAGGAAATTCAAAGAGGGACTTTTTTGAATGCGCTTCAACATCACGGTGTTTCCATCAAACAAAgac                                                  >  2:5702579/1‑91 (MQ=255)
                                       taAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTg                                                  >  2:1881321/1‑93 (MQ=255)
                                       taAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTg                                                  >  1:1125537/1‑93 (MQ=255)
                                       taAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTg                                                  >  1:1254167/1‑93 (MQ=255)
                                         aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgcg                                                >  2:308922/1‑93 (MQ=255)
                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCa                                              >  2:3313965/1‑93 (MQ=255)
                                            aGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCat                                             <  1:1531440/93‑1 (MQ=255)
                                             gCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCata                                            <  1:1186333/93‑1 (MQ=255)
                                               ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATg                                          <  1:1730506/93‑1 (MQ=255)
                                               ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATg                                          <  1:4710958/93‑1 (MQ=255)
                                                gCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGc                                         >  1:400491/1‑93 (MQ=255)
                                                gCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGc                                         <  2:2570589/93‑1 (MQ=255)
                                                gCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGc                                         <  2:3874554/93‑1 (MQ=255)
                                                 ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCt                                        <  1:143476/93‑1 (MQ=255)
                                                  caacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTc                                       <  2:6172218/93‑1 (MQ=255)
                                                   aacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCg                                      >  1:5812733/1‑93 (MQ=255)
                                                   aacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCg                                      >  1:1032608/1‑93 (MQ=255)
                                                    acaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgc                                     >  1:5325643/1‑93 (MQ=255)
                                                    acaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgc                                     <  1:2082712/93‑1 (MQ=255)
                                                     caaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgca                                    <  2:5751972/93‑1 (MQ=255)
                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  >  1:2651802/1‑93 (MQ=255)
                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  <  1:3330966/93‑1 (MQ=255)
                                                         aGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagc                                <  1:3505611/93‑1 (MQ=255)
                                                          gCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagca                               >  2:108084/1‑93 (MQ=255)
                                                            cAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagaaaa                             <  1:6101830/93‑1 (MQ=255)
                                                               gCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAAtt                          >  2:4695985/1‑93 (MQ=255)
                                                                 gAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTaa                        <  1:2691988/93‑1 (MQ=255)
                                                                 gAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTaa                        >  1:2239929/1‑93 (MQ=255)
                                                                  aGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAg                       >  2:1337428/1‑93 (MQ=255)
                                                                   gTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGc                      >  2:3351414/1‑93 (MQ=255)
                                                                   gTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGc                      <  2:1368967/93‑1 (MQ=255)
                                                                       aTTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAc                  <  1:3729704/93‑1 (MQ=255)
                                                                       aTTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAc                  <  1:4069306/93‑1 (MQ=255)
                                                                        ttCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAcc                 <  2:1126227/93‑1 (MQ=255)
                                                                          cAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCtt               >  1:3161653/1‑93 (MQ=255)
                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             <  1:5390472/93‑1 (MQ=255)
                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             <  1:5494916/93‑1 (MQ=255)
                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             <  1:2512224/93‑1 (MQ=255)
                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             <  1:4906645/93‑1 (MQ=255)
                                                                             agagGGTCTTTTTGGTAGGCGCTTCAACATCACGGCGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaat            >  1:192735/1‑92 (MQ=255)
                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           >  2:1033089/1‑93 (MQ=255)
                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           >  1:2824637/1‑93 (MQ=255)
                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           <  1:4569227/93‑1 (MQ=255)
                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          >  2:846310/1‑93 (MQ=255)
                                                                                  gTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCtt       <  2:517934/93‑1 (MQ=255)
                                                                                  gTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCtt       <  2:2732979/93‑1 (MQ=255)
                                                                                  gTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCtt       <  1:596843/93‑1 (MQ=255)
                                                                                    cTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa     <  2:5642880/93‑1 (MQ=255)
                                                                                    cTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa     <  2:2315030/93‑1 (MQ=255)
                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTCAGTTCTTAAg    >  2:5795844/1‑93 (MQ=255)
                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    >  2:2773996/1‑93 (MQ=255)
                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    >  1:862870/1‑93 (MQ=255)
                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGGGCGCATATGTTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    >  2:5080417/1‑93 (MQ=255)
                                                                                          |                                                                                         
ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241

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