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% 1929.8 / NA 554A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (297/257);  total (297/257)
*NC_0009134,007,7190CG99.8% 1930.5 / NA 553A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (296/256);  minor base T (0/1);  total (296/257)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 567 total aligned reads displayed.

ACCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCTT  >  NC_000913/4007631‑4007805
                                                                                       ||                                                                                      
aCCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCa                                                                                    <  1:525623/93‑1 (MQ=255)
 ccgccgCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGGAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAg                                                                                   >  2:40720/1‑93 (MQ=255)
  cgccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGt                                                                                  <  1:1031543/93‑1 (MQ=255)
   gccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                                 >  1:7156043/1‑93 (MQ=255)
   gccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                                 <  1:1855225/93‑1 (MQ=255)
   gccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                                 >  2:6126919/1‑93 (MQ=255)
   gccgcGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGGGCCTACCATGAGATCCTTTTTGAAAAGGAAGCAATGTGCTCAGTc                                                                                 >  2:15223/1‑93 (MQ=255)
   gccgcGGGCCACCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAATGGACGCAATGCGCTCAGTc                                                                                 <  1:4653974/93‑1 (MQ=255)
       cgGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCgcgc                                                                             <  2:2244258/93‑1 (MQ=255)
       cgGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCgcgc                                                                             <  1:2341977/93‑1 (MQ=255)
        gggCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCt                                                                            >  2:467150/1‑93 (MQ=255)
         ggCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTc                                                                           >  2:5001758/1‑93 (MQ=255)
         ggCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTc                                                                           >  1:5696195/1‑93 (MQ=255)
         ggCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTc                                                                           >  2:3086459/1‑93 (MQ=255)
           ccATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCa                                                                         <  2:273070/93‑1 (MQ=255)
            cATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                                        >  1:8185697/1‑93 (MQ=255)
                 tGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                                   <  1:120843/93‑1 (MQ=255)
                  gTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATc                                                                  >  2:3734920/1‑93 (MQ=255)
                  gTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATc                                                                  >  2:5521424/1‑93 (MQ=255)
                   tCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCg                                                                 >  2:1669674/1‑93 (MQ=255)
                   gcGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCg                                                                 <  1:1515620/92‑1 (MQ=255)
                    cGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGt                                                                <  2:1188190/93‑1 (MQ=255)
                     gAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                               <  1:7300982/93‑1 (MQ=255)
                     gAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                               <  1:1748647/93‑1 (MQ=255)
                      aaGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTg                                                              >  2:8349157/1‑93 (MQ=255)
                       aggaggACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGa                                                             >  2:8733925/1‑93 (MQ=255)
                         gaggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtt                                                           <  1:208668/93‑1 (MQ=255)
                            gaCGGCCGTTGGTAATTAAAGGTGCTTACATGAGATTCTTTTTGAAAAGGACGCAATGCGCTCAATCGCGCTCCACGCCCTCGTTGGtttttt                                                       >  1:1334676/1‑93 (MQ=255)
                              cGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTc                                                      >  2:6467058/1‑93 (MQ=255)
                               ggCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCa                                                     >  2:1959776/1‑93 (MQ=255)
                                 gcGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAAc                                                   <  2:5896851/92‑1 (MQ=255)
                                 ccGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAAc                                                   >  1:2648365/1‑93 (MQ=255)
                                 ccGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAAc                                                   <  2:1801862/93‑1 (MQ=255)
                                  cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTGAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                  <  1:1458499/93‑1 (MQ=255)
                                  cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                  <  1:1922106/93‑1 (MQ=255)
                                  cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                                  >  1:399868/1‑93 (MQ=255)
                                   gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                                 >  2:7670375/1‑93 (MQ=255)
                                    ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAgg                                                >  2:7171388/1‑93 (MQ=255)
                                    ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAgg                                                >  2:428541/1‑93 (MQ=255)
                                     tGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                                               >  2:5142771/1‑93 (MQ=255)
                                       gTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                                             <  2:6272879/93‑1 (MQ=255)
                                       gTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                                             <  2:7492497/93‑1 (MQ=255)
                                       gTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                                             <  2:1015743/93‑1 (MQ=255)
                                         aaTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATaa                                           <  1:3147379/93‑1 (MQ=255)
                                          aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                          >  1:3655872/1‑93 (MQ=255)
                                          aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                                          <  1:5349750/93‑1 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                         >  2:8019744/1‑93 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                         >  1:3797745/1‑93 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                         <  2:3338376/93‑1 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                         <  2:3874190/93‑1 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                         <  2:6194606/93‑1 (MQ=255)
                                           ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                                         >  1:6170249/1‑93 (MQ=255)
                                            tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                                        >  2:3543106/1‑93 (MQ=255)
                                            tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                                        >  2:3013402/1‑93 (MQ=255)
                                            tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                                        >  1:3444100/1‑93 (MQ=255)
                                             aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCa                                       <  1:5461208/93‑1 (MQ=255)
                                               aGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAcc                                     <  1:8030549/93‑1 (MQ=255)
                                                ggTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                                    <  1:2979904/93‑1 (MQ=255)
                                                 gtacTTAGCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCTCCCa                                   <  1:3596818/90‑1 (MQ=255)
                                                  tGCTTACCATGAGATCCTTTTTGAGAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAg                                  <  2:573095/93‑1 (MQ=255)
                                                   gCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGc                                 >  2:8526681/1‑93 (MQ=255)
                                                    cTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCg                                >  1:561633/1‑93 (MQ=255)
                                                    cTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCg                                <  1:3710468/93‑1 (MQ=255)
                                                    cTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCg                                <  1:216381/93‑1 (MQ=255)
                                                     ttACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCgg                               >  2:241377/1‑93 (MQ=255)
                                                     ttACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCgg                               >  2:6623007/1‑93 (MQ=255)
                                                       aCCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaa                             >  1:7761456/1‑93 (MQ=255)
                                                       aCCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaa                             <  1:4157211/93‑1 (MQ=255)
                                                        ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaaa                            <  2:2123609/93‑1 (MQ=255)
                                                         cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                           >  1:8742721/1‑93 (MQ=255)
                                                         cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                           >  1:830394/1‑93 (MQ=255)
                                                         cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                           <  1:7930833/93‑1 (MQ=255)
                                                         cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                           <  2:8586154/93‑1 (MQ=255)
                                                          aTGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc                          <  1:7421230/93‑1 (MQ=255)
                                                          aTGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc                          >  1:5864897/1‑93 (MQ=255)
                                                            gagaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGc                        >  1:6722846/1‑93 (MQ=255)
                                                             agaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGct                       <  2:7229404/93‑1 (MQ=255)
                                                              gaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctc                      <  1:8425835/93‑1 (MQ=255)
                                                               aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct                     <  2:7555167/93‑1 (MQ=255)
                                                               aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct                     <  1:6546476/93‑1 (MQ=255)
                                                               aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct                     <  1:7117331/93‑1 (MQ=255)
                                                                tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                    >  2:5381751/1‑93 (MQ=255)
                                                                tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                    >  2:5257410/1‑93 (MQ=255)
                                                                tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa                    >  2:2616795/1‑93 (MQ=255)
                                                                    ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACaga                >  1:6076998/1‑93 (MQ=255)
                                                                        gAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGtt            <  1:686624/93‑1 (MQ=255)
                                                                         aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGttt           <  2:205035/93‑1 (MQ=255)
                                                                          aaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTa          >  1:4513018/1‑93 (MQ=255)
                                                                           aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAATCACCCAGCGGAAACCGCTCTACAGAGGTTTaa         >  1:8493735/1‑93 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa        >  2:7771392/1‑93 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa        >  2:8163295/1‑93 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa        <  1:4827911/93‑1 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa        >  2:7366658/1‑93 (MQ=255)
                                                                            aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTaaa        >  1:5052480/1‑93 (MQ=255)
                                                                             ggACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAt       >  2:6839566/1‑93 (MQ=255)
                                                                             ggACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAt       >  1:3269445/1‑93 (MQ=255)
                                                                              gACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAtt      >  2:4590171/1‑93 (MQ=255)
                                                                              gACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAtt      >  1:4228308/1‑93 (MQ=255)
                                                                               aCGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAAttt     <  2:2928618/93‑1 (MQ=255)
                                                                                  cAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCCCCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCtt  <  1:6413077/93‑1 (MQ=255)
                                                                                       ||                                                                                      
ACCGCCGCGGGCCATCCTGTCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGTTTAAATTTCTT  >  NC_000913/4007631‑4007805

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