Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,878,872 G→C 6.0% Q221H (CAG→CAC mexF → multidrug RND transporter MexF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,878,8720GC 6.0% 365.6 / 10.0 114Q221H (CAG→CACmexFmultidrug RND transporter MexF
Reads supporting (aligned to +/- strand):  ref base G (50/57);  new base C (2/5);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 114 total aligned reads displayed.

GGCGTAGGCGACGTGCAGCTGTTTGGCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCGCCTGAAAGACATCGCCCGGG  >  NC_002947/3878726‑3879020
                                                                                                                                                    |                                                                                                                                                    
ggcgTAGGCGACGTGCAGCTGTTTGGCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTg                                                                                                                                                   >  1:1019250/1‑150 (MQ=255)
 gcgTAGGCGACGTGCAGCTGTTTGGCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGTCGGCCAGTGATGTGGTCGCGGCAATCCGCGAGAGATACTGCCAGCTGGCTGGCGGCCACCTgg                                                                                                                                                  <  1:80187/150‑1 (MQ=255)
 gcgTAGGCGACGTGCAGCTGTTTGGCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACACGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTgg                                                                                                                                                  <  2:2486841/150‑1 (MQ=255)
                   ggTTTGGCATGGGCACATAACCCCTGCGCGTCTGGCTGCACCAGAACAATACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAg                                                                                                                                <  2:807564/149‑1 (MQ=255)
                    gTTTGGCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAgg                                                                                                                               >  1:2477372/1‑150 (MQ=255)
                     tttGGCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGc                                                                                                                              >  1:2581464/1‑150 (MQ=255)
                       tGGCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTc                                                                                                                            >  2:1792530/1‑150 (MQ=255)
                         gCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGa                                                                                                                          <  1:2293520/150‑1 (MQ=255)
                         gCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGa                                                                                                                          <  1:1669516/150‑1 (MQ=255)
                            tGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCa                                                                                                                       >  1:2232509/1‑150 (MQ=255)
                            tGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCa                                                                                                                       >  2:2063692/1‑150 (MQ=255)
                              ggCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGc                                                                                                                     <  2:2273783/150‑1 (MQ=255)
                                cGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGAAGCCTTCCCCAGGCTCGACCAGCtt                                                                                                                   >  2:1478076/1‑150 (MQ=255)
                                 gACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGc                                                                                                                                     >  2:2895434/1‑131 (MQ=255)
                                 gACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGc                                                                                                                                     <  1:2895434/131‑1 (MQ=255)
                                 gACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTc                                                                                                                  >  1:1941070/1‑150 (MQ=255)
                                   ctactCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCa                                                                                                                >  2:725466/1‑150 (MQ=255)
                                    tactCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTggg                                                                                                                                                 <  2:1394063/116‑1 (MQ=255)
                                    tactCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAg                                                                                                               >  1:1060754/1‑150 (MQ=255)
                                    tactCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGCTGGCCGGCGGCCACCTggcggtt                                                                                                                                             >  1:1396713/1‑115 (MQ=255)
                                    tactCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGCTGGCCGGCGGCCACCTggcggtt                                                                                                                                             <  2:1396713/120‑6 (MQ=255)
                                     actCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGc                                                                                                              <  2:1941070/150‑1 (MQ=255)
                                      ctCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCCAGCTGGCCGGCGGCCACCTggcgg                                                                                                                                               >  1:1144977/1‑114 (MQ=38)
                                      ctCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCCAGCTGGCCGGCGGCCACCTggcgg                                                                                                                                               <  2:1144977/117‑4 (MQ=38)
                                      ctCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCt                                                                                                             >  2:2462095/1‑150 (MQ=255)
                                      ctCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGACAGACCAGCTTCCAGCt                                                                                                             >  2:2790644/1‑150 (MQ=255)
                                      ctCGCTGCGCCTGTGTCGTAACCTGATCAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCt                                                                                                             <  2:2477372/150‑1 (MQ=255)
                                       tCGCTGCGCGTCTGGCTGGACGCGAACAACACCGCTTCGCGCAACCTGCCGGCCAGTGATGTGGTCGCGGCAATCCACGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTg                                                                                                            >  2:1298910/1‑150 (MQ=255)
                                       tCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCCGGTGGCCGCCGGCCAGCAGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTg                                                                                                            >  1:32739/1‑150 (MQ=255)
                                       tCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTg                                                                                                            >  1:49327/1‑150 (MQ=255)
                                        cGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCCGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGt                                                                                                           <  1:2790644/150‑1 (MQ=255)
                                         gCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCAGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTc                                                                                                          <  2:32739/150‑1 (MQ=255)
                                           tGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATACGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGa                                                                                                        >  1:2181191/1‑150 (MQ=255)
                                            gcgcGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGAt                                                                                                       <  1:2937261/150‑1 (MQ=255)
                                             cgcgTCTGGCTGGACCCGCACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATc                                                                                                      <  1:2462095/150‑1 (MQ=255)
                                              gcgTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCa                                                                                                     <  2:1019250/150‑1 (MQ=255)
                                               ttccgatctggACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCaa                                                                                                    <  2:2983574/143‑1 (MQ=255)
                                                  ctggctggACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATAc                                                                                                 <  1:222263/150‑1 (MQ=255)
                                                     gctggACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCa                                                                                              <  1:2269305/150‑1 (MQ=255)
                                                     gctggACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCa                                                                                              >  1:1996101/1‑150 (MQ=255)
                                                      ctggACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCaa                                                                                             <  1:2976631/150‑1 (MQ=255)
                                                       tggACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAg                                                                                            <  2:2063408/150‑1 (MQ=255)
                                                       tggACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTAGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAg                                                                                            <  1:2690108/150‑1 (MQ=255)
                                                         gACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATAccc                                                                                               >  2:1030756/1‑145 (MQ=255)
                                                         gACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATAccc                                                                                               <  1:1030756/145‑1 (MQ=255)
                                                          aCCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGgcc                                                                                         >  1:1528384/1‑150 (MQ=255)
                                                            ccGAACAAGACTGCTTCGCGCAACCTGACGGCCGGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGTTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGgccgc                                                                                       >  1:1413197/1‑150 (MQ=255)
                                                            ccGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGCTGGCCGGCGGCCACCt                                                                                                                                                    <  1:1943613/89‑1 (MQ=255)
                                                              gAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCt                                                                                     >  2:619818/1‑150 (MQ=255)
                                                               aaCAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTg                                                                                    <  2:2581464/150‑1 (MQ=255)
                                                                aCAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTgg                                                                                   <  2:2495411/150‑1 (MQ=255)
                                                                    gACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCaa                                                                               <  1:2964493/150‑1 (MQ=255)
                                                                     aCCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAAc                                                                              <  1:1838705/150‑1 (MQ=255)
                                                                     aCCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAAc                                                                              >  2:1744036/1‑150 (MQ=255)
                                                                       cGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGa                                                                            <  1:1136330/150‑1 (MQ=255)
                                                                       cGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGa                                                                            <  1:2063692/150‑1 (MQ=255)
                                                                         cTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCTGGACACCTGGGCGCGCCGCCGGCGCGAGGCTCGACCAGCTTCCAGCTGTCGATCAATGCCCAAGGCCGCCTGGTCAACGAgg                                                                          <  1:1148125/150‑1 (MQ=255)
                                                                          ttCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGa                                                                         >  1:1511629/1‑150 (MQ=255)
                                                                             gcgcAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAaga                                                                      >  1:616519/1‑150 (MQ=255)
                                                                             gcgcAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAaga                                                                      <  1:337394/150‑1 (MQ=255)
                                                                              cgcAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAagag                                                                     >  2:1742159/1‑150 (MQ=255)
                                                                                 aaCCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTc                                                                  <  2:2546919/150‑1 (MQ=255)
                                                                                    cTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTg                                                                                                                                                   <  1:1118393/66‑1 (MQ=255)
                                                                                    cTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTg                                                                                                                                                   >  2:1118393/1‑66 (MQ=255)
                                                                                     tGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGaaa                                                              <  2:49327/150‑1 (MQ=255)
                                                                                       aCGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAAc                                                            <  1:725466/150‑1 (MQ=255)
                                                                                         ggCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAg                                                                                                                                >  1:2595021/1‑80 (MQ=255)
                                                                                         ggCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAg                                                                                                                                <  2:2595021/80‑1 (MQ=255)
                                                                                          gCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAg                                                                                                               >  2:2429911/1‑96 (MQ=255)
                                                                                          gCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAAcatc                                                         >  1:1768650/1‑150 (MQ=255)
                                                                                           ccAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATGAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAGcatca                                                        >  1:2179539/1‑150 (MQ=255)
                                                                                                 aTGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATcc                                                  <  1:889977/150‑1 (MQ=255)
                                                                                                     ggTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCgcgc                                              <  2:2232509/150‑1 (MQ=255)
                                                                                                       tCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCg                                            >  1:1653165/1‑150 (MQ=255)
                                                                                                        cgcgGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCgg                                           <  2:1511629/150‑1 (MQ=255)
                                                                                                          cgGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTg                                         <  1:889542/150‑1 (MQ=255)
                                                                                                              aaTCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGa                                     >  2:1865734/1‑150 (MQ=255)
                                                                                                              aaTCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGa                                     >  2:2207698/1‑150 (MQ=255)
                                                                                                                tCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGgcc                                                                                         >  1:2656879/1‑96 (MQ=255)
                                                                                                                tCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGgcc                                                                                         <  2:2656879/96‑1 (MQ=255)
                                                                                                                tCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATg                                   <  1:1792530/150‑1 (MQ=255)
                                                                                                                tCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATg                                   >  1:1266310/1‑150 (MQ=255)
                                                                                                                 ccGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGTCGATgg                                  <  2:1060754/150‑1 (MQ=255)
                                                                                                                  cgcgAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCTCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAg                                                                           >  1:59933/1‑108 (MQ=255)
                                                                                                                  cgcgAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAgg                                                                          <  2:927281/109‑1 (MQ=255)
                                                                                                                  cgcgAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAgg                                                                          >  1:927281/1‑109 (MQ=255)
                                                                                                                  cgcgAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAg                                                                           <  2:59933/108‑1 (MQ=255)
                                                                                                                   gcgAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAg                                                                           <  2:2791718/107‑1 (MQ=255)
                                                                                                                   gcgAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAg                                                                           >  1:2791718/1‑107 (MQ=255)
                                                                                                                      aGCAGAACCGCCAGGTGGCCGCCGGCCAGCGGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCgaga                             <  2:193063/150‑1 (MQ=255)
                                                                                                                         aGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCa                          <  2:141967/150‑1 (MQ=255)
                                                                                                                              ccGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCg                      >  1:2610634/1‑150 (MQ=255)
                                                                                                                                     ggTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCGCCTGaa                <  2:1768650/150‑1 (MQ=255)
                                                                                                                                      gTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGTTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCGCCTGgaa               <  1:2891892/150‑1 (MQ=255)
                                                                                                                                      gTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCGCCTGaaa               <  2:2048929/150‑1 (MQ=255)
                                                                                                                                      gTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCGCCTGaaa               <  1:619818/150‑1 (MQ=255)
                                                                                                                                       tGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCGCCTGAAAg              <  1:1478076/150‑1 (MQ=255)
                                                                                                                                          ccgccgGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCGCCTGAAAGACa           <  1:1744036/150‑1 (MQ=255)
                                                                                                                                              cgGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCGCCTGAAAGACATCGc       >  2:2663898/1‑150 (MQ=255)
                                                                                                                                                  cAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTACCGATGGCGAGATCACCCGCCTGAAAGACATCGCCCgg   >  2:2913648/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCGTAGGCGACGTGCAGCTGTTTGGCATGGGCGACTACTCGCTGCGCGTCTGGCTGGACCCGAACAAGACCGCTTCGCGCAACCTGACGGCCAGTGATGTGGTCGCGGCAATCCGCGAGCAGAACCGCCAGGTGGCCGCCGGCCAGCTGGGCGCGCCGCCTGCCCCAGGCTCGACCAGCTTCCAGCTGTCGATCAATACCCAAGGCCGCCTGGTCAACGAGGAAGAGTTCGAAAACATCATCATCCGCGCCGGTGCCGATGGCGAGATCACCCGCCTGAAAGACATCGCCCGGG  >  NC_002947/3878726‑3879020

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