Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,355 A→G H314R (CAC→CGC)  metL → aspartokinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,3550AG100.0% 589.5 / NA 186H314R (CAC→CGC) metLaspartokinase II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (95/91);  total (95/91)

Only 100 of 203 total aligned reads displayed.

AACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCG  >  NC_000913/4128265‑4128444
                                                                                          |                                                                                         
aaCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCcgcg                                                                                         >  2:2807998/1‑93 (MQ=255)
 aCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCga                                                                                        >  1:2382536/1‑93 (MQ=255)
 aCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCga                                                                                        <  2:2434349/93‑1 (MQ=255)
  cTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgat                                                                                       <  1:1448314/93‑1 (MQ=255)
    gcgcTGTAGCTACACGCGGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgatga                                                                                     <  1:371818/93‑1 (MQ=255)
    gcgcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgatga                                                                                     >  2:4345446/1‑93 (MQ=255)
    gcgcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgatga                                                                                     <  2:2226297/93‑1 (MQ=255)
    gcgcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgatga                                                                                     <  1:1800652/93‑1 (MQ=255)
      gcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGatgatg                                                                                   >  2:2465690/1‑93 (MQ=255)
       cTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGt                                                                                  <  1:1279471/93‑1 (MQ=255)
       cTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGt                                                                                  >  2:3616308/1‑93 (MQ=255)
         gTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCt                                                                                >  1:1321422/1‑93 (MQ=255)
         gTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCt                                                                                >  1:1978149/1‑93 (MQ=255)
         gTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCt                                                                                >  2:128971/1‑93 (MQ=255)
           aGCTACACGCCGGATCAAGGTTCGACGCGCATTGAACGCGTGCTGGCCTCTGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGt                                                                              <  1:2520538/93‑1 (MQ=255)
            gCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGtt                                                                             >  1:2308043/1‑93 (MQ=255)
            gCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGtt                                                                             >  1:4305691/1‑93 (MQ=255)
             cTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGttt                                                                            >  1:3913723/1‑93 (MQ=255)
             cTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGttt                                                                            <  2:4476561/93‑1 (MQ=255)
              tACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttg                                                                           >  2:1133575/1‑93 (MQ=255)
               acacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttga                                                                          <  2:395873/93‑1 (MQ=255)
               acacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttga                                                                          <  1:1343678/93‑1 (MQ=255)
                  cGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgattg                                                                       <  2:1658582/93‑1 (MQ=255)
                  cGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgattg                                                                       <  2:4665135/93‑1 (MQ=255)
                   gCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgattga                                                                      <  1:679276/93‑1 (MQ=255)
                    ccGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAg                                                                     <  1:3279276/93‑1 (MQ=255)
                    ccGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAg                                                                     <  1:3136344/93‑1 (MQ=255)
                     cGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGt                                                                    <  2:612167/93‑1 (MQ=255)
                     cGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGt                                                                    <  1:2535675/93‑1 (MQ=255)
                      ggATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGtt                                                                   >  1:1932894/1‑93 (MQ=255)
                      ggATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGtt                                                                   <  1:3703434/93‑1 (MQ=255)
                       gATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGttt                                                                  <  2:1051197/93‑1 (MQ=255)
                        aTCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTc                                                                 <  1:4523954/93‑1 (MQ=255)
                        aTCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTc                                                                 <  2:2509405/93‑1 (MQ=255)
                        aTCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGCATTGTCACCAGCCGCGATGATGCCTGTTTGATTGAGTTTc                                                                 >  2:3802199/1‑93 (MQ=255)
                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                >  1:509904/1‑93 (MQ=255)
                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                <  1:1020776/93‑1 (MQ=255)
                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                >  1:1369419/1‑93 (MQ=255)
                          cAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAg                                                               <  2:3207647/93‑1 (MQ=255)
                          cAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAg                                                               <  1:2468479/93‑1 (MQ=255)
                            aGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGt                                                             >  2:2770565/1‑93 (MQ=255)
                            aGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGt                                                             <  1:582541/93‑1 (MQ=255)
                            aGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGt                                                             <  2:913662/93‑1 (MQ=255)
                               ttCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGcc                                                          <  2:669804/93‑1 (MQ=255)
                               ttCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGcc                                                          >  2:1320003/1‑93 (MQ=255)
                                tCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGGTTCAGGTGccc                                                         >  2:1623997/1‑93 (MQ=255)
                                  cACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGc                                                       >  1:2455625/1‑93 (MQ=255)
                                   ccgcgCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGcc                                                      <  1:1591914/92‑1 (MQ=255)
                                   aCGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGcc                                                      >  1:280787/1‑93 (MQ=255)
                                   aCGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGcc                                                      <  2:3369701/93‑1 (MQ=255)
                                    cgcgCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCGCCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCa                                                     <  1:1398906/93‑1 (MQ=255)
                                    cgcgCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCa                                                     >  1:3415153/1‑93 (MQ=255)
                                      cgcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAgc                                                   >  2:4532743/1‑92 (MQ=255)
                                      cgcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGt                                                   >  1:631474/1‑93 (MQ=255)
                                       gcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTc                                                  >  1:182879/1‑93 (MQ=255)
                                        cATTGAACGCCGCCTGGCCTCCGGTACTGGTGCGCGTATTGTCGCCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCa                                                 <  1:2170148/93‑1 (MQ=255)
                                         aTTGAGCGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAg                                                <  2:2038630/93‑1 (MQ=255)
                                         aTTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAg                                                >  2:1038231/1‑93 (MQ=255)
                                          ttGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAgg                                               <  2:3952841/93‑1 (MQ=255)
                                           tGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGa                                              >  2:4232501/1‑93 (MQ=255)
                                             aaCGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAtt                                            >  1:967742/1‑93 (MQ=255)
                                                  gTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCaaa                                       <  2:2364086/93‑1 (MQ=255)
                                                  gTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCaaa                                       >  1:818556/1‑93 (MQ=255)
                                                    gCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACt                                     <  1:3093885/93‑1 (MQ=255)
                                                     ttcgCCCCCGGTACTGGGGCGCGTATTGGCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTg                                    <  1:1145558/90‑1 (MQ=255)
                                                     cTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTg                                    >  1:78402/1‑93 (MQ=255)
                                                      tGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTgg                                   >  2:1860634/1‑93 (MQ=255)
                                                      tGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTgg                                   >  1:4093531/1‑93 (MQ=255)
                                                      tGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTgg                                   >  1:1752448/1‑93 (MQ=255)
                                                      tGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTgg                                   >  2:771444/1‑93 (MQ=255)
                                                       ggCCTCCGGTCGTGGTGCGTGTATTGCAACCAGCCGCGATGATGTCTGTTTAATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgc                                  <  2:4276998/93‑1 (MQ=255)
                                                       ggCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgc                                  >  1:2145970/1‑93 (MQ=255)
                                                       gcccTCCGGTACTGGGGCGCGTATGGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgc                                  <  1:2970261/91‑1 (MQ=255)
                                                          cTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCa                               <  1:1846685/93‑1 (MQ=255)
                                                            ccGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATa                             <  2:1583613/93‑1 (MQ=255)
                                                            ccGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATa                             <  1:2065250/93‑1 (MQ=255)
                                                            ccGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATa                             <  1:573470/93‑1 (MQ=255)
                                                            ccGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATa                             <  2:1934317/93‑1 (MQ=255)
                                                             cGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATaa                            <  2:3697754/93‑1 (MQ=255)
                                                              ggTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATaaa                           <  2:1356937/93‑1 (MQ=255)
                                                              ggTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATaaa                           <  1:1064571/93‑1 (MQ=255)
                                                              ggTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATaaa                           >  1:4263759/1‑93 (MQ=255)
                                                               gTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAag                          >  2:1877225/1‑93 (MQ=255)
                                                                tACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAaga                         <  1:3047902/93‑1 (MQ=255)
                                                                tACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAaga                         >  2:2703160/1‑93 (MQ=255)
                                                                       gcgcGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGAc                  <  1:3879159/93‑1 (MQ=255)
                                                                       gcgcGTATGGTCACCAGCCGGGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGAc                  <  2:3298497/93‑1 (MQ=255)
                                                                         gcgTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCCAACTGGCGCATAACGAGATAGACCa                >  1:234284/1‑93 (MQ=255)
                                                                            tATTGTCACCAGCCGCGATGATGTGTGTTTGATGGAGTTGCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAAt             <  2:4050640/93‑1 (MQ=255)
                                                                            tATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAAt             >  1:1720892/1‑93 (MQ=255)
                                                                            tATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAAt             >  1:647111/1‑93 (MQ=255)
                                                                            tATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAAt             >  2:4585957/1‑93 (MQ=255)
                                                                              ttGTCGCCAGCCGCGAGGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATcc           <  2:1100851/93‑1 (MQ=255)
                                                                              ttGTCACCACCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCGGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATcc           <  1:2517421/93‑1 (MQ=255)
                                                                               tGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCt          <  2:3157692/93‑1 (MQ=255)
                                                                                 tCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGa        >  2:2636211/1‑93 (MQ=255)
                                                                                  cACCAGCCGCGATGATGTCTGTTTGATTGAGGTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGGGCATAAAGAGATCGACCCAATCCTGaa       >  2:417773/1‑93 (MQ=255)
                                                                                    ccAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAc     <  2:1785447/93‑1 (MQ=255)
                                                                                      aGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgc   >  2:2186643/1‑93 (MQ=255)
                                                                                       gCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgcg  >  2:3188629/1‑93 (MQ=255)
                                                                                          |                                                                                         
AACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCG  >  NC_000913/4128265‑4128444

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