Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,007,718 2 bp→CG coding (943‑944/1023 nt) pldB → lysophospholipase L(2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,7180GC99.6% 767.7 / NA 237A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (136/100);  total (137/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,007,7190CG100.0% 771.7 / NA 236A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (136/100);  total (136/100)

Only 100 of 243 total aligned reads displayed.

ACCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCTT  >  NC_000913/4007631‑4007805
                                                                                       ||                                                                                      
   gccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                                 <  2:2313121/93‑1 (MQ=255)
   gccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                                 >  2:5051122/1‑93 (MQ=255)
   gccgcGGGCCATCCTGTCGAAGGAGGACGGCCGCTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                                 >  2:49076/1‑93 (MQ=255)
    ccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGtgg                                                                                >  1:4863847/1‑91 (MQ=255)
    ccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTcg                                                                                >  1:1916834/1‑93 (MQ=255)
        gggCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCt                                                                            >  2:537858/1‑93 (MQ=255)
         ggCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTc                                                                           <  1:3026019/93‑1 (MQ=255)
            cATCCTGGTGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                                        <  2:3491071/93‑1 (MQ=255)
             aTCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACg                                                                       >  2:4120017/1‑93 (MQ=255)
             aTCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACg                                                                       >  1:2128854/1‑93 (MQ=255)
               ccTGTCGAAGGAGGACGGCCGTTGGCAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGcc                                                                     <  1:6087473/93‑1 (MQ=255)
                cTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCa                                                                    <  2:3305087/93‑1 (MQ=255)
                 tGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                                   <  2:5881832/93‑1 (MQ=255)
                 tGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                                   <  1:5160366/93‑1 (MQ=255)
                     gAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                               >  1:1429368/1‑93 (MQ=255)
                     gAAGGAGGACGGCCGTTGGTAATTAAAGGAGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                               >  1:158269/1‑93 (MQ=255)
                     gAAGGAGGACGGCCGTTGGGAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                               >  1:5729117/1‑93 (MQ=255)
                        ggaggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAt                                                            >  2:4483253/1‑93 (MQ=255)
                          aggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttt                                                          <  1:835073/93‑1 (MQ=255)
                           ggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATACTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttt                                                         >  1:6294372/1‑93 (MQ=255)
                            gaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttttt                                                        >  1:1330617/1‑93 (MQ=255)
                              cGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTc                                                      <  1:3546847/93‑1 (MQ=255)
                              cGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTc                                                      >  2:6020319/1‑93 (MQ=255)
                                gCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCaa                                                    >  2:20776/1‑93 (MQ=255)
                                gCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCaa                                                    <  1:5179516/93‑1 (MQ=255)
                                 ccGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAAc                                                   >  1:5042593/1‑93 (MQ=255)
                                  cgctGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                  <  2:4746979/90‑1 (MQ=255)
                                  cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                  >  2:2382596/1‑93 (MQ=255)
                                  cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                  <  2:3008346/93‑1 (MQ=255)
                                  cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                  >  1:3239425/1‑93 (MQ=255)
                                   gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                                 >  2:2937751/1‑93 (MQ=255)
                                   gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                                 <  1:886723/93‑1 (MQ=255)
                                   gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                                 <  1:3161099/93‑1 (MQ=255)
                                    ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAgg                                                >  2:1957497/1‑93 (MQ=255)
                                    ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAgg                                                >  1:686434/1‑93 (MQ=255)
                                     tGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                                               >  2:1290880/1‑93 (MQ=255)
                                     tGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                                               <  1:1157305/93‑1 (MQ=255)
                                     tGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                                               <  1:5257609/93‑1 (MQ=255)
                                      ggTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCa                                              >  2:5938370/1‑93 (MQ=255)
                                      ggTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCa                                              >  2:1726593/1‑93 (MQ=255)
                                       gTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                                             <  2:1609010/93‑1 (MQ=255)
                                        tAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATa                                            <  2:3894545/93‑1 (MQ=255)
                                         aaTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATaa                                           >  2:1149022/1‑93 (MQ=255)
                                         aaTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATaa                                           >  2:1059865/1‑93 (MQ=255)
                                          aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                          <  2:2645810/93‑1 (MQ=255)
                                          aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                          >  1:1523833/1‑93 (MQ=255)
                                          aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                          <  2:6151726/93‑1 (MQ=255)
                                          aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                          <  1:478084/93‑1 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                         <  2:3728009/93‑1 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                         <  1:4439601/93‑1 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                         >  2:5815243/1‑93 (MQ=255)
                                            tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                                        >  2:4607143/1‑93 (MQ=255)
                                             aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCa                                       <  2:4347775/93‑1 (MQ=255)
                                              aaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                                      >  1:5308785/1‑93 (MQ=255)
                                              aaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                                      >  1:4469209/1‑93 (MQ=255)
                                              aaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                                      <  1:253460/93‑1 (MQ=255)
                                              aaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                                      >  1:4104197/1‑93 (MQ=255)
                                                ggTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                                    >  1:617190/1‑93 (MQ=255)
                                                ggTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                                    >  1:6073239/1‑93 (MQ=255)
                                                ggTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                                    >  1:4064028/1‑93 (MQ=255)
                                                ggTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                                    <  2:2773454/93‑1 (MQ=255)
                                                 gTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCa                                   <  2:2561647/93‑1 (MQ=255)
                                                  tGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAg                                  >  2:3978002/1‑93 (MQ=255)
                                                    cTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCg                                <  1:6039978/93‑1 (MQ=255)
                                                     ttACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCgg                               >  2:6165402/1‑93 (MQ=255)
                                                     ttACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCgg                               >  1:2852839/1‑93 (MQ=255)
                                                         cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                           >  2:5779542/1‑93 (MQ=255)
                                                         cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                           <  2:1781255/93‑1 (MQ=255)
                                                         cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                           >  1:5102631/1‑93 (MQ=255)
                                                          aTGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc                          >  2:4261195/1‑93 (MQ=255)
                                                            gagaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGc                        <  2:4967077/93‑1 (MQ=255)
                                                             agaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGct                       >  1:4396369/1‑93 (MQ=255)
                                                             agaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGct                       <  2:120782/93‑1 (MQ=255)
                                                               aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct                     <  1:2056403/93‑1 (MQ=255)
                                                               aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct                     <  2:1808405/93‑1 (MQ=255)
                                                                tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctctc                    >  2:480369/1‑92 (MQ=255)
                                                                tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                    <  1:5084002/93‑1 (MQ=255)
                                                                tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                    >  2:4893363/1‑93 (MQ=255)
                                                                tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                    >  1:2475309/1‑93 (MQ=255)
                                                                 ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTAc                   >  1:1098619/1‑93 (MQ=255)
                                                                 ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTAc                   >  1:2789801/1‑93 (MQ=255)
                                                                  cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                  <  2:702919/93‑1 (MQ=255)
                                                                   tttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACag                 >  2:3054275/1‑93 (MQ=255)
                                                                   tttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACag                 >  1:4253948/1‑93 (MQ=255)
                                                                    ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACata                >  2:2800875/1‑91 (MQ=255)
                                                                    ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACaga                <  1:5563672/93‑1 (MQ=255)
                                                                    ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACaga                >  2:5546987/1‑93 (MQ=255)
                                                                     tttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACagag               <  1:3418738/93‑1 (MQ=255)
                                                                     tttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACagag               <  2:5182044/93‑1 (MQ=255)
                                                                         aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGttt           >  1:3706056/1‑93 (MQ=255)
                                                                          aaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTa          >  2:3204036/1‑93 (MQ=255)
                                                                           aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaa         <  1:710468/93‑1 (MQ=255)
                                                                           aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaa         >  2:840699/1‑93 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa        >  2:3650817/1‑93 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa        >  1:3778312/1‑93 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa        <  1:5621999/93‑1 (MQ=255)
                                                                             ggACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTACAt       >  1:3667465/1‑93 (MQ=255)
                                                                             ggACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAt       >  1:2331640/1‑93 (MQ=255)
                                                                                cGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTc    <  1:5891562/93‑1 (MQ=255)
                                                                                  cAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCtt  <  1:5399643/93‑1 (MQ=255)
                                                                                       ||                                                                                      
ACCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCTT  >  NC_000913/4007631‑4007805

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