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,7180GC100.0% 786.1 / NA 240A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (134/106);  total (134/106)
*NC_0009134,007,7190CG99.6% 787.0 / NA 240A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (134/105);  minor base A (0/1);  total (134/106)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

CCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCTTATGTAC  >  NC_000913/4007632‑4007811
                                                                                      ||                                                                                            
ccgccgCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAg                                                                                         >  2:3712764/1‑93 (MQ=255)
   ccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTcg                                                                                      >  2:4730430/1‑93 (MQ=255)
        ggCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTc                                                                                 <  2:2932130/93‑1 (MQ=255)
         gacATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTcc                                                                                <  1:6282048/91‑1 (MQ=255)
           cATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCACTCAGTCGCGCTCCAc                                                                              <  1:1126629/93‑1 (MQ=255)
            aTCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACg                                                                             <  1:5702275/93‑1 (MQ=255)
            aTCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACg                                                                             >  1:2812949/1‑93 (MQ=255)
             tCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGc                                                                            >  1:1102108/1‑93 (MQ=255)
             tCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGc                                                                            >  2:2421753/1‑93 (MQ=255)
             tCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGc                                                                            >  2:5057186/1‑93 (MQ=255)
                tGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                                         >  2:3883491/1‑93 (MQ=255)
                   cGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGt                                                                      <  1:881812/93‑1 (MQ=255)
                    gAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                                     <  1:1999022/93‑1 (MQ=255)
                     aaGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTg                                                                    >  1:5860380/1‑93 (MQ=255)
                     aaGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTg                                                                    <  2:4491234/93‑1 (MQ=255)
                     aaGGAGGACGGCCGGTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTg                                                                    >  1:4063695/1‑93 (MQ=255)
                      aggaggACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGa                                                                   >  1:6207585/1‑93 (MQ=255)
                      aggaggACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGa                                                                   <  2:5315668/93‑1 (MQ=255)
                      aggaggACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGa                                                                   <  1:471363/93‑1 (MQ=255)
                       ggaggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAt                                                                  <  2:2588576/93‑1 (MQ=255)
                        gaggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtt                                                                 <  1:2357945/93‑1 (MQ=255)
                        gaggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtt                                                                 >  2:2806699/1‑93 (MQ=255)
                         aggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttt                                                                >  1:1920335/1‑93 (MQ=255)
                          ggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttt                                                               >  2:2256115/1‑93 (MQ=255)
                           gaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttttt                                                              >  2:2089865/1‑93 (MQ=255)
                            aCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttttt                                                             <  2:1795503/93‑1 (MQ=255)
                             cGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTc                                                            <  2:5983271/93‑1 (MQ=255)
                              ggCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttttttta                                                           >  1:2462650/1‑91 (MQ=255)
                                 cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                        <  1:527410/93‑1 (MQ=255)
                                  gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                                       >  2:5908580/1‑93 (MQ=255)
                                  gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                                       <  2:2101302/93‑1 (MQ=255)
                                  gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                                       <  1:3669111/93‑1 (MQ=255)
                                   ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAgg                                                      >  1:4174169/1‑93 (MQ=255)
                                      gTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                                                   >  2:3008479/1‑93 (MQ=255)
                                       tAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATa                                                  <  2:2645705/93‑1 (MQ=255)
                                         aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                                <  1:3533651/93‑1 (MQ=255)
                                         aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                                >  1:1817169/1‑93 (MQ=255)
                                         aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                                >  2:4786086/1‑93 (MQ=255)
                                         aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                                >  1:2734024/1‑93 (MQ=255)
                                         aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                                >  2:5609021/1‑93 (MQ=255)
                                          ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                               >  1:4978259/1‑93 (MQ=255)
                                          ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCAATCGTTGATTTTTTCAACAGGCATAACt                                               >  1:1595148/1‑93 (MQ=255)
                                           tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                                              >  2:377306/1‑93 (MQ=255)
                                           tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                                              >  2:6431178/1‑93 (MQ=255)
                                            aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCa                                             <  2:3231162/93‑1 (MQ=255)
                                             aaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                                            >  2:80893/1‑93 (MQ=255)
                                             aaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                                            <  2:2344116/93‑1 (MQ=255)
                                              aGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAcc                                           >  1:4385135/1‑93 (MQ=255)
                                                gTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCa                                         <  1:3496519/93‑1 (MQ=255)
                                                  gCTTTCGATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCGCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGc                                       <  1:3213221/93‑1 (MQ=255)
                                                  gCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGc                                       <  1:1375971/93‑1 (MQ=255)
                                                   cTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCg                                      <  2:953444/93‑1 (MQ=255)
                                                     tACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGa                                    >  2:4769481/1‑93 (MQ=255)
                                                     tACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGa                                    <  1:3138797/93‑1 (MQ=255)
                                                       ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaaa                                  <  1:1925963/93‑1 (MQ=255)
                                                       ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaaa                                  >  1:6580099/1‑93 (MQ=255)
                                                        cgtGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                                 >  1:6411930/3‑93 (MQ=255)
                                                        cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                                 <  1:904657/93‑1 (MQ=255)
                                                        cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                                 >  2:3180128/1‑93 (MQ=255)
                                                        cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                                 <  1:1493640/93‑1 (MQ=255)
                                                         aTGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc                                >  2:1883377/1‑93 (MQ=255)
                                                           gagaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGc                              <  2:1541087/93‑1 (MQ=255)
                                                             gaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctc                            <  2:3621185/93‑1 (MQ=255)
                                                              aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct                           <  1:4928375/93‑1 (MQ=255)
                                                              aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCCGAAACCGctct                           >  1:2922633/1‑93 (MQ=255)
                                                               tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGCAACCGCTCTa                          >  1:2135315/1‑93 (MQ=255)
                                                               tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                          >  2:3382353/1‑93 (MQ=255)
                                                               tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                          >  1:4564881/1‑93 (MQ=255)
                                                                ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTAc                         >  1:1008986/1‑93 (MQ=255)
                                                                ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTAc                         >  1:5490552/1‑93 (MQ=255)
                                                                 cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                        <  1:3985224/93‑1 (MQ=255)
                                                                 cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                        >  2:3344971/1‑93 (MQ=255)
                                                                 cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATAGTTGATTTTTTCAACAGGCATAACTCACCCAGCCGAAACCGCTCTACa                        >  2:6630071/1‑93 (MQ=255)
                                                                  tttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACag                       <  1:6088068/93‑1 (MQ=255)
                                                                  tttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACag                       >  2:1673339/1‑93 (MQ=255)
                                                                  tttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACag                       >  1:5111091/1‑93 (MQ=255)
                                                                  tttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAGCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGGAACCGCTCTACag                       >  2:3561562/1‑93 (MQ=255)
                                                                   ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACaga                      >  1:2698690/1‑93 (MQ=255)
                                                                    tttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACagag                     <  2:4609953/93‑1 (MQ=255)
                                                                    tttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACagag                     >  1:5822/1‑93 (MQ=255)
                                                                      tGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGt                   <  2:1405455/93‑1 (MQ=255)
                                                                      tGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGt                   <  1:226361/93‑1 (MQ=255)
                                                                      tGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGt                   <  2:6080732/93‑1 (MQ=255)
                                                                        aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGttt                 <  1:699623/93‑1 (MQ=255)
                                                                        aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGttt                 >  1:3808585/1‑93 (MQ=255)
                                                                        aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGttt                 <  1:3753604/93‑1 (MQ=255)
                                                                         aaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTa                >  1:3337497/1‑93 (MQ=255)
                                                                          aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaa               >  1:4822358/1‑93 (MQ=255)
                                                                          aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaa               <  2:4144851/93‑1 (MQ=255)
                                                                           aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa              >  1:6043064/1‑93 (MQ=255)
                                                                           aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa              >  2:5528895/1‑93 (MQ=255)
                                                                           aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa              >  2:6184386/1‑93 (MQ=255)
                                                                            ggACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAt             >  2:3968938/1‑93 (MQ=255)
                                                                            ggACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAt             >  2:2770508/1‑93 (MQ=255)
                                                                             gACGTAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAtt            <  2:1152604/93‑1 (MQ=255)
                                                                             gACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAtt            >  2:2901629/1‑93 (MQ=255)
                                                                              aCGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAttt           >  2:6260488/1‑93 (MQ=255)
                                                                                gCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCt         <  1:3574714/93‑1 (MQ=255)
                                                                                gCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCt         <  1:4319080/93‑1 (MQ=255)
                                                                                       ccTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCTTATGTAc  <  1:5376466/93‑1 (MQ=255)
                                                                                      ||                                                                                            
CCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCTTATGTAC  >  NC_000913/4007632‑4007811

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