Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,007,907 +GG intergenic (+109/‑11) pldB → / → b4555 lysophospholipase L(2)/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,9071.G96.5% 299.0 / 5.0 114intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (60/50);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,007,9072.G96.4% 294.9 / 5.0 113intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (59/49);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

AGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTA  >  NC_000913/4007832‑4007978
                                                                              ||                                                                       
aGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgc‑‑gcgc                                                                           <  2:5257995/76‑1 (MQ=255)
aGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgc‑‑gcgc                                                                           <  2:693761/76‑1 (MQ=255)
 gATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgc‑‑gcgcg                                                                          >  1:4936989/1‑75 (MQ=255)
 gATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgc‑‑gcgcg                                                                          <  1:3648541/76‑2 (MQ=255)
          ttACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACt                                                                 >  2:3135064/1‑76 (MQ=255)
           tACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACtt                                                                >  2:3508880/1‑76 (MQ=255)
            aCTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACttt                                                               >  1:2277270/1‑76 (MQ=255)
             cTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTCtg                                                              >  2:1044842/1‑76 (MQ=255)
              tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTtgt                                                             >  1:4728416/1‑76 (MQ=255)
               ttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTtgtg                                                            >  1:5711080/1‑76 (MQ=255)
               ttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTtgtg                                                            >  2:3110760/1‑76 (MQ=255)
                tctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTgtgt                                                           >  1:1161463/1‑76 (MQ=255)
                 ctcCCGACCATACGTTATCCCCTTACGCAAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGtt                                                          <  1:4146484/76‑1 (MQ=255)
                  tcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGttt                                                         <  2:4491958/76‑1 (MQ=255)
                    ccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGc                                                       >  2:5609480/1‑76 (MQ=255)
                    ccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGc                                                       <  1:1672332/76‑1 (MQ=255)
                    ccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGc                                                       >  2:5167519/1‑76 (MQ=255)
                    ccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGc                                                       >  1:3239327/1‑76 (MQ=255)
                    ccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGc                                                       <  2:4281528/76‑1 (MQ=255)
                    ccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGc                                                       <  1:3612559/76‑1 (MQ=255)
                    ccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGc                                                       <  1:3832651/76‑1 (MQ=255)
                     cGACCATACGTTATCCCCTTACGCCAAAGAAGCTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCg                                                      <  2:4725290/76‑1 (MQ=255)
                     cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCg                                                      <  1:1172553/76‑1 (MQ=255)
                     cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTATTTGCg                                                      <  1:1888208/76‑1 (MQ=255)
                      gACCATACGTTATCCCCTTACGCCAACGAACCTATGCAGCCGCTCACCGC‑‑GCGCGGCACACACTTTGTGTCTGCGa                                                     >  2:708115/1‑76 (MQ=255)
                      gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgc‑‑gcgcg                                                                          >  1:3142135/1‑54 (MQ=255)
                      gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgc‑‑gcgcg                                                                          <  2:3142135/55‑2 (MQ=255)
                      gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGa                                                     >  2:3249576/1‑76 (MQ=255)
                      gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGa                                                     >  2:705416/1‑76 (MQ=255)
                       aCCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGAc                                                    <  1:1013165/76‑1 (MQ=255)
                        ccATACGTTATCCCCTTGCGCCAAAGAAACTCTGAAGCTTCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGAcc                                                   <  1:1803390/76‑1 (MQ=255)
                        ccATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGAcc                                                   <  2:2587772/76‑1 (MQ=255)
                         cATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCg                                                  <  2:638879/76‑1 (MQ=255)
                         cATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCg                                                  >  1:2717372/1‑76 (MQ=255)
                          aTACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCgg                                                 >  2:4452639/1‑76 (MQ=255)
                            aCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGgtc                                               >  2:1607623/1‑76 (MQ=255)
                             ctttATCCCCTAACACCAAAGAAACCCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGgtcg                                              <  1:1181590/74‑1 (MQ=255)
                             cGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGgtcg                                              <  1:4691496/76‑1 (MQ=255)
                             cGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGgtcg                                              >  2:923065/1‑76 (MQ=255)
                             cGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGgtcg                                              <  1:5613866/76‑1 (MQ=255)
                              gTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGgtcgt                                             >  2:3228970/1‑76 (MQ=255)
                               tcaTCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGgtcgtc                                            <  1:3879513/74‑1 (MQ=255)
                                tATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGACATCAACTTTGTGTTTGCGACCGGTCGTca                                           <  1:5025726/76‑1 (MQ=255)
                                 aTCCCGTTACTCCAAATAAACTCGGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcac                                          <  1:891211/76‑1 (MQ=255)
                                 aTCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcac                                          >  2:2740586/1‑76 (MQ=255)
                                  tCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacc                                         <  1:1555264/76‑1 (MQ=255)
                                    cccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcaccac                                       <  1:3904894/76‑1 (MQ=255)
                                    cccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcaccac                                       >  1:4808835/1‑76 (MQ=255)
                                    cccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcaccac                                       <  2:4540476/76‑1 (MQ=255)
                                     ccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCCCCCCg                                      >  2:5571746/1‑76 (MQ=255)
                                     ccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACg                                      <  1:5616636/76‑1 (MQ=255)
                                       ttACGCCACAGAAGCGCTGAAGCTGCTCACCGC‑‑TCGCGGCATCCACTTTGTGTTTTCGACCGGTCGTCACCACGtt                                    <  1:1692217/76‑1 (MQ=255)
                                       ttACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGtt                                    <  1:2985346/76‑1 (MQ=255)
                                       ttACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGtt                                    <  1:4951167/76‑1 (MQ=255)
                                        tACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACCTTg                                   >  2:5595691/1‑76 (MQ=255)
                                           gCCAAAGAAACTCTGAGGCTGCTCCCCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtg                                <  2:4718824/76‑1 (MQ=255)
                                           gCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtg                                >  2:4310697/1‑76 (MQ=255)
                                            ccAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGCCGTCACCCCGTTGAcgt                               >  2:75312/1‑76 (MQ=255)
                                             cAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgtg                              >  2:3293098/1‑76 (MQ=255)
                                               aaGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                            <  1:4452640/76‑1 (MQ=255)
                                               aaGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                            >  1:5753779/1‑76 (MQ=255)
                                                aGAAACTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTgggg                           <  1:3245886/76‑1 (MQ=255)
                                                  aaaCTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCa                         >  2:5768071/1‑76 (MQ=255)
                                                  aaaCTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCa                         >  1:82223/1‑76 (MQ=255)
                                                   aaCTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaa                        >  2:3608170/1‑76 (MQ=255)
                                                    aCTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                       >  2:2651185/1‑76 (MQ=255)
                                                    aCTCTGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGCTGATCTGGGGCaaa                       >  2:400044/1‑76 (MQ=255)
                                                     ctctGAAGCTGCTCCCCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAt                      <  2:3846632/76‑1 (MQ=255)
                                                     ctctGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAt                      >  1:233913/1‑76 (MQ=255)
                                                     ctctGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAt                      <  2:1332806/76‑1 (MQ=255)
                                                      tctGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtt                     >  2:3585187/1‑76 (MQ=255)
                                                      tctGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtt                     <  2:4281118/76‑1 (MQ=255)
                                                       ctGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTc                    <  1:717836/76‑1 (MQ=255)
                                                       ctGAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTc                    <  2:2188430/76‑1 (MQ=255)
                                                       ctGAACCTCCTCGCCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCCCCACGTTGATGTGGGGCAAATTc                    <  1:417847/76‑1 (MQ=255)
                                                         gAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGTGGCAAATTcgc                  >  2:3528659/1‑76 (MQ=255)
                                                         gAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgc                  >  1:597732/1‑76 (MQ=255)
                                                         gAAGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgc                  <  1:2160434/76‑1 (MQ=255)
                                                         gAACCTGCTCACCGC‑‑GCGCGGCATCAGCTTTGTGTTTGCGACCGGTCGACACCACGTTGATGTGGGGCAAATTcgc                  <  1:4244678/76‑1 (MQ=255)
                                                          aaGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgcg                 >  1:4773220/1‑76 (MQ=255)
                                                           aGCTGCTCCCCGCCGGCGG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGa                <  2:500039/76‑1 (MQ=255)
                                                           aGCTGCTCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGa                >  1:2575441/1‑76 (MQ=255)
                                                            gctgctCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCCAATTCGCGac               >  2:94958/1‑75 (MQ=255)
                                                            gctgctCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGAt               >  2:3546002/1‑76 (MQ=255)
                                                              tgctCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATaa             >  2:3883095/1‑76 (MQ=255)
                                                              tgctCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTACTCACCACGGTGATGAGGGGCAAATTCGCGATaa             >  1:2944984/1‑76 (MQ=255)
                                                               gctCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAAt            >  2:2219302/1‑76 (MQ=255)
                                                               gctCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAAt            >  1:3615301/1‑76 (MQ=255)
                                                                ctCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTgggg                           >  1:5499438/1‑60 (MQ=255)
                                                                ctCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATc           >  1:281797/1‑76 (MQ=255)
                                                                ctCACCGC‑‑GCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATc           >  1:2123268/1‑76 (MQ=255)
                                                                   cgcgc‑‑gcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGgag      <  1:4628938/74‑1 (MQ=255)
                                                                    gcgc‑‑gcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGgaga     >  2:3267445/2‑76 (MQ=255)
                                                                     cgc‑‑gcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAt    <  1:5768072/76‑1 (MQ=255)
                                                                      gc‑‑gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGCGAtt   >  2:4207108/1‑76 (MQ=255)
                                                                      gc‑‑gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAtt   <  2:2575441/76‑1 (MQ=255)
                                                                      gc‑‑gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAtt   <  1:3207490/76‑1 (MQ=255)
                                                                      gc‑‑gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCAACACGTTGATGTGCGGCAAATTCGCGATAATCTGGAGAt  >  2:4610795/1‑75 (MQ=255)
                                                                       c‑‑gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa  >  2:5489158/1‑76 (MQ=255)
                                                                       c‑‑gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa  >  1:2542497/1‑76 (MQ=255)
                                                                              ||                                                                       
AGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGC‑‑GCGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTA  >  NC_000913/4007832‑4007978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: