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.5% 759.0 / NA 209A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (105/103);  minor base A (0/1);  total (105/104)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,007,7190CG100.0% 768.1 / NA 209A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (105/104);  total (105/104)

Only 100 of 214 total aligned reads displayed.

ACCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCTTATGTA  >  NC_000913/4007631‑4007810
                                                                                       ||                                                                                           
aCCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCa                                                                                         >  1:545551/1‑93 (MQ=255)
   gccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                                      >  2:1153031/1‑93 (MQ=255)
   gccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                                      >  1:5424844/1‑93 (MQ=255)
    ccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTcg                                                                                     >  2:5341368/1‑93 (MQ=255)
     cgcgGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTcgc                                                                                    <  2:3134279/93‑1 (MQ=255)
          gCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTcc                                                                               <  1:4287167/93‑1 (MQ=255)
           ccATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCa                                                                              <  2:191370/93‑1 (MQ=255)
           ccATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCa                                                                              <  2:4162202/93‑1 (MQ=255)
            cATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                                             <  2:3898619/93‑1 (MQ=255)
             aTCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACg                                                                            >  1:1865818/1‑93 (MQ=255)
              tCCTGTCGAAGGAGGGCGGCCGTTGGTAATTAAAGGTGCTTAGCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGc                                                                           >  1:5533324/1‑93 (MQ=255)
              tCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGc                                                                           >  1:4069627/1‑93 (MQ=255)
              tCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGc                                                                           >  1:3546764/1‑93 (MQ=255)
                 tGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                                        <  2:5027051/93‑1 (MQ=255)
                 tGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                                        <  1:2796579/93‑1 (MQ=255)
                   tCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCg                                                                      >  1:704179/1‑93 (MQ=255)
                    cGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGt                                                                     >  2:3451312/1‑93 (MQ=255)
                    cGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGt                                                                     <  2:5894524/93‑1 (MQ=255)
                     gAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                                    <  1:4507105/93‑1 (MQ=255)
                      aaGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTg                                                                   >  2:3746565/1‑93 (MQ=255)
                      aaGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTg                                                                   >  2:5357342/1‑93 (MQ=255)
                      aaGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTg                                                                   >  1:2111790/1‑93 (MQ=255)
                        ggaggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAt                                                                 >  1:5557182/1‑93 (MQ=255)
                           ggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttt                                                              <  2:4321350/93‑1 (MQ=255)
                           ggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttt                                                              <  2:5210760/93‑1 (MQ=255)
                           ggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttt                                                              >  2:4747206/1‑93 (MQ=255)
                            gaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttttt                                                             <  1:1459835/93‑1 (MQ=255)
                              cGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTc                                                           <  2:1683135/93‑1 (MQ=255)
                                gCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCaa                                                         >  1:1004634/1‑93 (MQ=255)
                                gCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCaa                                                         <  1:1813081/93‑1 (MQ=255)
                                gCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCaa                                                         <  1:3066525/93‑1 (MQ=255)
                                  cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCGACGCCATCGTTGATTTTTTCAACa                                                       <  1:2628986/93‑1 (MQ=255)
                                  cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                       <  1:1514125/93‑1 (MQ=255)
                                  cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                       >  1:3298218/1‑93 (MQ=255)
                                   gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                                      <  1:3640743/93‑1 (MQ=255)
                                   gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                                      <  2:2488999/93‑1 (MQ=255)
                                    ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAgg                                                     >  2:1351722/1‑93 (MQ=255)
                                    ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAgg                                                     <  1:1595467/93‑1 (MQ=255)
                                    ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAgg                                                     >  2:2585481/1‑93 (MQ=255)
                                     tGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                                                    <  1:4676948/93‑1 (MQ=255)
                                     tGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                                                    <  2:5444371/93‑1 (MQ=255)
                                     tGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                                                    <  1:1306670/93‑1 (MQ=255)
                                      ggTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCa                                                   <  1:1905727/93‑1 (MQ=255)
                                       gTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                                                  <  2:5977087/93‑1 (MQ=255)
                                          aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                               >  1:1078968/1‑93 (MQ=255)
                                          aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                               <  2:2893757/93‑1 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                              <  1:369039/93‑1 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                              <  2:3695823/93‑1 (MQ=255)
                                            tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                                             <  1:455164/93‑1 (MQ=255)
                                            tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                                             >  1:4150729/1‑93 (MQ=255)
                                            tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                                             <  1:3430087/93‑1 (MQ=255)
                                             aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCa                                            <  1:87930/93‑1 (MQ=255)
                                             aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCa                                            <  1:4098064/93‑1 (MQ=255)
                                              aaGGTGCTTACCATGAGATCCTTTTTGATAAGGCCGCAATGAGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                                           <  2:5483173/93‑1 (MQ=255)
                                                ggTGCTTCCCAGGAGATCCTTTGTGAAAAGGACGCAATGCGCTCAGGCGCGCTCCCCGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                                         <  1:2611545/93‑1 (MQ=255)
                                                ggTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                                         >  2:803575/1‑93 (MQ=255)
                                                ggTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                                         >  1:3885573/1‑93 (MQ=255)
                                                   gCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGc                                      >  2:3083372/1‑93 (MQ=255)
                                                     ttACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCgg                                    >  2:2513464/1‑93 (MQ=255)
                                                        ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaaa                                 <  2:682737/93‑1 (MQ=255)
                                                         cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                                <  1:1229445/93‑1 (MQ=255)
                                                         cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                                <  2:3733125/93‑1 (MQ=255)
                                                          aTGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc                               <  2:3407339/93‑1 (MQ=255)
                                                           tGAGATCCTTTTTGAAACGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGCTTTTTTCAACAGGCCTAACTCACCCAGCGGCAACCg                              >  1:656761/1‑93 (MQ=255)
                                                           tGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCg                              >  1:160920/1‑93 (MQ=255)
                                                             agaTCCTTTTTGGAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCGTAACTCACCCAGCGGAAACCGct                            >  1:3494662/1‑93 (MQ=255)
                                                             agaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGct                            >  2:1001325/1‑93 (MQ=255)
                                                              gaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctc                           >  2:4699703/1‑93 (MQ=255)
                                                              gaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctc                           <  2:6076882/93‑1 (MQ=255)
                                                               aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCCGCGGAAACCGctct                          >  1:6121412/1‑93 (MQ=255)
                                                               aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct                          >  1:812952/1‑93 (MQ=255)
                                                               aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct                          >  2:2409349/1‑93 (MQ=255)
                                                               aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct                          <  2:737498/93‑1 (MQ=255)
                                                                tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                         >  2:1025695/1‑93 (MQ=255)
                                                                tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                         <  1:5885238/93‑1 (MQ=255)
                                                                 ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTAc                        <  1:4339859/93‑1 (MQ=255)
                                                                 ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTAc                        <  1:573636/93‑1 (MQ=255)
                                                                  cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                       <  1:5216711/93‑1 (MQ=255)
                                                                  cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                       <  2:2487074/93‑1 (MQ=255)
                                                                  cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                       <  1:2187484/93‑1 (MQ=255)
                                                                  cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                       <  2:5305239/93‑1 (MQ=255)
                                                                  cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                       >  2:3007212/1‑93 (MQ=255)
                                                                  cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                       >  2:4805540/1‑93 (MQ=255)
                                                                  cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACa                       >  1:3834027/1‑93 (MQ=255)
                                                                   tttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACag                      >  1:4805740/1‑93 (MQ=255)
                                                                   tttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACag                      <  1:2881953/93‑1 (MQ=255)
                                                                    ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACaga                     >  2:4627457/1‑93 (MQ=255)
                                                                    ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACaga                     <  1:3121424/93‑1 (MQ=255)
                                                                    ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACaga                     >  1:2296746/1‑93 (MQ=255)
                                                                    ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACaga                     <  2:1294899/93‑1 (MQ=255)
                                                                     tttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACagag                    <  2:2177162/93‑1 (MQ=255)
                                                                        gAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGtt                 <  2:4888208/93‑1 (MQ=255)
                                                                         aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGttt                <  2:2563029/93‑1 (MQ=255)
                                                                          aaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTGCAGAGGTTTa               >  2:3878378/1‑93 (MQ=255)
                                                                           aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaa              >  2:5102267/1‑93 (MQ=255)
                                                                           aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaa              >  2:2708949/1‑93 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa             <  2:4511000/93‑1 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa             >  2:1874920/1‑93 (MQ=255)
                                                                              gACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAtt           >  2:4465207/1‑93 (MQ=255)
                                                                                  cAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCtt       <  1:1968689/93‑1 (MQ=255)
                                                                                       ||                                                                                           
ACCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCTTATGTA  >  NC_000913/4007631‑4007810

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