Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,376,742 G→C 6.6% intergenic (+60/‑145) 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,7420GC 6.6% 218.0 / 4.6 71intergenic (+60/‑145)PP_2083/PP_2084lipase/ribonuclease activity regulator
Reads supporting (aligned to +/- strand):  ref base G (26/40);  new base C (3/2);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCTGGTGGAGCGCTGGCTGCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTTGTGACTGACAGGAGCTT  >  NC_002947/2376598‑2376883
                                                                                                                                                |                                                                                                                                              
gCCGCTGGTGGAGCGCTGGCTGCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTc                                                                                                                                           >  1:238896/1‑150 (MQ=255)
gCCGCTGGTGGAGCGCTGGCTGCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTc                                                                                                                                           >  2:1280404/1‑150 (MQ=255)
     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)
                            ccGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCCCGTTCCtgca                                                                                                                                      >  2:1061946/1‑125 (MQ=255)
                            ccGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCCCGTTCCtgca                                                                                                                                      <  1:1061946/127‑3 (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                                                                                                                          <  1:360821/126‑4 (MQ=39)
                                        gTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCCCGTTCCTGCACAGGAACgcgc                                                                                                                          >  2:360821/1‑123 (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                                                                                                      <  2:650863/128‑1 (MQ=255)
                                                          gTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCgaga                                                                                                      >  1:650863/1‑128 (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                                                                                                                            <  1:1012754/73‑2 (MQ=25)
                                                                                           tGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGCACAGGAACgc                                                                                                                            >  2:1012754/1‑72 (MQ=25)
                                                                                            gTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAAt                                                                                      >  2:1112065/1‑110 (MQ=255)
                                                                                            gTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAAt                                                                                      <  1:1112065/110‑1 (MQ=255)
                                                                                            gTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGgtg                                              >  1:801900/1‑150 (MQ=255)
                                                                                                cTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAAc                                                                                                                              <  1:1580402/66‑1 (MQ=255)
                                                                                                cTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAAc                                                                                                                              >  2:1580402/1‑66 (MQ=255)
                                                                                                cTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGGCAGGATGAATACAAATGACTGCCCGTTCCCGGATGACTACGACAGCCTCCCAGGTGTCGc                                          >  2:1163362/1‑150 (MQ=255)
                                                                                                             ttCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCga                                                                                                        <  1:1643363/75‑1 (MQ=255)
                                                                                                             ttCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCga                                                                                                        >  2:1643363/1‑75 (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                                                                                    >  1:1149184/1‑83 (MQ=255)
                                                                                                                         cccGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGa                                                                                    <  2:1149184/83‑1 (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)
                                                                                                                                                |                                                                                                                                              
GCCGCTGGTGGAGCGCTGGCTGCACAACCCGCTGCTGCCGGTGCATGCGTCCACTGTCGTTGCGGAACCTGTGCCTGCCAGCTGAGTTTGCTGTACCTGTCATGGCCTCTTCGCGGGTGAACCCGCTCCTACACTGGATAACGCGCGTTCCTGTGCAGGAACGGGTCATGCATGGAAGGGCCGAGACAGGATGAATACAAATGACTGCCCGGTCCCGGATGACTACGACAGCCTCCAAGGTGTAGCCTTGGCCGATACCCGCTTTCGTTGTGACTGACAGGAGCTT  >  NC_002947/2376598‑2376883

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: