Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,629,559 T→C Y420H (TAC→CAC)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,629,5590TC84.8% 675.0 / 101.5 263Y420H (TAC→CAC) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base C (107/116);  total (126/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 266 total aligned reads displayed.

AAGAGTATGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGGCGCGGTAT  >  NC_000913/4629467‑4629648
                                                                                            |                                                                                         
    gTATGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCTAGTTGATGTACTg                                                                                       >  1:463762/1‑93 (MQ=255)
    gTATGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTg                                                                                       <  1:5008990/93‑1 (MQ=255)
     tATGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTgg                                                                                      >  1:3301752/1‑93 (MQ=255)
      aTGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTGAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGa                                                                                     <  2:3561616/93‑1 (MQ=255)
      aTGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGa                                                                                     <  2:3286927/93‑1 (MQ=255)
       tGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGaa                                                                                    >  2:624187/1‑93 (MQ=255)
          gCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAAtct                                                                                 >  2:4714005/1‑93 (MQ=255)
           cTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAAtctc                                                                                <  2:4536540/93‑1 (MQ=255)
            tGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCg                                                                               >  2:4094430/1‑93 (MQ=255)
            tGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGCCGCGCGGGCGCGAGTTGATGCACTGTAATCTCg                                                                               >  2:2562182/1‑93 (MQ=255)
              aaGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGat                                                                             >  1:2993303/1‑93 (MQ=255)
               aGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGata                                                                            >  1:3199233/1‑93 (MQ=255)
                gATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGAATCTCGataa                                                                           >  1:5739248/1‑93 (MQ=255)
                gATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGataa                                                                           >  2:3079291/1‑93 (MQ=255)
                gATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGataa                                                                           >  1:5875425/1‑93 (MQ=255)
                 aTTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGataat                                                                          <  1:4125527/93‑1 (MQ=255)
                 aTTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGataat                                                                          <  2:6527305/93‑1 (MQ=255)
                 aTTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGataat                                                                          <  2:1375754/93‑1 (MQ=255)
                 aTTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGataat                                                                          >  1:3337826/1‑93 (MQ=255)
                  ttGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGTCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGataata                                                                         <  1:529957/93‑1 (MQ=255)
                  ttGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGataata                                                                         >  2:5504433/1‑93 (MQ=255)
                   tGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGAATCTCGATAATAc                                                                        <  1:5418101/93‑1 (MQ=255)
                   tGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATAc                                                                        <  1:1370946/93‑1 (MQ=255)
                    gATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATAcc                                                                       <  2:5808467/93‑1 (MQ=255)
                      tAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACcgc                                                                     >  2:5107168/1‑93 (MQ=255)
                      tAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACcgc                                                                     >  1:2494988/1‑93 (MQ=255)
                      tAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACcgc                                                                     <  1:6632467/93‑1 (MQ=255)
                       aaaGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACcgcg                                                                    >  1:6261203/1‑93 (MQ=255)
                        aaGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACcgcgc                                                                   >  1:1057452/1‑93 (MQ=255)
                         aGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGAATCTCGATAATACcgcgcg                                                                  >  1:6049977/1‑93 (MQ=255)
                          gcgcCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGt                                                                 >  1:196740/1‑93 (MQ=255)
                            gccgcAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGGGATGGCGCGCGTGCGCGAGTTGATGCACTGGCATCTCGATAATACCGCGCGCAg                                                               >  1:5591010/1‑93 (MQ=255)
                            gccgcAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGAATCTCGATAATACCGCGCGTAg                                                               <  2:2899252/93‑1 (MQ=255)
                            gccgcAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAg                                                               <  2:5929203/93‑1 (MQ=255)
                            gccgcAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAg                                                               <  2:2992248/93‑1 (MQ=255)
                             ccgcAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGc                                                              >  2:752471/1‑93 (MQ=255)
                              cgcAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCg                                                             >  2:3813391/1‑93 (MQ=255)
                                cAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAg                                                           <  1:1707951/93‑1 (MQ=255)
                                 aGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGt                                                          <  1:4693870/93‑1 (MQ=255)
                                 aGAATGTTGACAGCGCCCTGACTCAGGGAGCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGt                                                          >  1:3226753/1‑93 (MQ=255)
                                  gAATGTTGCCAGCGCCCTGAATCGGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTg                                                         <  1:3910842/93‑1 (MQ=255)
                                  gAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTg                                                         <  2:1608058/93‑1 (MQ=255)
                                   aaTGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTgg                                                        <  2:5147157/93‑1 (MQ=255)
                                   aaTGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTgg                                                        >  1:5784386/1‑93 (MQ=255)
                                     tGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGc                                                      <  2:4957401/93‑1 (MQ=255)
                                      gTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGcc                                                     <  1:5363033/93‑1 (MQ=255)
                                      gTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGcc                                                     >  1:3653149/1‑93 (MQ=255)
                                        tGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGCGGGCCaa                                                   <  2:2954536/93‑1 (MQ=255)
                                         gACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAAt                                                  <  1:3721553/93‑1 (MQ=255)
                                          aCAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATc                                                 <  2:64408/93‑1 (MQ=255)
                                          aCAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATc                                                 <  2:2647852/93‑1 (MQ=255)
                                           cAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCt                                                <  1:491257/93‑1 (MQ=255)
                                            aGCGCGCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTg                                               <  2:5597582/93‑1 (MQ=255)
                                             gcgcCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTgg                                              <  2:5306958/93‑1 (MQ=255)
                                             gcgcCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTgg                                              >  2:3198599/1‑93 (MQ=255)
                                             gcgcCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTgg                                              <  2:4838551/93‑1 (MQ=255)
                                                cccTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGa                                           >  2:4660208/1‑93 (MQ=255)
                                                cccTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGa                                           >  2:4943500/1‑93 (MQ=255)
                                                cccTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGa                                           >  2:2139186/1‑93 (MQ=255)
                                                cccTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGa                                           >  1:3550327/1‑93 (MQ=255)
                                                cccTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGa                                           >  1:4391040/1‑93 (MQ=255)
                                                 ccTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGaa                                          >  2:3737292/1‑93 (MQ=255)
                                                 ccTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGaa                                          >  2:3389861/1‑93 (MQ=255)
                                                 ccTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGaa                                          >  2:5962738/1‑93 (MQ=255)
                                                  cTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAag                                         >  1:2095364/1‑93 (MQ=255)
                                                  cTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAag                                         <  2:4253856/93‑1 (MQ=255)
                                                   tGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAaga                                        <  2:5362457/93‑1 (MQ=255)
                                                   tGAATCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAaga                                        >  1:1952151/1‑93 (MQ=255)
                                                    gACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAagag                                       >  1:1289434/1‑93 (MQ=255)
                                                    gACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAagag                                       >  1:2400844/1‑93 (MQ=255)
                                                      cTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCa                                     >  1:1365131/1‑93 (MQ=255)
                                                       tCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCaa                                    <  1:2811374/93‑1 (MQ=255)
                                                         aGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGt                                  >  1:3618203/1‑93 (MQ=255)
                                                           ggggCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATCCCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCa                                <  1:5904614/93‑1 (MQ=255)
                                                            ggCCGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCaa                               >  1:6383989/1‑93 (MQ=255)
                                                             gCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCaaa                              >  1:2726632/1‑93 (MQ=255)
                                                               cGGAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCaaaaa                            <  2:3890152/93‑1 (MQ=255)
                                                                ggAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAAc                           <  2:4148727/93‑1 (MQ=255)
                                                                ggAGATGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAAc                           >  1:6216312/1‑93 (MQ=255)
                                                                  agaTGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACag                         >  1:3988432/1‑93 (MQ=255)
                                                                   gaTGGCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACaga                        <  1:1515207/93‑1 (MQ=255)
                                                                      ggCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCa                     <  2:1156229/93‑1 (MQ=255)
                                                                      ggCGCGCGTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCa                     <  1:5620639/93‑1 (MQ=255)
                                                                        cgcgcgTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAgg                   >  2:4454941/1‑93 (MQ=255)
                                                                            cgTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCa               <  2:3963366/93‑1 (MQ=255)
                                                                            cgTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCa               <  2:1263426/93‑1 (MQ=255)
                                                                             gTGCGCGAGTTGATGCACTGGAGTCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCaa              <  2:6650384/93‑1 (MQ=255)
                                                                             gTGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCaa              <  1:6333779/93‑1 (MQ=255)
                                                                              tGCGCGAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAAc             >  1:5231033/1‑93 (MQ=255)
                                                                                 gcgAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTgg          >  1:3116392/1‑93 (MQ=255)
                                                                                 gcgAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTgg          <  1:2673501/93‑1 (MQ=255)
                                                                                  cgAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGgc         >  1:1668771/1‑93 (MQ=255)
                                                                                   gAGTTGATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGgcg        >  2:1865293/1‑93 (MQ=255)
                                                                                    aGTTGATGTACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGgcgc       <  1:4750778/93‑1 (MQ=255)
                                                                                        gATGTACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGGCGCGGTa   <  1:929939/93‑1 (MQ=255)
                                                                                        gATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCACAAACAGAGCAGGCTCAACTGGCGCGGTa   >  2:2627594/1‑93 (MQ=255)
                                                                                        gATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGGCGCGGTa   <  1:3836383/93‑1 (MQ=255)
                                                                                        gATGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGGCGCGGTa   <  1:1139103/93‑1 (MQ=255)
                                                                                         aTGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGGCGCGGtgt  <  1:3465337/93‑3 (MQ=255)
                                                                                         aTGCACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGGCGCGGTAt  <  2:6256120/93‑1 (MQ=255)
                                                                                            |                                                                                         
AAGAGTATGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTGACTCAGGGGCCGGAGATGGCGCGCGTGCGCGAGTTGATGTACTGGAATCTCGATAATACCGCGCGTAGCGAGTGGGCCAATCTGGTGAAGAGCAAGTCAAAAACAGAGCAGGCTCAACTGGCGCGGTAT  >  NC_000913/4629467‑4629648

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