Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,777,693 T→C 5.2% Q540R (CAG→CGG)  pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,777,6930TC 5.2% 334.7 / 4.1 109Q540R (CAG→CGG) pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (48/55);  new base C (3/3);  total (51/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GGGTAGTCCGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCTGATTCCACTGG  >  NC_002947/4777549‑4777842
                                                                                                                                                |                                                                                                                                                     
gggTAGTCCGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTg                                                                                                                                                  <  1:2086420/150‑1 (MQ=255)
  gTAGTCCGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGc                                                                                                                                                <  1:2083951/150‑1 (MQ=255)
  gTAGTCCGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGc                                                                                                                                                <  2:635695/150‑1 (MQ=255)
     gTCCGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTg                                                                                                                                             <  1:850062/150‑1 (MQ=255)
     gTCCGGGTCGAGCGGCACGTAAGCGCCCCCTGACAAGAATATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTg                                                                                                                                             >  1:1716691/1‑150 (MQ=255)
      tCCGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTgg                                                                                                                                            >  2:1744173/1‑150 (MQ=255)
      tCCGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACCAGTTGGcccac                                                                                                                                            <  1:464368/150‑4 (MQ=255)
       ccGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCAACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGc                                                                                                                                           <  2:1690804/150‑1 (MQ=255)
        cGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCg                                                                                                                                          >  1:1819461/1‑150 (MQ=255)
          ggTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGt                                                                                                                                        <  1:1675253/150‑1 (MQ=255)
           gTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGtt                                                                                                                                                   >  1:1100119/1‑138 (MQ=255)
           gTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGtt                                                                                                                                                   <  2:1100119/138‑1 (MQ=255)
           gTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGtt                                                                                                                                       >  2:1717872/1‑150 (MQ=255)
           gTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGtt                                                                                                                                       >  1:2188925/1‑150 (MQ=255)
             cGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGa                                                                                                                                     <  2:1830955/150‑1 (MQ=255)
                  ggCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGAAGGATTTCCGGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCg                                                                                                                                <  2:1716691/150‑1 (MQ=255)
                   gCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCgg                                                                                                                               >  2:718857/1‑150 (MQ=255)
                   gCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCgg                                                                                                                               <  1:718857/150‑1 (MQ=255)
                       gTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTa                                                                                                                           >  2:893380/1‑150 (MQ=255)
                         aaGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGc                                                                                                                         >  1:939338/1‑150 (MQ=255)
                             gcCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCaa                                                                                                                     >  1:1821135/1‑150 (MQ=255)
                              cccccTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAg                                                                                                                    <  2:939338/150‑1 (MQ=255)
                              cccccTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAg                                                                                                                    <  2:2092141/150‑1 (MQ=255)
                                cccTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGt                                                                                                                  >  1:1209221/1‑150 (MQ=255)
                                cccTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGt                                                                                                                  >  1:1463555/1‑150 (MQ=255)
                                cccTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGt                                                                                                                  <  2:1209221/150‑1 (MQ=255)
                                cccTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGAGTGGAACAACCGCATAGCCTGGCTGTCTAGGTCGGCGTCGCTCAAggc                                                                                                                  >  2:22950/1‑149 (MQ=255)
                                 ccTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTc                                                                                                                 >  1:1954144/1‑150 (MQ=255)
                                  cTGCCTTGAAGATCGCCAGCAGGCCTACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCt                                                                                                                <  1:2062752/150‑1 (MQ=255)
                                    gCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGc                                                                                                              <  2:1222993/150‑1 (MQ=255)
                                      cTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTc                                                                                                            <  2:1956497/150‑1 (MQ=255)
                                       ttGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCg                                                                                                           >  2:2275152/1‑150 (MQ=255)
                                            ctcagcaagcagGCCGACCAGCACTCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCaaa                                                                                                      >  2:1441570/8‑150 (MQ=255)
                                             aTCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCaaaa                                                                                                     <  1:2166913/150‑1 (MQ=255)
                                                 ccAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACCAGTTGGccca                                                                                                                                             >  1:1349893/1‑104 (MQ=255)
                                                 ccAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACCAGTTGGccca                                                                                                                                             <  2:1349893/106‑3 (MQ=255)
                                                  cagcagGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTgg                                                                                                                                            <  1:1627749/106‑1 (MQ=255)
                                                  cagcagGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTgg                                                                                                                                            >  2:1627749/1‑106 (MQ=255)
                                                  cagcagGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGTAg                                                                                                <  1:1035431/150‑1 (MQ=255)
                                                  cagcagGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAg                                                                                                <  1:1988935/150‑1 (MQ=255)
                                                       ggCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCa                                                                                           >  2:2310623/1‑150 (MQ=255)
                                                       ggCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTAGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGGAGCGCCa                                                                                           >  1:926723/1‑150 (MQ=255)
                                                        gCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCAc                                                                                          <  1:949902/150‑1 (MQ=255)
                                                          ccaccaGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACgg                                                                                        >  1:1936974/1‑150 (MQ=255)
                                                                gCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGACCCGGCATCgg                                                                                  <  2:2329980/150‑1 (MQ=255)
                                                                  aTGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGc                                                                                >  2:375727/1‑150 (MQ=255)
                                                                   tGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCg                                                                               <  1:568858/150‑1 (MQ=255)
                                                                    gcgcGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGt                                                                              <  2:1819461/150‑1 (MQ=255)
                                                                      gcgGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGc                                                                            >  1:1978627/1‑150 (MQ=255)
                                                                      gcgGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGc                                                                            >  2:1927110/1‑150 (MQ=255)
                                                                         gCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTg                                                                                                                                                  <  2:486223/77‑1 (MQ=255)
                                                                         gCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTg                                                                                                                                                  >  1:486223/1‑77 (MQ=255)
                                                                             gcgcgGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCa                                                                                           <  1:1783725/128‑1 (MQ=255)
                                                                             gcgcgGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGacacac                                                                     <  2:767391/150‑1 (MQ=255)
                                                                                 ggCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGa                                                                          <  2:486265/141‑1 (MQ=255)
                                                                                 ggCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGa                                                                          <  2:486421/141‑1 (MQ=255)
                                                                                 ggCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGa                                                                          >  1:486421/1‑141 (MQ=255)
                                                                                 ggCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGa                                                                          >  1:486265/1‑141 (MQ=255)
                                                                                 ggCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACctgc                                                                 <  1:2275152/150‑1 (MQ=255)
                                                                                     tCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGtt                                                             <  1:2310623/150‑1 (MQ=255)
                                                                                     tCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGtt                                                             <  2:1954144/150‑1 (MQ=255)
                                                                                     tCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGtt                                                             <  1:1371916/150‑1 (MQ=255)
                                                                                      cGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTc                                                            <  1:1717872/150‑1 (MQ=255)
                                                                                        ccACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAg                                                                                                                    <  1:503568/92‑1 (MQ=255)
                                                                                        ccACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAg                                                                                                                    >  2:503568/1‑92 (MQ=255)
                                                                                          acacCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAAc                                                        >  1:698822/1‑150 (MQ=255)
                                                                                          acacCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCCCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGacc                                                        >  2:2113012/1‑148 (MQ=255)
                                                                                           cacccaccAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGCTCGGCTCGCGAAAACGCAGCGCCCCGGCATCGGGCGTGCGACACACCTGCTATTCGAACa                                                       >  2:437785/1‑150 (MQ=255)
                                                                                            acccaccAGCACTTCAGGGCCGACGCCACGCTCGCGCAGTGGGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGTGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAg                                                      <  1:22950/150‑1 (MQ=255)
                                                                                            acccaccAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATc                                                                                    <  2:2228193/120‑1 (MQ=255)
                                                                                            acccaccAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGACGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATc                                                                                    >  1:2228193/1‑120 (MQ=255)
                                                                                             cccaccAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGACAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAgg                                                     <  2:1651526/150‑1 (MQ=255)
                                                                                             cccaccAGCACTTCAGGGCCGACGCCACGCTCGCCTCGGGCATTGGCCAACTGGTTGGCCTGGCGGTTGAGTGAGGGTTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAgg                                                     <  1:1103232/150‑1 (MQ=255)
                                                                                              ccaccaGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTa                                                                                                                           >  1:1353533/1‑79 (MQ=255)
                                                                                              ccaccaGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCCTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTa                                                                                                                           <  2:1353533/79‑1 (MQ=255)
                                                                                                accaGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGa                                                  >  1:853674/1‑150 (MQ=255)
                                                                                                 ccaGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTC‑AGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATa                                                <  2:1019369/150‑1 (MQ=255)
                                                                                                   aGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAg                                               >  2:1458011/1‑150 (MQ=255)
                                                                                                      aCTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACg                                            <  1:1744173/150‑1 (MQ=255)
                                                                                                               ccGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTtg                                   <  2:1559403/150‑1 (MQ=255)
                                                                                                                 gACGCCGCGCCCGCGCAGGGCGTGGGCCAACTGGTTGGCGTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCAAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCtgtgg                                 <  2:1338249/150‑5 (MQ=255)
                                                                                                                     ccACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAATGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCCGCTTGTCGAACAGGCGATAGACGCCGAGCTTGGTGGGa                             <  1:550988/150‑1 (MQ=255)
                                                                                                                     ccACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGa                             <  2:853674/150‑1 (MQ=255)
                                                                                                                      cacaCTCGCGCAGGGCGGGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGCGCTCTGTGGGAt                            >  2:633519/1‑150 (MQ=255)
                                                                                                                        cGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAg                          <  2:525489/150‑1 (MQ=255)
                                                                                                                        cGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAg                          <  2:1936974/150‑1 (MQ=255)
                                                                                                                             gcgcAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCg                                                                                                                                <  2:2135363/43‑1 (MQ=255)
                                                                                                                             gcgcAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCg                                                                                                                                >  1:2135363/1‑43 (MQ=255)
                                                                                                                                     gccaggCCAACCAGTTGGCCCGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATAGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCt             >  1:2216160/5‑150 (MQ=255)
                                                                                                                                       tGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCTGa           <  2:2188925/150‑1 (MQ=255)
                                                                                                                                       caggCCAACCAGTTGGCCCGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATAGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCTGa           <  2:2216160/148‑1 (MQ=255)
                                                                                                                                          gCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCCTCAGGGGTGCGACACACCTGCTGTTCGAACAGGCGCTAGATTCCGATCAGGGTGGGATAGGCCATCTGCGTCTGATTc        >  1:1102551/1‑150 (MQ=255)
                                                                                                                                          gCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCTGATTc        >  1:232958/1‑150 (MQ=255)
                                                                                                                                            cAACTGGTTGGCCTGTCGGTTGAGTTCGGCGTAGCTCAAGGGCTGCGCGCCAAAACGCAGCGCAAAGGCATCGGGCATGCGACACACCTGCTGTTCGAACAGGCGATAAACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCTGATGCCa      >  2:376278/1‑150 (MQ=255)
                                                                                                                                            cAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCGATCTGCGTGTGATTCCa      >  1:577218/1‑150 (MQ=255)
                                                                                                                                            cAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCTGATTCCa      <  1:656311/150‑1 (MQ=255)
                                                                                                                                            cAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCTGATTCCa      >  2:656311/1‑150 (MQ=255)
                                                                                                                                                tGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGAt                                                 <  2:1987053/103‑1 (MQ=255)
                                                                                                                                                tGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGAt                                                 >  1:1987053/1‑103 (MQ=255)
                                                                                                                                                tGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCTGATTCCACTgg  >  2:1586320/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGGTAGTCCGGGTCGAGCGGCACGTAAGCGCCCCCTGCCTTGAAGATCGCCAGCAGGCCCACCAGCATGCGCGGCCCGCGCGGCATCGCCACACCCACCAGCACTTCAGGGCCGACGCCACGCTCGCGCAGGGCGTGGGCCAACTGGTTGGCCTGGCGGTTGAGTTCGGCGTAGCTCAAGGTCTGCTCGCCAAAACGCAGCGCCACGGCATCGGGCGTGCGACACACCTGCTGTTCGAACAGGCGATAGACGCCGAGCTTGGTGGGATAGGCCATCTGCGTCTGATTCCACTGG  >  NC_002947/4777549‑4777842

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: ‑