Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,821,673 G→A 6.0% Q3173* (CAA→TAA)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,821,6730GA 6.0% 200.9 / 2.3 66Q3173* (CAA→TAA) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base A (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CCAGCGGGAAGTCCGATGGGGTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAG  >  NC_002947/4821526‑4821821
                                                                                                                                                    |                                                                                                                                                    
ccAGCGGGAAGTCCGATGGGGTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGc                                                                                                                                                     <  5:225551/149‑1 (MQ=255)
 cAGCGGGAAGTCCGATGGGGTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTg                                                                                                                                                      >  3:321553/1‑147 (MQ=255)
  aGCGGGAAGTCCGATGGGGTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCag                                                                                                                                                   <  4:17295/149‑1 (MQ=255)
      ggAAGTCCGATGGGGTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCa                                                                                                                                               >  7:60535/1‑149 (MQ=255)
       gAAGTCCGATGGGGTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCAc                                                                                                                                              >  8:92405/1‑149 (MQ=255)
       gAAGTCCGATGGGGTGAAGCTGCCATTGCCATCGGCAAGGCAATGAGCGACCAGCGCCAGCAGCTCCTGACGATAGACCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTACAGAGCAc                                                                                                                                              >  5:240654/1‑149 (MQ=255)
               aTGGGGTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATa                                                                                                                                       >  5:80278/1‑148 (MQ=255)
                  ggggGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGGTCCAGCGCAGTTGCAGAGCACCGCCATacac                                                                                                                                    <  1:233441/149‑1 (MQ=255)
                    gTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTg                                                                                                                                 <  3:364295/149‑1 (MQ=255)
                    gTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTg                                                                                                                                 >  8:16678/1‑149 (MQ=255)
                       aaGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGcc                                                                                                                              <  7:183892/149‑1 (MQ=255)
                        aGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGcc                                                                                                                              <  4:222556/148‑1 (MQ=255)
                                     atCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAAc                                                                                                                 >  6:30133/1‑148 (MQ=255)
                                      tCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCa                                                                                                               <  8:60535/149‑1 (MQ=255)
                                          ccAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGCGCACCGCCATACACCTGGCCATCTACGCTCAACCAGTTg                                                                                                           >  7:130962/1‑149 (MQ=255)
                                          ccAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTg                                                                                                           <  1:124572/149‑1 (MQ=255)
                                          ccAGGCAATCCGCGACCAGCGCCAGCAGCTCCATAACATATCCCTCAGCCAGACGAGCCAAGGTATACTGGTCGTATCCCTCGGAGCTCAACCTCCAGCGCAGTTGCACAGCACATACAAAAACCAGCCCATAGACCCTCAAACAGt                                                                                                             >  4:148788/1‑147 (MQ=255)
                                           cAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCCGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATac                                                                                                                                      >  7:243894/1‑121 (MQ=255)
                                           cAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCCGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATac                                                                                                                                      <  8:243894/121‑1 (MQ=255)
                                            aGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTggg                                                                                                         >  2:29185/1‑149 (MQ=255)
                                                aaTGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAAc                                                                                                     >  1:247400/1‑149 (MQ=255)
                                                        aCCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATc                                                                                             <  7:92405/149‑1 (MQ=255)
                                                        aCCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATc                                                                                             <  5:30133/149‑1 (MQ=255)
                                                          cAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATc                                                                                             >  8:309430/1‑147 (MQ=255)
                                                           aGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAg                                                                                          <  3:148788/149‑1 (MQ=255)
                                                           aGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAg                                                                                          <  4:337864/149‑1 (MQ=255)
                                                           aGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAg                                                                                          <  3:192267/149‑1 (MQ=255)
                                                            gcgcCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAgg                                                                                         >  6:141197/1‑149 (MQ=255)
                                                                cagcagGTCCTGACGATAGCCCTCAGCCAGTCGCGGGATCGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTACAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCTACGGCGCATCGAGGTCAt                                                                                     <  6:240654/149‑1 (MQ=255)
                                                                cagcagCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCAt                                                                                     >  2:301274/1‑149 (MQ=255)
                                                                   cagcTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATgcg                                                                                  >  1:323342/1‑149 (MQ=255)
                                                                          tGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATgcgc                                                                                 <  6:110669/143‑1 (MQ=255)
                                                                          tGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATgcgc                                                                                 >  5:110669/1‑143 (MQ=255)
                                                                             cGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGc                                                                        <  6:80278/149‑1 (MQ=255)
                                                                                tAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGa                                                                     <  4:149805/149‑1 (MQ=255)
                                                                                   cccTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTACAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGCCATGCGCCAACCCGGCCGAgg                                                                   >  7:237662/1‑148 (MQ=255)
                                                                                    ccTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAg                                                                    <  4:218734/146‑1 (MQ=255)
                                                                                    ccTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGc                                                                  >  3:218734/1‑148 (MQ=255)
                                                                                     cTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAAc                                                                                                                 <  8:163579/100‑1 (MQ=255)
                                                                                     cTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAAc                                                                                                                 >  7:163579/1‑100 (MQ=255)
                                                                                        aGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTc                                                               >  6:233799/1‑147 (MQ=255)
                                                                                         gCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAg                                                            <  1:301274/149‑1 (MQ=255)
                                                                                          ccAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGGATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAgg                                                           <  6:219205/149‑1 (MQ=255)
                                                                                          ccAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAgg                                                           <  7:309430/149‑1 (MQ=255)
                                                                                          ccAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAgg                                                           >  3:17430/1‑149 (MQ=255)
                                                                                                gcgcgATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGTGCCAACCCGGCCGAGGCCTCCAGGGGCTg                                                      >  2:272436/1‑148 (MQ=255)
                                                                                                  gcgATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGaa                                                   <  3:297273/149‑1 (MQ=255)
                                                                                                   cgATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAAc                                                  <  2:323342/149‑1 (MQ=255)
                                                                                                     aTGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAAc                                                  >  7:250201/1‑147 (MQ=255)
                                                                                                         tCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGcc                                            <  6:164692/149‑1 (MQ=255)
                                                                                                         tCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGcc                                            >  3:197397/1‑149 (MQ=255)
                                                                                                                gTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTg                                      <  8:250201/148‑1 (MQ=255)
                                                                                                                   gTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCg                                   >  3:89190/1‑148 (MQ=255)
                                                                                                                    tAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGaa                                 >  1:207027/1‑149 (MQ=255)
                                                                                                                    tAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATAGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGaa                                 >  2:313292/1‑149 (MQ=255)
                                                                                                                       cgcTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCa                                         <  3:197597/138‑1 (MQ=255)
                                                                                                                       cgcTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCa                                         >  4:197597/1‑138 (MQ=255)
                                                                                                                       cgcTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCAACGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAActg                              >  8:133869/1‑149 (MQ=255)
                                                                                                                        gcTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAActgc                             >  5:383913/1‑149 (MQ=255)
                                                                                                                         cTCGGCGCTGAAGGTCCAGCGCAGTTGCATAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAActgct                            <  3:285669/149‑1 (MQ=255)
                                                                                                                                    aGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTgg                 <  4:197397/149‑1 (MQ=255)
                                                                                                                                      ggTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAActgc                             >  7:205466/1‑136 (MQ=255)
                                                                                                                                      ggTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAActgc                             <  8:205466/136‑1 (MQ=255)
                                                                                                                                        tCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGccc              <  4:17430/149‑1 (MQ=255)
                                                                                                                                          cAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAg            <  8:303177/149‑1 (MQ=255)
                                                                                                                                             cgcAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTa         <  2:247400/149‑1 (MQ=255)
                                                                                                                                                aGTTACAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGCCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGtt      <  8:237662/149‑1 (MQ=255)
                                                                                                                                                    gCATAGCACCGCCATACTCCTGGCCATCGTCGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCTTGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGaa   <  4:321553/148‑1 (MQ=255)
                                                                                                                                                    gCATAGCACCGCCATACACCTGGCCAGCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAg  <  8:130962/149‑1 (MQ=255)
                                                                                                                                                    gCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAg  >  3:168592/1‑149 (MQ=255)
                                                                                                                                                    gCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAg  >  1:340788/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGCGGGAAGTCCGATGGGGTGAAGCTGCCATTGCCATCGGCCAGGCAATGCGCGACCAGCGCCAGCAGCTCCTGACGATAGCCCTCAGCCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAG  >  NC_002947/4821526‑4821821

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