Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,990,664 A→G S59S (AGT→AGC pgsA ← phosphatidylglycerophosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,990,6640AG98.6% 428.7 / ‑2.8 138S59S (AGT→AGCpgsAphosphatidylglycerophosphate synthetase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (92/44);  total (94/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 142 total aligned reads displayed.

CAGCTGTGATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCA  >  NC_000913/1990575‑1990753
                                                                                         |                                                                                         
caGCTGTGATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGctt                                                                                        >  1:279716/1‑92 (MQ=255)
 aGCTGTGATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGcttt                                                                                       >  2:2292182/1‑91 (MQ=255)
 aGCTGTGATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGcttt                                                                                       <  2:4525030/93‑3 (MQ=255)
  gCTGTGATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTg                                                                                      >  1:3018897/1‑93 (MQ=255)
   cTGTGATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTgg                                                                                     <  2:4189484/93‑1 (MQ=255)
    tgtgATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGt                                                                                    >  1:412937/1‑93 (MQ=255)
      tgATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGATCTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTc                                                                                  <  2:2764815/93‑1 (MQ=255)
      tcaTAATGCTCGGTTACCAGCACCATGGCGATAGCCGCGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAGCCGGGTGCTTTGGTTc                                                                                  <  1:2228554/91‑1 (MQ=255)
       taTAATGCTCGGTTACCAGCACCATGGCGATATCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTcc                                                                                 <  1:1056818/92‑1 (MQ=255)
       gATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTcc                                                                                 <  1:1056817/93‑1 (MQ=255)
        aTAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCa                                                                                >  2:1515694/1‑93 (MQ=255)
        aTAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCa                                                                                >  1:3466997/1‑93 (MQ=255)
         tAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAg                                                                               >  1:923231/1‑93 (MQ=255)
               tCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGt                                                                         >  1:882722/1‑93 (MQ=255)
               tCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGt                                                                         >  1:231027/1‑93 (MQ=255)
               tCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGt                                                                         <  1:3127776/93‑1 (MQ=255)
                cGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTg                                                                        >  1:454397/1‑93 (MQ=255)
                  gTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGcc                                                                      >  2:4136139/1‑93 (MQ=255)
                    tACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAg                                                                    >  1:1927075/1‑93 (MQ=255)
                    tACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAg                                                                    >  2:4739118/1‑93 (MQ=255)
                        aGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGcaaa                                                                >  2:4592173/1‑93 (MQ=255)
                        aGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGaaaa                                                                <  2:2790682/93‑1 (MQ=255)
                        aGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGaaaa                                                                >  2:4185329/1‑93 (MQ=255)
                         gCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAAc                                                               >  1:532014/1‑93 (MQ=255)
                           aCCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCa                                                             >  2:732644/1‑93 (MQ=255)
                           aCCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCa                                                             >  2:3578799/1‑93 (MQ=255)
                            ccATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCGAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCAt                                                            >  1:4678653/1‑93 (MQ=255)
                            ccATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCAt                                                            <  1:2133802/93‑1 (MQ=255)
                              aTGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCg                                                          >  2:3273109/1‑93 (MQ=255)
                               tGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGa                                                         >  2:2224823/1‑93 (MQ=255)
                                ggCGATAGCCCCGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGaa                                                        <  2:3613387/93‑1 (MQ=255)
                                 gCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAAc                                                       >  1:2647512/1‑93 (MQ=255)
                                   gATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACca                                                     <  2:3761804/93‑1 (MQ=255)
                                   gATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACca                                                     <  2:2350993/93‑1 (MQ=255)
                                    aTAGCCACGAGAATTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTCGGATCCATCGGCGTGCCAGAAAATCATCGAACcag                                                    >  2:4194761/1‑93 (MQ=255)
                                    aTAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcag                                                    >  2:1200499/1‑93 (MQ=255)
                                     tAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagt                                                   >  1:1431527/1‑93 (MQ=255)
                                     tAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagt                                                   >  1:1431526/1‑93 (MQ=255)
                                     tAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCAAACcagt                                                   >  2:2213821/1‑93 (MQ=255)
                                      aGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagtc                                                  <  1:2839356/93‑1 (MQ=255)
                                      aGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagtc                                                  <  2:3678227/93‑1 (MQ=255)
                                      aGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagtc                                                  >  1:1354457/1‑93 (MQ=255)
                                      aGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagtc                                                  >  1:2642431/1‑93 (MQ=255)
                                        ccACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagtcag                                                >  1:1830701/1‑93 (MQ=255)
                                         cACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagtcagt                                               <  1:2516753/93‑1 (MQ=255)
                                         cACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagtcagt                                               <  2:1311615/93‑1 (MQ=255)
                                          aCGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCAGCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACcagtcagta                                              >  2:3191500/1‑92 (MQ=255)
                                          aCGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCagtcagtc                                              >  1:3188960/1‑93 (MQ=255)
                                          aCGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCagtcagtc                                              >  1:137327/1‑93 (MQ=255)
                                          aCGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCagtcagtc                                              >  2:412117/1‑93 (MQ=255)
                                           cGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAgtcagtca                                             >  2:4161886/1‑93 (MQ=255)
                                           cGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAgtcagtca                                             <  2:4309838/93‑1 (MQ=255)
                                             agaACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCAcc                                           <  2:4699693/93‑1 (MQ=255)
                                               aaCTTTATCTGCCACAGGGTCAAGGAACGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGCACCAGTCAGTCAccgc                                         >  2:2842033/1‑93 (MQ=255)
                                               aaCTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCAccgc                                         >  1:2735247/1‑93 (MQ=255)
                                                aCTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCAccgcc                                        >  1:2031521/1‑93 (MQ=255)
                                                 cTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCAccgccg                                       >  2:2322917/1‑93 (MQ=255)
                                                  tttATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACcgccgc                                      >  2:3071385/1‑93 (MQ=255)
                                                  tttATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACcgccgc                                      >  2:430752/1‑93 (MQ=255)
                                                  tttATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACcgccgc                                      >  1:4738977/1‑93 (MQ=255)
                                                   ttATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCcgcg                                     >  1:3782924/1‑93 (MQ=255)
                                                    tATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGa                                    <  1:4724907/93‑1 (MQ=255)
                                                       cTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGc                                 >  2:4006200/1‑93 (MQ=255)
                                                          ccACAGGGTCAAGGAAAGCACCAGACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGa                              <  1:1730286/93‑1 (MQ=255)
                                                           cacaGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGaa                             >  1:3995955/1‑93 (MQ=255)
                                                           cacaGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGaa                             >  1:116794/1‑93 (MQ=255)
                                                            acaGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGaaa                            >  2:4436065/1‑93 (MQ=255)
                                                             caGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCATTAACCGCCGCGACGCAGaaaa                           >  1:1669103/1‑93 (MQ=255)
                                                             caGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGaaaa                           <  2:3275956/93‑1 (MQ=255)
                                                              aGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAAt                          >  2:481961/1‑93 (MQ=255)
                                                              aGGGTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAAt                          <  2:3757799/93‑1 (MQ=255)
                                                               gggTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATg                         >  1:2521654/1‑93 (MQ=255)
                                                               gggTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATg                         >  1:44623/1‑93 (MQ=255)
                                                                 gTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAg                       <  1:2117282/93‑1 (MQ=255)
                                                                 gTCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAg                       >  2:169372/1‑93 (MQ=255)
                                                                  tcgaGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAgc                      >  2:3941232/4‑93 (MQ=255)
                                                                  tCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGACCCAGTCAGTCACCGCCGCGACGCAGAAAATGAgc                      <  1:4348967/93‑1 (MQ=255)
                                                                  tCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAgc                      <  2:67625/93‑1 (MQ=255)
                                                                  tCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAgc                      >  2:3250957/1‑93 (MQ=255)
                                                                  tCAAGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAgc                      >  1:501450/1‑93 (MQ=255)
                                                                    aaGGAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAgcgc                    <  1:1509507/93‑1 (MQ=255)
                                                                      ggAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCgcgg                  >  1:36791/1‑93 (MQ=255)
                                                                      ggAAAGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCgcgg                  >  2:2312818/1‑93 (MQ=255)
                                                                        aaaGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCgcggcg                >  1:4233820/1‑93 (MQ=255)
                                                                         aaGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGcggcgg               <  1:2974610/93‑1 (MQ=255)
                                                                         aaGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGcggcgg               >  2:4085741/1‑93 (MQ=255)
                                                                         aaGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGcggcgg               <  2:2598743/93‑1 (MQ=255)
                                                                          aGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAGCCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCggcggc              >  1:2341994/1‑93 (MQ=255)
                                                                          aGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCggcggc              >  1:205341/1‑93 (MQ=255)
                                                                          aGCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCggcggc              >  1:4672466/1‑93 (MQ=255)
                                                                           gCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCa             >  2:3499707/1‑93 (MQ=255)
                                                                           gCACCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCa             >  2:3065819/1‑93 (MQ=255)
                                                                             aCCAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCaaa           >  1:1746564/1‑93 (MQ=255)
                                                                              ccAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAAc          <  2:4049652/93‑1 (MQ=255)
                                                                              ccAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAAc          <  2:3414812/93‑1 (MQ=255)
                                                                               cAAACCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACg         <  1:3960222/93‑1 (MQ=255)
                                                                                 aaCCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGc       <  2:1905542/93‑1 (MQ=255)
                                                                                  aCCGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCg      >  2:461283/1‑93 (MQ=255)
                                                                                    cGGGTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGAc    <  1:613293/93‑1 (MQ=255)
                                                                                     gggTGCTTTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGAcc   <  2:3118403/93‑1 (MQ=255)
                                                                                         |                                                                                         
CAGCTGTGATAATGCTCGGTTACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACCGGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCA  >  NC_000913/1990575‑1990753

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: