Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,376,751 T→C 9.8% intergenic (+69/‑136) PP_2083 → / → PP_2084 lipase/ribonuclease activity regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,376,7510TC 9.8% 187.0 / 8.0 70intergenic (+69/‑136)PP_2083/PP_2084lipase/ribonuclease activity regulator
Reads supporting (aligned to +/- strand):  ref base T (25/37);  new base C (4/3);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

TGGTGGAGCGCTGGCTGCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTTGTGACTGACAGGAGCTTCCCATGCAGCATTA  >  NC_002947/2376603‑2376897
                                                                                                                                                    |                                                                                                                                                   
tggtggAGCGCTGGCTGCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCtgtg                                                                                                                                                    <  2:525439/150‑1 (MQ=255)
       gcgcTGGCTGCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAAc                                                                                                                                            <  1:1555461/150‑1 (MQ=255)
       gaactggctgCACAACCCGCTGCTGCAGGTACATGCGTCCACTGTCGTGGCGGAACCTGTGTCTGCCAGCTTAGTTTGCTGTACCTGACACGGCCTCTTCGAGGGCAAACCCGCTCCTACACTGGAGAACGCGGGGTCCTGTGCAGGAAc                                                                                                                                            >  2:365783/4‑150 (MQ=255)
         gctggctgCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACgg                                                                                                                                          <  1:1608492/150‑1 (MQ=255)
         gctggctgCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACgg                                                                                                                                          <  1:1389596/150‑1 (MQ=255)
         gctggctgCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACgg                                                                                                                                          <  1:794970/150‑1 (MQ=255)
          ctggctgCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGTAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACggg                                                                                                                                         <  2:949245/150‑1 (MQ=255)
          ctggctgCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACggg                                                                                                                                         <  2:949250/150‑1 (MQ=255)
          ctggctgCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACggg                                                                                                                                         <  2:1243922/150‑1 (MQ=255)
            ggctgCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAGCCCGCCACTAAACTCGATACCGCGCGTTCCTGTGCAGGAACGGGTc                                                                                                                                       <  1:450532/150‑1 (MQ=255)
                  acaACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTcatgcat                                                                                                                                 <  2:395738/150‑1 (MQ=255)
                       ccGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAg                                                                                                                            >  2:664779/1‑150 (MQ=255)
                       ccGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAg                                                                                                                            <  2:266531/150‑1 (MQ=255)
                          ctgctgCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGc                                                                                                                         <  2:1187613/150‑1 (MQ=255)
                               gcCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCgaga                                                                                                                    >  1:424636/1‑150 (MQ=255)
                                 cGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGCGTTTGCGTTACCTGTCATGGCCTCTTCGCGGGTGACCCCGCTCCTACACTGGATAACGCGTGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACa                                                                                                                  <  1:919382/150‑1 (MQ=255)
                                  ggTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAg                                                                                                                 <  1:965585/150‑1 (MQ=255)
                                   gTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCCCGTTCCTGCACAGGAACgcgc                                                                                                                                        >  2:360821/1‑123 (MQ=39)
                                   gTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCCCGTTCCTGCACAGGAACgcgc                                                                                                                                        <  1:360821/126‑4 (MQ=39)
                                         gCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAAt                                                                                                          <  1:799233/150‑1 (MQ=255)
                                           gTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATAc                                                                                                        <  2:238896/150‑1 (MQ=255)
                                            tCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCGCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACa                                                                                                       >  2:918024/1‑150 (MQ=255)
                                                 tGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGCTAAGGCGCGTTACTGTGCAGGAACGGGTCATGCATGGCAGGGCCGACACAGGCTGAATACAAATgg                                                                                                  >  2:259387/1‑149 (MQ=255)
                                                  gTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGAc                                                                                                 <  2:353131/150‑1 (MQ=255)
                                                     gTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCgaga                                                                                                                    >  1:650863/1‑128 (MQ=255)
                                                     gTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCgaga                                                                                                                    <  2:650863/128‑1 (MQ=255)
                                                        gCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGACTATTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCg                                                                                           <  1:848500/150‑1 (MQ=255)
                                                          ggAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTcat                                                                                                                                     >  1:1523403/1‑106 (MQ=255)
                                                          ggAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTcat                                                                                                                                     <  2:1523403/106‑1 (MQ=255)
                                                           gAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGc                                                                                              <  1:328188/144‑1 (MQ=255)
                                                           gAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGc                                                                                              >  2:328188/1‑144 (MQ=255)
                                                             aCCTGTGCCTGCCAGCTGAGTTTGCTGTACCTAACAGGGCCTGCTCGGGGGTGAACCCGTTACTCCCCTCGGTAGCGCCCGGTCTGCGCCGCCCACCGGTCCTGCATGGAAGGCCCGAGACAGGATGAATACCAACCACTGCCCGGTccc                                                                                      >  1:1606824/1‑150 (MQ=255)
                                                                tgtgCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGa                                                                                   <  2:68779/150‑1 (MQ=255)
                                                                tgtgCCTGCCAGCTGAGTTTGCTGTACCTGCCATGGGCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGa                                                                                   <  1:259387/150‑1 (MQ=255)
                                                                 gtgCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGAt                                                                                  <  1:1280404/150‑1 (MQ=255)
                                                                  tgcctgccAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCtgtg                                                                                                                                                    >  2:1000129/1‑84 (MQ=255)
                                                                  tgcctgccAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCtgtg                                                                                                                                                    <  1:1000129/84‑1 (MQ=255)
                                                                       gccAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTAc                                                                            <  1:664779/150‑1 (MQ=255)
                                                                            cTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAg                                                                       >  1:260609/1‑150 (MQ=255)
                                                                             tGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGc                                                                      <  2:1279983/150‑1 (MQ=255)
                                                                                      tGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGCACAGGAACgc                                                                                                                                          >  2:1012754/1‑72 (MQ=25)
                                                                                      tGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGCACAGGAACgc                                                                                                                                          <  1:1012754/73‑2 (MQ=25)
                                                                                       gTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAAt                                                                                                    >  2:1112065/1‑110 (MQ=255)
                                                                                       gTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAAt                                                                                                    <  1:1112065/110‑1 (MQ=255)
                                                                                       gTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGgtg                                                            >  1:801900/1‑150 (MQ=255)
                                                                                           cTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAAc                                                                                                                                            >  2:1580402/1‑66 (MQ=255)
                                                                                           cTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAAc                                                                                                                                            <  1:1580402/66‑1 (MQ=255)
                                                                                           cTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGGCAGGATGAATACAAATGACTGCCCGTTCCCGGATGACTACGACAGCCTCCCAGGTGTCGc                                                        >  2:1163362/1‑150 (MQ=255)
                                                                                                        ttCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCga                                                                                                                      >  2:1643363/1‑75 (MQ=255)
                                                                                                        ttCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCga                                                                                                                      <  1:1643363/75‑1 (MQ=255)
                                                                                                           gcgGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGc                                        >  1:601213/1‑150 (MQ=255)
                                                                                                                  aaCGCGGTCCTACCCGGGAGAACGCGCGGTGCTGGGCAGGGACGGGGCATGCAGGGAAGGGCCGAGACAGGATGAAGAGAAATGACGGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGtt                                 <  1:540983/150‑1 (MQ=255)
                                                                                                                  aaCCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGtt                                 <  2:1606824/150‑1 (MQ=255)
                                                                                                                   aCCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTc                                    <  2:500320/146‑1 (MQ=255)
                                                                                                                   aCCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTc                                    >  1:500320/1‑146 (MQ=255)
                                                                                                                    cccGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGa                                                                                                  <  2:1149184/83‑1 (MQ=255)
                                                                                                                    cccGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGa                                                                                                  >  1:1149184/1‑83 (MQ=255)
                                                                                                                    cccGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTtgt                               <  1:1483719/150‑1 (MQ=255)
                                                                                                                    cccGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTtgt                               >  2:1060545/1‑150 (MQ=255)
                                                                                                                      cGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTTGtga                             <  1:1163362/150‑1 (MQ=255)
                                                                                                                        cTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTTGtgact                           >  2:33795/1‑150 (MQ=255)
                                                                                                                         tCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGcc                                                                     >  1:547546/1‑107 (MQ=255)
                                                                                                                         tCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGcc                                                                     <  2:547546/107‑1 (MQ=255)
                                                                                                                           cTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGgtgt                                                           <  1:536071/115‑1 (MQ=255)
                                                                                                                           cTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGgtgt                                                           >  2:536071/1‑115 (MQ=255)
                                                                                                                               acTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTTGTGACTGACAGGa                    <  2:1410002/150‑1 (MQ=255)
                                                                                                                                   gATAACGAGCGTTACTGTGCAGGAACGGGTCATGCATGGAAGGGGGGAGAGCAGATGAATACAAATGACTGCCCGGTGCCGGGTGACTACGACAGACACTAAGGCGTAGCCTTGGCCGATACCGACTTTGGTTGTGAGTGACAGGAGCtt                >  2:1464066/1‑150 (MQ=255)
                                                                                                                                          cccgTTCCTGCACAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTg                                                    <  1:1429343/105‑1 (MQ=255)
                                                                                                                                          cccgTTCCTGCACAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTg                                                    >  2:1429343/3‑107 (MQ=255)
                                                                                                                                              ttCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTTGTGACTGACAGGAGCTTCCCATgcagca     >  2:397924/1‑150 (MQ=255)
                                                                                                                                                ccTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGa                                                                                   >  2:1578779/1‑70 (MQ=255)
                                                                                                                                                ccTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGa                                                                                   <  1:1578779/70‑1 (MQ=255)
                                                                                                                                                 cTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCGCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTTGTGACTGACAGGAGCTTCCCATGCAGCATTa  >  1:921662/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGTGGAGCGCTGGCTGCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTTGTGACTGACAGGAGCTTCCCATGCAGCATTA  >  NC_002947/2376603‑2376897

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