Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,829,283 C→G E2301Q (GAA→CAA)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,829,2830CG99.7% 1357.2 / ‑5.0 385E2301Q (GAA→CAA) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (186/198);  total (187/198)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 385 total aligned reads displayed.

GCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCGGTACATAGGCGC  >  NC_002947/4829136‑4829429
                                                                                                                                                   |                                                                                                                                                  
ccATGGCGAGTTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCa                                                                                                                                                  >  2:6870042/2‑150 (MQ=255)
 cATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAg                                                                                                                                                 >  1:5363525/1‑150 (MQ=255)
    ggCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCt                                                                                                                                              >  2:1850550/1‑150 (MQ=255)
       gAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCa                                                                                                                                           >  1:7914728/1‑150 (MQ=255)
        aaCTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAg                                                                                                                                          <  1:5151233/150‑1 (MQ=255)
          cTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCt                                                                                                                                        >  2:7059291/1‑150 (MQ=255)
           tCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTg                                                                                                                                       >  2:8015540/1‑150 (MQ=255)
           tCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTg                                                                                                                                       >  1:8618827/1‑150 (MQ=255)
             gCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGc                                                                                                                                     <  1:5340260/150‑1 (MQ=255)
               cATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCa                                                                                                                                   <  2:4007959/150‑1 (MQ=255)
                aTGGCTGACCACCCCGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGGAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCAc                                                                                                                                  <  2:7741580/150‑1 (MQ=255)
                     tGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCtt                                                                                                                             >  2:4370242/1‑150 (MQ=255)
                     tGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCtt                                                                                                                             <  2:4461321/150‑1 (MQ=255)
                       accaccACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCa                                                                                                                           <  1:8049266/150‑1 (MQ=255)
                       accaccACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCa                                                                                                                           <  2:6770341/150‑1 (MQ=255)
                       accaccACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCa                                                                                                                           <  1:5061182/150‑1 (MQ=255)
                       accaccACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCa                                                                                                                           >  1:872047/1‑150 (MQ=255)
                       accaccACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCa                                                                                                                           <  1:1004366/150‑1 (MQ=255)
                         caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAgg                                                                                                                         >  1:7137154/1‑150 (MQ=255)
                         caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAgg                                                                                                                         >  1:410453/1‑150 (MQ=255)
                          accacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGc                                                                                                                        >  2:8370858/1‑150 (MQ=255)
                                  cTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGcc                                                                                                                <  1:594134/150‑1 (MQ=255)
                                    tGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCa                                                                                                              <  2:7923430/150‑1 (MQ=255)
                                    tGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCa                                                                                                              <  2:5263219/150‑1 (MQ=255)
                                    tGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCa                                                                                                              <  2:6948233/150‑1 (MQ=255)
                                      ggTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAgg                                                                                                            <  1:4435707/150‑1 (MQ=255)
                                      ggTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAgg                                                                                                            <  2:4303167/150‑1 (MQ=255)
                                      ggTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAgg                                                                                                            <  2:4855572/150‑1 (MQ=255)
                                      ggTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAgg                                                                                                            <  1:852432/150‑1 (MQ=255)
                                        tttGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGgt                                                                                                          <  1:218418/150‑1 (MQ=255)
                                        tttGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCACTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGgt                                                                                                          <  2:799807/150‑1 (MQ=255)
                                          tGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGgtgt                                                                                                        <  1:1234781/150‑1 (MQ=255)
                                           gCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGtgtg                                                                                                       >  1:1777232/1‑150 (MQ=255)
                                               gTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTgcagc                                                                                                   >  2:2183002/1‑150 (MQ=255)
                                                   aGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAgcg                                                                                               <  1:7370213/150‑1 (MQ=255)
                                                    gCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAgcgc                                                                                              <  2:1974450/150‑1 (MQ=255)
                                                     ccTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCa                                                                                             >  2:6634885/1‑150 (MQ=255)
                                                     ccTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCa                                                                                             <  1:8483684/150‑1 (MQ=255)
                                                       tGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACg                                                                                           >  2:5444408/1‑150 (MQ=255)
                                                         aGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGGCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgct                                                                                         <  2:6350507/150‑1 (MQ=255)
                                                          ggTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctg                                                                                        >  1:304077/1‑150 (MQ=255)
                                                            tgtAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGgt                                                                                                          <  1:6754943/130‑1 (MQ=255)
                                                            tgtAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctggc                                                                                      <  1:3481264/150‑1 (MQ=255)
                                                            ggtAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctggc                                                                                      <  2:1261513/149‑1 (MQ=255)
                                                             gtAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctggct                                                                                     <  1:1992645/150‑1 (MQ=255)
                                                              tAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctggctg                                                                                    >  2:8718867/1‑150 (MQ=255)
                                                                  tCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctg                                                                                        <  2:5036066/142‑1 (MQ=255)
                                                                  tCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGagca                                                                                >  1:292787/1‑150 (MQ=255)
                                                                  tCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGCGTGCGCTGCTGCACACcctggg                                                                                      >  1:387148/1‑143 (MQ=255)
                                                                  tCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAACTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGagca                                                                                >  1:2479961/1‑150 (MQ=255)
                                                                    aGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGagcagc                                                                              <  1:376449/150‑1 (MQ=255)
                                                                      gTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCaa                                                                            >  1:7721165/1‑150 (MQ=255)
                                                                        aaGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCAGCCAGCGTGCGCTGCTGCACACCCTGGGCGagctgcctga                                                                          <  2:1528807/150‑5 (MQ=25)
                                                                          gCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCg                                                                        >  2:7324340/1‑150 (MQ=255)
                                                                           ccTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGa                                                                       >  1:3817647/1‑150 (MQ=255)
                                                                             tGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATg                                                                     <  2:3439938/150‑1 (MQ=255)
                                                                              ggTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGc                                                                    <  2:5691640/150‑1 (MQ=255)
                                                                                tgctgcGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCgg                                                                  <  1:1459600/150‑1 (MQ=255)
                                                                                  ctgcGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctggc                                                                                      <  1:4950285/128‑1 (MQ=255)
                                                                                  ctgcGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCGGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCt                                                                <  1:480442/150‑1 (MQ=255)
                                                                                    gcGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGt                                                              <  2:2739677/150‑1 (MQ=255)
                                                                                        cAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTc                                                          <  1:5816366/150‑1 (MQ=255)
                                                                                         aGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCg                                                         >  2:4921451/1‑150 (MQ=255)
                                                                                           gCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCTCACGATGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGAt                                                       <  1:6664270/150‑1 (MQ=255)
                                                                                           gCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGCCTTCGAt                                                       >  2:7539124/1‑150 (MQ=255)
                                                                                             tGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCa                                                     >  1:8314930/1‑150 (MQ=255)
                                                                                             tGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCa                                                     <  2:3805503/150‑1 (MQ=255)
                                                                                             tGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCa                                                     >  1:6312202/1‑150 (MQ=255)
                                                                                              gAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCAt                                                    >  2:1464448/1‑150 (MQ=255)
                                                                                               aGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATg                                                   >  1:6248129/1‑150 (MQ=255)
                                                                                                  ttGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGAt                                                       >  1:6876934/1‑143 (MQ=255)
                                                                                                    ctacAGCGATGCACCGCTGGAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAgcg                                              >  2:3149212/4‑150 (MQ=255)
                                                                                                     tCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAgcgc                                             <  2:6712467/150‑1 (MQ=255)
                                                                                                        aGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctggc                                                                                      <  2:8593718/106‑1 (MQ=255)
                                                                                                        aGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAgg                                          >  2:8265538/1‑150 (MQ=255)
                                                                                                         gCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGc                                         >  1:2617156/1‑150 (MQ=255)
                                                                                                             tGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTc                                     <  1:6119488/150‑1 (MQ=255)
                                                                                                             tGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGACTGTCTTCGATCATGTAGCGCAGGCGGTc                                     <  1:6223016/150‑1 (MQ=255)
                                                                                                               taTCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGa                                   <  2:8517320/149‑1 (MQ=255)
                                                                                                                 tCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAgg                                                                                                            >  2:2702170/1‑75 (MQ=255)
                                                                                                                 tCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCa                                            <  1:7837491/139‑1 (MQ=255)
                                                                                                                 tCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGc                                 >  1:3372555/1‑150 (MQ=255)
                                                                                                                 tCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGc                                 >  2:6421778/1‑150 (MQ=255)
                                                                                                                 tCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCAAGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGc                                 <  2:2821933/150‑1 (MQ=255)
                                                                                                                   gCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCgg                               <  1:3639191/150‑1 (MQ=255)
                                                                                                                       tAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTAc                           >  1:1112157/1‑150 (MQ=255)
                                                                                                                       tAGGCCGCCAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTAc                           <  1:6510728/150‑1 (MQ=255)
                                                                                                                            cgccAGGCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCggg                      <  1:928393/150‑1 (MQ=255)
                                                                                                                              ccAGCCTGGCGGCCTCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGagca                                                                                <  2:3320652/90‑1 (MQ=255)
                                                                                                                                 ggCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGa                 <  2:3818600/150‑1 (MQ=255)
                                                                                                                                 ggCTGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGa                 <  2:7835696/150‑1 (MQ=255)
                                                                                                                                    ttgcggcATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCg              >  1:7586823/3‑150 (MQ=255)
                                                                                                                                    tGGCGGCATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCg              <  2:1737802/150‑1 (MQ=255)
                                                                                                                                      gcggcCTCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCGGt            <  2:5613449/150‑1 (MQ=255)
                                                                                                                                         gcATCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGCGCGGGTACTCGGGGTCGAGCGGTAcc         >  2:3547635/1‑149 (MQ=255)
                                                                                                                                           aTCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCGGTACATa       <  2:1562682/150‑1 (MQ=255)
                                                                                                                                           aTCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCGGTACATa       >  1:2702456/1‑150 (MQ=255)
                                                                                                                                            tCGTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCGGTACATAg      <  1:1540500/150‑1 (MQ=255)
                                                                                                                                              gTCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCGGTACATAGgc    >  2:2551774/1‑150 (MQ=255)
                                                                                                                                               tCTTGCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCGGTACATAGgcg   >  1:4543457/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGAGCGGTACATAGGCGC  >  NC_002947/4829136‑4829429

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