Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,170,208 T→C G191G (GGT→GGC)  murB → UDP‑N‑acetylenolpyruvoylglucosamine reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,170,2080TC99.6% 768.8 / ‑5.4 227G191G (GGT→GGC) murBUDP‑N‑acetylenolpyruvoylglucosamine reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (113/112);  total (115/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 247 total aligned reads displayed.

TATTTTTAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCC  >  NC_000913/4170118‑4170297
                                                                                          |                                                                                         
tATTTTTAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGa                                                                                         >  2:6744409/1‑93 (MQ=255)
tATTTTTAAACATGAATACCAGGACCGCTTCGCTATTGGAGCCGTAGGTCTGCGGCTGCCAAAAGAGAGGCAACCTGCACTAACGTAGGGCGa                                                                                         >  1:6684831/1‑93 (MQ=255)
    tttAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTa                                                                                     <  1:4565947/93‑1 (MQ=255)
     ttAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTaa                                                                                    >  2:4605184/1‑93 (MQ=255)
        aaCATGAATCCCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTc                                                                                 <  2:8421652/93‑1 (MQ=255)
         aCATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCg                                                                                >  1:4490794/1‑93 (MQ=255)
         aCATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCg                                                                                >  2:6044507/1‑93 (MQ=255)
          cATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGt                                                                               <  2:3125558/93‑1 (MQ=255)
          cATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGCGGCAACCTGTACTAACGTATGGCGACTTAACTCGt                                                                               >  2:3014614/1‑93 (MQ=255)
           aTGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTc                                                                              >  2:4938300/1‑93 (MQ=255)
            tGAATACCAGTACCGCTTCGCTATGGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCACCCTGTACTAACGTATGGCGACTTAACTCGTCt                                                                             <  2:2391140/93‑1 (MQ=255)
             gAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTg                                                                            <  2:5764300/93‑1 (MQ=255)
              aaTACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTgg                                                                           <  2:375857/93‑1 (MQ=255)
               aTACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGa                                                                          >  1:3626692/1‑93 (MQ=255)
               aTACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGa                                                                          <  1:2025189/93‑1 (MQ=255)
                tACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGAt                                                                         >  2:4382020/1‑93 (MQ=255)
                  ccAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATcc                                                                       <  1:8123301/93‑1 (MQ=255)
                   cAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCt                                                                      >  2:786559/1‑93 (MQ=255)
                   cAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCt                                                                      >  1:1131641/1‑93 (MQ=255)
                     ggACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTac                                                                    <  1:3703984/93‑1 (MQ=255)
                     ggACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTac                                                                    >  2:2732350/1‑93 (MQ=255)
                     ggACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTac                                                                    <  2:7614049/93‑1 (MQ=255)
                     ggACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTac                                                                    >  1:7611450/1‑93 (MQ=255)
                         cGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTacaaca                                                                <  2:1902115/93‑1 (MQ=255)
                         cGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTaaaaca                                                                >  2:1521852/1‑93 (MQ=255)
                          gCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAg                                                               <  2:665254/93‑1 (MQ=255)
                          gCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAg                                                               >  1:678193/1‑93 (MQ=255)
                          gCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGAAACCTGTACTCACGTATGGCGACTTAACCGGTCTGGATCCTccaacat                                                               >  1:1273186/1‑92 (MQ=255)
                           cTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGt                                                              >  1:3492476/1‑93 (MQ=255)
                           cTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGt                                                              <  1:8759149/93‑1 (MQ=255)
                           cTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGt                                                              >  2:2863999/1‑93 (MQ=255)
                           cTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGt                                                              >  2:7260567/1‑93 (MQ=255)
                            ttCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTa                                                             >  2:2765618/1‑93 (MQ=255)
                             tCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTaa                                                            <  1:4000996/93‑1 (MQ=255)
                             tCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTaa                                                            <  1:1567776/93‑1 (MQ=255)
                              cGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAAc                                                           <  2:2521303/93‑1 (MQ=255)
                               gCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACg                                                          >  1:5246238/1‑93 (MQ=255)
                                cTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGc                                                         >  1:1441511/1‑93 (MQ=255)
                                cTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGc                                                         >  2:8200092/1‑93 (MQ=255)
                                  aTTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCca                                                       >  2:4077166/1‑93 (MQ=255)
                                     ttAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAg                                                    <  1:7166045/92‑1 (MQ=255)
                                     gTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAg                                                    >  2:2113498/1‑93 (MQ=255)
                                       aGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCa                                                  <  1:1976718/93‑1 (MQ=255)
                                        gCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCaa                                                 >  2:1741166/1‑93 (MQ=255)
                                         ccGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAg                                                >  2:4747303/1‑93 (MQ=255)
                                          cGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGt                                               <  1:826640/93‑1 (MQ=255)
                                          cGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGt                                               <  2:3932702/93‑1 (MQ=255)
                                            tAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTAt                                             <  2:7369851/93‑1 (MQ=255)
                                              ggTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTAttt                                           >  1:3686212/1‑93 (MQ=255)
                                              agtctgcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTAttt                                           <  1:2182555/92‑1 (MQ=255)
                                                tctgcgtctgcCAGAGGAGTGGCCGCCTGTACTGACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTACCGCCACAGCAAGTAttttg                                         <  1:7802328/93‑3 (MQ=255)
                                                tctgcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTaa                                         >  2:7122125/1‑93 (MQ=255)
                                                tctgcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATCTaa                                         <  2:5124557/93‑1 (MQ=255)
                                                  tgcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATg                                       <  2:5395753/93‑1 (MQ=255)
                                                  tgcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATg                                       >  2:7231150/1‑93 (MQ=255)
                                                   gcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGc                                      >  2:7944752/1‑93 (MQ=255)
                                                   gcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGc                                      <  1:8501222/93‑1 (MQ=255)
                                                   gcgtctgcCAAAAGAGAGGCAACCTGTACTAACGTATGGCGTCTTGACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGc                                      <  1:7845082/93‑1 (MQ=255)
                                                    cgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCg                                     <  2:656304/93‑1 (MQ=255)
                                                    cgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCg                                     <  2:3570464/93‑1 (MQ=255)
                                                     gtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCgg                                    <  1:5574128/93‑1 (MQ=255)
                                                     gtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCgg                                    <  1:5387433/93‑1 (MQ=255)
                                                      tctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGgt                                   <  2:1215316/93‑1 (MQ=255)
                                                      tctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGgt                                   >  1:2910619/1‑93 (MQ=255)
                                                       ctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGgtg                                  >  1:1090445/1‑93 (MQ=255)
                                                       ctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGgtg                                  <  2:622443/93‑1 (MQ=255)
                                                          ccAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGgtgtgt                               >  2:2325341/1‑93 (MQ=255)
                                                           cAAAAGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTc                              <  1:2436266/93‑1 (MQ=255)
                                                            aaaaGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCa                             >  1:8617497/1‑93 (MQ=255)
                                                            aaaaGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCa                             >  1:5680311/1‑93 (MQ=255)
                                                             aaaGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCat                            <  1:8021002/93‑1 (MQ=255)
                                                              aaGAGTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCata                           <  1:5874054/93‑1 (MQ=255)
                                                               agagTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCatat                          <  2:1109588/93‑1 (MQ=255)
                                                                gagTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATg                         >  1:6368663/1‑93 (MQ=255)
                                                                gagTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATg                         >  1:6026427/1‑93 (MQ=255)
                                                                  gTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATgcg                       <  1:3849622/93‑1 (MQ=255)
                                                                  gTGGCAACCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATgcg                       >  2:5619608/1‑93 (MQ=255)
                                                                    ggCAACCTGTACTAACGTATGGGGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGca                     >  1:7118234/1‑93 (MQ=255)
                                                                       aaCCTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGcacca                  <  1:5605605/93‑1 (MQ=255)
                                                                        aCCTGTACTAACATATGGCGACGTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGcacccc                 <  2:8671529/93‑3 (MQ=255)
                                                                         ccTGTCCTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCaccacc                <  2:132270/93‑1 (MQ=255)
                                                                         ccTGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCaccacc                >  2:2902571/1‑93 (MQ=255)
                                                                           tGTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCaa              <  2:6827334/93‑1 (MQ=255)
                                                                            gTACTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCaaa             <  2:8144780/93‑1 (MQ=255)
                                                                               cTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTc          <  2:5705643/93‑1 (MQ=255)
                                                                               cTAACGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTc          <  2:4477579/93‑1 (MQ=255)
                                                                                 aaCGTATGGCGACTTAACTCGTCTGGATCCTACAAGAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTccc        >  1:2269703/1‑93 (MQ=255)
                                                                                 aaCGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTccc        <  2:3494919/93‑1 (MQ=255)
                                                                                 aaCGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTccc        >  2:5198499/1‑93 (MQ=255)
                                                                                 aaCGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTccc        <  2:7728421/93‑1 (MQ=255)
                                                                                 aaCGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTccc        <  1:4201852/93‑1 (MQ=255)
                                                                                 aaCGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTccc        >  2:8650581/1‑93 (MQ=255)
                                                                                  aCGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTACCGCCACAGCAAGTATTTAATGCGGTGTGTCATATCCGCACCCCCAAACTCCCt       <  2:2420756/93‑1 (MQ=255)
                                                                                  aCGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCt       >  2:8338487/1‑93 (MQ=255)
                                                                                  aCGTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCt       <  1:6975586/93‑1 (MQ=255)
                                                                                    gTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGa     <  1:2469852/93‑1 (MQ=255)
                                                                                    gTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGa     >  1:2840866/1‑93 (MQ=255)
                                                                                    gTATGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGa     >  2:5643330/1‑93 (MQ=255)
                                                                                      aTGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATc   >  2:6143862/1‑93 (MQ=255)
                                                                                       tGGCGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAACAAGCATTTAATGCGGTGTGTCATATACGCACCACCAAACTCCCCGCTcc  >  1:7506250/1‑93 (MQ=255)
                                                                                          |                                                                                         
TATTTTTAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCC  >  NC_000913/4170118‑4170297

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