Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,880,378 C→G 6.7% A738G (GCC→GGC)  mexF → multidrug RND transporter MexF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,880,3780CG 6.7% 228.2 / 6.0 75A738G (GCC→GGC) mexFmultidrug RND transporter MexF
Reads supporting (aligned to +/- strand):  ref base C (37/33);  new base G (2/3);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

CAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAGGCCGAGCAGCAGTTC  >  NC_002947/3880239‑3880526
                                                                                                                                           |                                                                                                                                                    
cAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCg                                                                                                                                             >  1:44747/1‑149 (MQ=255)
 aaaTCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActtg                                                                                                                                            <  2:337104/149‑4 (MQ=255)
 aaaTCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCgg                                                                                                                                            <  5:130689/149‑1 (MQ=255)
   aTCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActtgtg                                                                                                                                          >  1:173832/1‑144 (MQ=255)
   aTCGAAGACCGGGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActtgtg                                                                                                                                          <  3:304256/149‑6 (MQ=255)
      gAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGaa                                                                                                                                         >  8:357030/1‑147 (MQ=255)
       aaGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActtgt                                                                                                                                           >  6:354222/1‑140 (MQ=255)
       aaGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActtgt                                                                                                                                           <  5:354222/144‑5 (MQ=255)
              tGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCgg                                                                                                                                            <  8:7867/136‑1 (MQ=255)
              tGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCgg                                                                                                                                            >  7:7867/1‑136 (MQ=255)
               ggCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGAc                                                                                                                              >  1:78169/1‑149 (MQ=255)
                ggaaCCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACAAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGAcc                                                                                                                             >  5:223109/3‑149 (MQ=255)
                  aaCCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATc                                                                                                                                                  <  6:173501/126‑1 (MQ=255)
                  aaCCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATc                                                                                                                                                  >  5:173501/1‑126 (MQ=255)
                  aaCCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCa                                                                                                                           >  5:67646/1‑149 (MQ=255)
                   aCCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCAc                                                                                                                          >  5:306590/1‑149 (MQ=255)
                       gggTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGgtg                                                                                                                      >  2:5671/1‑149 (MQ=255)
                         gtgaCGAAGCGCTGGACAAGGAAACCCAGAACATCATCGCCAAGATCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAACTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCGACGGtgtg                                                                                                                    <  4:336343/146‑1 (MQ=255)
                         gTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGtgtg                                                                                                                    <  1:5671/149‑1 (MQ=255)
                            aCGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGc                                                                                                                  >  1:261853/1‑148 (MQ=255)
                             cGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATTCCGCCATCGACCGGGTAAAGGCCAAGACCCACGGtgtg                                                                                                                    >  8:61180/1‑145 (MQ=255)
                             cGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGTAAAGGCCAAGACCCACGGtgtg                                                                                                                    <  7:61180/145‑1 (MQ=255)
                              gAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGAt                                                                                                               >  6:341776/1‑149 (MQ=255)
                                       tACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACa                                                                                                       >  4:311679/1‑148 (MQ=255)
                                           aGGAAACCCAGAACATCATCGCCAAGAGCCTCAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTc                                                                                                  >  1:33084/1‑149 (MQ=255)
                                           aGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTc                                                                                                  <  2:44747/149‑1 (MQ=255)
                                            ggAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTc                                                                                                  >  3:206434/1‑148 (MQ=255)
                                            ggAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCg                                                                                                 <  6:30631/149‑1 (MQ=255)
                                                   cggAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCt                                                                                          >  1:122760/3‑149 (MQ=255)
                                                    aGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGAccc                                                                                                                            <  7:116936/114‑1 (MQ=255)
                                                    aGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGAccc                                                                                                                            >  8:116936/1‑114 (MQ=255)
                                                       aCATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAg                                                                                      <  5:341776/149‑1 (MQ=255)
                                                        catcatCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAgg                                                                                     >  5:15034/1‑149 (MQ=255)
                                                               gCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCt                                                                              <  2:122760/149‑1 (MQ=255)
                                                                ccAAGGGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTg                                                                             <  2:38390/149‑1 (MQ=255)
                                                                ccAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTg                                                                             >  7:139959/1‑149 (MQ=255)
                                                                 cAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTgg                                                                            >  2:186402/1‑149 (MQ=255)
                                                                  aaGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTggg                                                                           >  1:185555/1‑149 (MQ=255)
                                                                    gagCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCt                                                                         <  2:78169/149‑1 (MQ=255)
                                                                        cacaACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTgg                                                                                                                   <  4:396046/103‑1 (MQ=255)
                                                                        cacaACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTgg                                                                                                                   >  3:396046/1‑103 (MQ=255)
                                                                         acaACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTg                                                                    >  2:351491/1‑149 (MQ=255)
                                                                         acaACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTg                                                                    >  2:368367/1‑149 (MQ=255)
                                                                           aaCGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTa                                                                  >  6:88925/1‑149 (MQ=255)
                                                                                 ccAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCa                                                                                                                           <  5:179390/86‑1 (MQ=255)
                                                                                 ccAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCa                                                                                                                           >  6:179390/1‑86 (MQ=255)
                                                                                  cAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAAc                                                           <  4:206434/149‑1 (MQ=255)
                                                                                  cAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAAc                                                           <  5:88925/149‑1 (MQ=255)
                                                                                  cAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAAc                                                           <  1:368367/149‑1 (MQ=255)
                                                                                              gCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCCGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGc                                               >  4:256263/1‑149 (MQ=255)
                                                                                                cTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCtt                                             <  6:306590/149‑1 (MQ=255)
                                                                                                     cACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCaa                                                            >  4:400796/1‑129 (MQ=255)
                                                                                                     cACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCaa                                                            <  3:400796/129‑1 (MQ=255)
                                                                                                        cAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGCTCCCCGACATCTTCGACACCCTGCAGGTTTACATGGGCTCGCTGTACACCAACGACTTCAACCGCTATGGCCGCa                                     >  7:379788/1‑149 (MQ=255)
                                                                                                        cAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCa                                     >  8:232661/1‑149 (MQ=255)
                                                                                                         aGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGc                                      <  4:20357/147‑1 (MQ=255)
                                                                                                         aGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCATTGGCCGc                                      >  3:20357/1‑147 (MQ=255)
                                                                                                              ccAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTAc                                <  1:351491/148‑1 (MQ=255)
                                                                                                               cAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCa                              <  2:33084/149‑1 (MQ=255)
                                                                                                                aGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAg                             <  7:357030/149‑1 (MQ=255)
                                                                                                                aGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCCGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCACACTTACCAg                             <  3:256263/149‑1 (MQ=255)
                                                                                                                 ggTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAgg                            <  6:67646/149‑1 (MQ=255)
                                                                                                                  gTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCg                                                                       >  6:237679/1‑105 (MQ=255)
                                                                                                                  gTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCg                                                                       <  5:237679/105‑1 (MQ=255)
                                                                                                                     aaCGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCaa                        >  7:278928/1‑149 (MQ=255)
                                                                                                                     aaCGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCaa                        >  7:330229/1‑149 (MQ=255)
                                                                                                                     aaCGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCa                         <  8:139959/148‑1 (MQ=255)
                                                                                                                      aCGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCaa                                                   >  5:409748/1‑121 (MQ=255)
                                                                                                                      aCGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCaa                                                   <  6:409748/121‑1 (MQ=255)
                                                                                                                      aCGTGCCACAAGTCGATGCCACCCTCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAAc                       <  8:128149/149‑1 (MQ=255)
                                                                                                                               aaGTCTATGCCGCCATCGACCGGTAAAAGGACAAGACCCACGGTGTGGCTATCTCGGACATCTTCGTCAGCCTGCAGGTTTACCTGGGCTCGCTGTACTCCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAAc                       <  6:383946/140‑1 (MQ=255)
                                                                                                                               aaGTCGATGCCGCCATCGACCGGGAAAAGTCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTATACCTGTGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAAc                       >  5:383946/1‑140 (MQ=255)
                                                                                                                                 gTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAGGcc              >  7:367258/1‑147 (MQ=255)
                                                                                                                                 gTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAGGCCGa            <  7:232661/149‑1 (MQ=255)
                                                                                                                                    gATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAGGCCGagc          <  2:185555/148‑1 (MQ=255)
                                                                                                                                           ccATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAGGCCGAGCAGCAGTTc  >  6:167756/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAGGCCGAGCAGCAGTTC  >  NC_002947/3880239‑3880526

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