Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,641,506 C→G 6.4% S373C (TCC→TGC)  PP_4951 → EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,641,5060CG 6.4% 376.9 / 13.4 119S373C (TCC→TGC) PP_4951EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (60/50);  new base G (3/5);  total (64/55)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 122 total aligned reads displayed.

ACCGCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTAC  >  NC_002947/5641357‑5641653
                                                                                                                                                     |                                                                                                                                                   
   gCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAt                                                                                                                                                  >  2:2099405/1‑150 (MQ=255)
    cTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATg                                                                                                                                                 >  2:1001191/1‑150 (MQ=255)
       aTCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAg                                                                                                                                              <  2:2485634/150‑1 (MQ=255)
       aTCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAg                                                                                                                                              <  1:277598/150‑1 (MQ=255)
        tCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTCCTCCAACTGCTGGTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGt                                                                                                                                             >  2:3049537/1‑150 (MQ=255)
         cGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGtt                                                                                                                                            <  1:2996827/150‑1 (MQ=255)
             ctgTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACt                                                                                                                                        <  1:619797/148‑1 (MQ=255)
             cGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCATGGAGTTCACt                                                                                                                                        <  2:2001794/150‑1 (MQ=255)
             cGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCATGGAGTTCACt                                                                                                                                        <  1:1644446/150‑1 (MQ=255)
             cGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACt                                                                                                                                        <  2:2709457/150‑1 (MQ=255)
             cGGTTCGGCTATCGCATCATCCTCACCGGTAATACATTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCATGGAGTTCACt                                                                                                                                        <  1:939804/150‑1 (MQ=255)
                ttCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCg                                                                                                                                     >  2:1814741/1‑150 (MQ=255)
                ttCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCg                                                                                                                                     <  1:1814741/150‑1 (MQ=255)
                 tCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGa                                                                                                                                    <  2:821963/150‑1 (MQ=255)
                 tCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGa                                                                                                                                    <  1:458538/150‑1 (MQ=255)
                 tCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATACTGCTGGCAAGCCTGGTCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGa                                                                                                                                    <  2:872036/150‑1 (MQ=255)
                  cGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGa                                                                                                                                    <  2:1746110/149‑1 (MQ=255)
                  cGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGa                                                                                                                                    >  1:1746110/1‑149 (MQ=255)
                   ggCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATg                                                                                                                                  <  2:1093531/150‑1 (MQ=255)
                   ggCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATg                                                                                                                                  <  1:854970/150‑1 (MQ=255)
                   ggCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATg                                                                                                                                  <  1:667551/150‑1 (MQ=255)
                         cGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAAcaca                                                                                                                            <  2:2727671/150‑1 (MQ=255)
                            atcatcCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTg                                                                                                                         >  1:169998/1‑150 (MQ=255)
                             tcatcCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGa                                                                                                                        <  1:740375/150‑1 (MQ=255)
                             tcatcCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGa                                                                                                                        <  1:1795048/150‑1 (MQ=255)
                             tcatcCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGa                                                                                                                        <  1:2779905/150‑1 (MQ=255)
                              gatcCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGAc                                                                                                                       >  2:2108446/2‑150 (MQ=255)
                                 ccTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCt                                                                                                                    >  1:2960042/1‑150 (MQ=255)
                                   tCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCa                                                                                                                  >  2:2591497/1‑150 (MQ=255)
                                    cACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCAt                                                                                                                 <  1:1792794/150‑1 (MQ=255)
                                     aCCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATc                                                                                                                >  1:2276496/1‑150 (MQ=255)
                                          tAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGAtct                                                                                                           <  2:2276496/150‑1 (MQ=255)
                                           aaTACCTTGCTGTTGGGCGTTCTGCTGGCGACCGGTGTCCTGGCGGACGAGCAGAACCCTTATTGACTTGGCGGCGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGAtctc                                                                                                          <  1:2361850/150‑1 (MQ=255)
                                             tACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTcga                                                                                                        <  1:1176600/150‑1 (MQ=255)
                                             tACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTcga                                                                                                        <  2:1749747/150‑1 (MQ=255)
                                             tACCTTGCTGTTGGGCAGCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGTCCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTcga                                                                                                        <  1:1459673/150‑1 (MQ=255)
                                              aCCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTTACTGCGATGAACACAGTGACGCTCATCGATCTcgac                                                                                                       >  1:805731/1‑150 (MQ=255)
                                              aCCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTcgac                                                                                                       <  1:232006/150‑1 (MQ=255)
                                                cTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCAc                                                                                                                                         <  1:2079822/114‑1 (MQ=255)
                                                cTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAg                                                                                                                           <  2:1610147/128‑1 (MQ=255)
                                                     tGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGACTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgcca                                                                                                >  2:1590493/1‑150 (MQ=255)
                                                     tGTTGGGCATATTGCTGGCAAGCCTGGGCCTGTTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgcca                                                                                                <  2:2960042/150‑1 (MQ=255)
                                                      gTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGa                                                                                                                                 >  2:2153537/1‑116 (MQ=255)
                                                      gTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGa                                                                                                                                 <  1:2153537/116‑1 (MQ=255)
                                                      gTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgccag                                                                                               >  1:1043871/1‑150 (MQ=255)
                                                      gTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGCCAGAAGGGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgccag                                                                                               >  2:2397690/1‑150 (MQ=255)
                                                       ttGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGAtc                                                                                                            <  2:2890009/136‑1 (MQ=255)
                                                       ttGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGAtc                                                                                                            >  1:2890009/1‑136 (MQ=255)
                                                       ttGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgccagc                                                                                              >  1:101053/1‑150 (MQ=255)
                                                             aTCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTACATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACgccagcgccagc                                                                                        >  2:2716910/1‑150 (MQ=255)
                                                                 tgctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCg                                                                                    >  2:190128/1‑150 (MQ=255)
                                                                 tgctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCg                                                                                    <  2:2771987/150‑1 (MQ=255)
                                                                      gCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAAc                                                                               >  1:2821938/1‑150 (MQ=255)
                                                                         aGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAaca                                                                                                                              >  2:135403/1‑100 (MQ=255)
                                                                         aGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGGGTTAGCAGCGGTAACAGc                                                                            <  2:2432073/150‑1 (MQ=255)
                                                                          gcctgggcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGcc                                                                           >  2:2933890/1‑150 (MQ=255)
                                                                            ctgggcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCctg                                                                         >  1:454491/1‑150 (MQ=255)
                                                                                gcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGc                                                                                                                                      >  1:1475215/1‑85 (MQ=255)
                                                                                gcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGc                                                                                                                                      <  2:1475215/85‑1 (MQ=255)
                                                                                gcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGt                                                                     >  1:2816801/1‑150 (MQ=255)
                                                                                   tggtggACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGACTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCgg                                                                  <  1:1590493/150‑1 (MQ=255)
                                                                                        gACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCATGGAGTTCACt                                                                                                                                        <  2:260112/75‑1 (MQ=38)
                                                                                        gACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCATGGAGTTCACt                                                                                                                                        >  1:260112/1‑75 (MQ=38)
                                                                                           gAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTcgcg                                                          >  1:2169291/1‑150 (MQ=255)
                                                                                           gAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTcgcg                                                          >  1:2734124/1‑150 (MQ=255)
                                                                                             gCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAAATGCCTACAGTTCACAAATATTTAGACAGGTACTCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCa                                                        >  1:1908659/1‑150 (MQ=255)
                                                                                              cAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTc                                                                                                                                                     >  2:308957/1‑56 (MQ=255)
                                                                                                 aCCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTg                                                    >  2:2448326/1‑150 (MQ=255)
                                                                                                 aCCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTg                                                    >  2:815500/1‑150 (MQ=255)
                                                                                                 aCCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTg                                                    >  1:1359888/1‑150 (MQ=255)
                                                                                                   cccTTATTGGCTGCTGTTGGTGCAGCTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTc                                                  >  2:2026751/1‑150 (MQ=255)
                                                                                                        aTTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGa                                             <  2:2821938/150‑1 (MQ=255)
                                                                                                            gctgctGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCt                                         <  1:2933890/150‑1 (MQ=255)
                                                                                                                ctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGgt                                     <  1:952828/150‑1 (MQ=255)
                                                                                                                ctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGgt                                     <  2:2169291/150‑1 (MQ=255)
                                                                                                                 tgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGgtg                                    >  2:567949/1‑150 (MQ=255)
                                                                                                                 tgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGgtg                                    >  1:2727285/1‑150 (MQ=255)
                                                                                                                 tgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGgtg                                    <  1:1838563/150‑1 (MQ=255)
                                                                                                                     ggTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGc                                <  2:2535175/150‑1 (MQ=255)
                                                                                                                      gTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGtgtg                                  <  1:2243244/147‑1 (MQ=255)
                                                                                                                      gTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGcc                               >  2:1690529/1‑150 (MQ=255)
                                                                                                                      gTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGcc                               >  1:607058/1‑150 (MQ=255)
                                                                                                                      gTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGcc                               >  1:2471899/1‑150 (MQ=255)
                                                                                                                       tGCAACTTGGCCTTCTGGGCGCAGTGAACTCGCTGCAGTTCACTGCGATGAACACAGTGAAGCGCATCAAGCTCTAGGACGCCATCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGCCCt                              >  2:1884928/1‑150 (MQ=255)
                                                                                                                       tGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCATTGGTAATAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCt                              <  1:376383/150‑1 (MQ=255)
                                                                                                                       tGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCCCAACCATCGCTGAACCTGGGTGTGGCCt                              >  2:230323/1‑150 (MQ=255)
                                                                                                                        gCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTg                             <  1:1001191/150‑1 (MQ=255)
                                                                                                                        gCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTg                             >  1:1443962/1‑150 (MQ=255)
                                                                                                                            cTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGcc                         >  1:2908072/1‑150 (MQ=255)
                                                                                                                              tGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCgg                       >  1:1015991/1‑150 (MQ=255)
                                                                                                                              tGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCgg                       <  1:3049537/150‑1 (MQ=255)
                                                                                                                                 ccTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGgcgc                    >  2:2902955/1‑150 (MQ=255)
                                                                                                                                  cTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCa                   >  2:2695198/1‑150 (MQ=255)
                                                                                                                                          cgcAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGTACTGTTggg           >  2:1993160/1‑150 (MQ=255)
                                                                                                                                            cAGTGAACTGCATGGAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGa                                             >  2:1678649/1‑114 (MQ=255)
                                                                                                                                            cAGTGAACTGCATGGAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGa                                             <  1:1678649/114‑1 (MQ=255)
                                                                                                                                               tGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCt      >  2:1831000/1‑150 (MQ=255)
                                                                                                                                                 aaCTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCttt    >  1:3053473/1‑150 (MQ=255)
                                                                                                                                                   cTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGGGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTAc  >  2:422132/1‑150 (MQ=255)
                                                                                                                                                   cTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTCTCGGTGGTCGCGCAACTGTCACTGAGACTAGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTAc  >  2:793598/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACCGCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTAC  >  NC_002947/5641357‑5641653

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