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% 1693.1 / NA 524L16P (CTA→CCA) yjdBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (288/236);  total (288/236)

Only 100 of 557 total aligned reads displayed.

ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241
                                                                                          |                                                                                         
acCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggt                                                                                         <  1:53112/93‑1 (MQ=255)
 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCACAGAGGGTCTTTTtggta                                                                                        >  2:4243100/1‑93 (MQ=255)
 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        <  2:8375040/93‑1 (MQ=255)
   tGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAgg                                                                                      >  1:9355221/1‑93 (MQ=255)
   tGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAgg                                                                                      <  1:5002248/93‑1 (MQ=255)
     tttAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGgcg                                                                                    >  1:6642978/1‑93 (MQ=255)
         aaaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTc                                                                                >  2:236107/1‑93 (MQ=255)
           aaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCaa                                                                              >  1:7049304/1‑93 (MQ=255)
           aaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCaa                                                                              <  2:953106/93‑1 (MQ=255)
             aGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACa                                                                            >  1:1430894/1‑93 (MQ=255)
             aGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACa                                                                            >  2:9667424/1‑93 (MQ=255)
             aGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACa                                                                            <  2:9260518/93‑1 (MQ=255)
              ggCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACAt                                                                           >  1:4279019/1‑93 (MQ=255)
                 aaTATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCAc                                                                        >  2:5856380/1‑93 (MQ=255)
                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      >  2:9178014/1‑93 (MQ=255)
                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      >  2:4471585/1‑93 (MQ=255)
                    atTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgt                                                                     <  1:468809/93‑1 (MQ=255)
                    atTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTATTTGGTAGGCGCTTCAACATGACGgt                                                                     >  1:6909565/1‑93 (MQ=255)
                      tCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtgt                                                                   <  1:2964936/93‑1 (MQ=255)
                       cAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGtt                                                                  >  1:2063578/1‑93 (MQ=255)
                        aGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGttt                                                                 <  2:4876819/93‑1 (MQ=255)
                         ggCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTc                                                                >  2:3501423/1‑93 (MQ=255)
                         ggCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTc                                                                >  2:6497525/1‑93 (MQ=255)
                          gCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTcc                                                               >  1:6549572/1‑93 (MQ=255)
                            aGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCAt                                                             <  2:9841707/93‑1 (MQ=255)
                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            >  1:612934/1‑93 (MQ=255)
                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            >  2:3274632/1‑93 (MQ=255)
                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            >  2:1401596/1‑93 (MQ=255)
                               tAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGa                                                          <  2:3097350/93‑1 (MQ=255)
                               tAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGa                                                          >  2:6458232/1‑93 (MQ=255)
                                agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGaa                                                         >  2:3678209/1‑93 (MQ=255)
                                agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGaa                                                         <  2:2660848/93‑1 (MQ=255)
                                agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGaa                                                         >  2:7040946/1‑93 (MQ=255)
                                agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGaa                                                         <  1:3211981/93‑1 (MQ=255)
                                  agaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCTCGGTGTTTCCATCGAACa                                                       <  2:6024733/93‑1 (MQ=255)
                                   gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACaa                                                      >  2:5295170/1‑93 (MQ=255)
                                    atatAAAAAGCGGCAAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTGTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACaaa                                                     <  1:7742681/93‑1 (MQ=255)
                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    <  1:9145170/93‑1 (MQ=255)
                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    >  1:8098950/1‑93 (MQ=255)
                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    <  1:1564072/93‑1 (MQ=255)
                                       taAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTg                                                  <  2:1907884/93‑1 (MQ=255)
                                        aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgc                                                 >  2:7571264/1‑93 (MQ=255)
                                        aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgc                                                 >  1:188149/1‑93 (MQ=255)
                                         aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgcg                                                <  1:7613965/93‑1 (MQ=255)
                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCGTCGAACAAAGTGCGCa                                              <  2:7661512/93‑1 (MQ=255)
                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCa                                              <  2:1163401/93‑1 (MQ=255)
                                            aGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCat                                             >  1:9851319/1‑93 (MQ=255)
                                             gCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCata                                            <  1:10015909/93‑1 (MQ=255)
                                              cGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCatat                                           <  2:1320594/93‑1 (MQ=255)
                                               ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATg                                          <  2:8875190/93‑1 (MQ=255)
                                               ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATg                                          >  1:7932474/1‑93 (MQ=255)
                                               ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATg                                          <  1:4904860/93‑1 (MQ=255)
                                               ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTGCCATCGAACAAAGTGCGCATATg                                          >  2:5518417/1‑93 (MQ=255)
                                                gCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGc                                         <  1:3532474/93‑1 (MQ=255)
                                                 ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCt                                        <  1:3778139/93‑1 (MQ=255)
                                                   aacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCg                                      >  1:3398843/1‑93 (MQ=255)
                                                    acaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgc                                     >  1:9560537/1‑93 (MQ=255)
                                                     caaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgca                                    >  2:8686284/1‑93 (MQ=255)
                                                     caaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgca                                    >  1:3898329/1‑93 (MQ=255)
                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  <  1:8676590/93‑1 (MQ=255)
                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 <  1:1130188/93‑1 (MQ=255)
                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 <  1:5905849/93‑1 (MQ=255)
                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 <  1:6487085/93‑1 (MQ=255)
                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 >  2:397979/1‑93 (MQ=255)
                                                          gCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagca                               >  2:6141991/1‑93 (MQ=255)
                                                            cAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagcaaa                             >  1:4524705/1‑93 (MQ=255)
                                                            cAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagcaaa                             >  1:7312205/1‑93 (MQ=255)
                                                             aGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCaaaa                            >  2:5165191/1‑93 (MQ=255)
                                                             aGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCaaaa                            >  1:1617956/1‑93 (MQ=255)
                                                             aGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCaaaa                            <  1:5207524/93‑1 (MQ=255)
                                                               gCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAAtt                          <  1:9448778/93‑1 (MQ=255)
                                                                 gAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTaa                        <  2:8098949/93‑1 (MQ=255)
                                                                 gAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTaa                        <  1:2627049/93‑1 (MQ=255)
                                                                  aGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCCTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAg                       <  1:5116463/93‑1 (MQ=255)
                                                                   gTAAATTCAAAGAGTGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGAATATGCTCGCAAAGCAAAATTAAGc                      <  1:8811665/93‑1 (MQ=255)
                                                                   gTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGc                      <  2:4347024/93‑1 (MQ=255)
                                                                    tAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGcc                     >  2:341615/1‑93 (MQ=255)
                                                                    tAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGcc                     >  2:20339/1‑93 (MQ=255)
                                                                      aaTTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCaa                   <  2:4566732/93‑1 (MQ=255)
                                                                       aTTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAc                  >  2:7739847/1‑93 (MQ=255)
                                                                        ttCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAcc                 >  1:5677614/1‑93 (MQ=255)
                                                                         tCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCt                >  1:8008791/1‑93 (MQ=255)
                                                                          cAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCtt               >  1:5518519/1‑93 (MQ=255)
                                                                          cAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCtt               >  2:565605/1‑93 (MQ=255)
                                                                          cAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCtt               <  1:7549397/93‑1 (MQ=255)
                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             >  2:7460245/1‑93 (MQ=255)
                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             <  1:7451705/93‑1 (MQ=255)
                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           >  1:9643967/1‑93 (MQ=255)
                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           >  2:413081/1‑93 (MQ=255)
                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           >  2:2430677/1‑93 (MQ=255)
                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           >  2:2923908/1‑93 (MQ=255)
                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          <  2:7150670/93‑1 (MQ=255)
                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          >  1:9007508/1‑93 (MQ=255)
                                                                                gggTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATTCTCGCAAAGCAAAATTAAGCCAACCTTAAGATc         >  2:8950069/1‑93 (MQ=255)
                                                                                gggTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTc         >  2:6853280/1‑93 (MQ=255)
                                                                                gggTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTc         <  2:6643271/93‑1 (MQ=255)
                                                                                  gTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCtt       >  2:4404110/1‑93 (MQ=255)
                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    >  1:2504943/1‑93 (MQ=255)
                                                                                      ttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAgg   >  2:734051/1‑93 (MQ=255)
                                                                                       tttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt  <  1:8239118/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241

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