Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,050,684 A→G 9.3% V239V (GTA→GTG)  oprB‑III → carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,050,6840AG 9.3% 173.1 / 7.0 64V239V (GTA→GTG) oprB‑IIIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base A (32/26);  new base G (3/3);  total (35/29)
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.71e-01

CCAACGCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTGGTCGACCAGGCG  >  NC_002947/4050536‑4050831
                                                                                                                                                    |                                                                                                                                                   
caaaCGCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGAGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCAGGGGCCCAAGCACACCACCGGCACCGTGGTa                                                                                                                                                     >  6:249668/3‑149 (MQ=255)
       ccACCTCGGCGTGCGCGTGAAGTACGACCTGCCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCgg                                                                                                                                                 >  6:27337/1‑146 (MQ=255)
         aCCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGaa                                                                                                                                            >  3:187843/1‑149 (MQ=255)
          ccTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGTACCACGGTGCCGGtggtgt                                                                                                                                           <  7:105618/149‑6 (MQ=255)
              ggCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTgg                                                                                                                                        <  5:315903/148‑1 (MQ=255)
               gCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGta                                                                                                                                      >  4:141048/1‑149 (MQ=255)
                           aGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAg                                                                                                                          <  3:188443/149‑1 (MQ=255)
                           aGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAg                                                                                                                          <  2:66836/149‑1 (MQ=255)
                            gTACGACCTGCCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAgg                                                                                                                         <  5:27337/149‑1 (MQ=255)
                              aCGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAggg                                                                                                                        <  3:245641/148‑1 (MQ=255)
                                     gTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAg                                                                                                                                <  8:81664/133‑1 (MQ=255)
                                     gTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAg                                                                                                                                >  7:81664/1‑133 (MQ=255)
                                            gcctgGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGc                                                                                                           >  1:197271/1‑147 (MQ=255)
                                             cctgGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGa                                                                                                         <  3:337095/148‑1 (MQ=255)
                                                         gcgTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGt                                                                                            >  2:171773/1‑149 (MQ=255)
                                                           gTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTg                                                                                              <  2:186698/145‑1 (MQ=255)
                                                           gTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTg                                                                                              >  1:186698/1‑145 (MQ=255)
                                                               cggcgTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtactac                                                                                      >  2:188795/1‑149 (MQ=255)
                                                                ggcgTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtactac                                                                                      >  2:18242/1‑148 (MQ=255)
                                                                    ttCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGa                                                                                 >  8:104263/1‑149 (MQ=255)
                                                                      cAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGa                                                                                 >  4:286920/1‑147 (MQ=255)
                                                                      cAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGa                                                                                 >  8:9514/1‑147 (MQ=255)
                                                                      cAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCACGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGa                                                                                 >  1:225390/1‑147 (MQ=255)
                                                                       aaCGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtactactac                                                                                   >  4:131526/1‑144 (MQ=255)
                                                                       aaCGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtactactac                                                                                   <  3:131526/144‑1 (MQ=255)
                                                                         cGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGct                                                                            <  3:141048/149‑1 (MQ=255)
                                                                         cGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGct                                                                            <  2:197271/149‑1 (MQ=255)
                                                                         cGTCCACCCCGAAACCAACGGCAACTCCACCCGCCCCTGCCACCTGCGGCCCAAGCACACCACCGGCACCGAGCTACCCCTAAAACACGTAAACAAGCAGCAGGGTGAACTGCCAGGCCAGTACAAGCATGCTTACAACAACCACAGct                                                                            >  8:244229/1‑149 (MQ=255)
                                                                          gTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGctc                                                                           >  7:223836/1‑149 (MQ=255)
                                                                          gTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGctc                                                                           >  4:313892/1‑149 (MQ=255)
                                                                              aCCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGAt                                                                       <  8:223836/149‑1 (MQ=255)
                                                                                 ccGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGAtgt                                                                     >  3:54622/1‑148 (MQ=255)
                                                                                  cGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGAtgtg                                                                    >  3:56450/1‑148 (MQ=255)
                                                                                    aaaGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGaaa                                                                 <  1:18242/149‑1 (MQ=255)
                                                                                         aaCGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCAt                                                            >  4:311228/1‑149 (MQ=255)
                                                                                         aaCGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCAt                                                            >  7:241821/1‑149 (MQ=255)
                                                                                             gCAACTCCAGCCGCGCCTGGCACATGGGGCCCAAGCACACCACCGGCACCCTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACAACTACGACAGCTCTGATGTGAAACGCATCgg                                                         >  2:110929/1‑148 (MQ=255)
                                                                                                aCTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGc                                                     <  4:54622/149‑1 (MQ=255)
                                                                                                    cAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCgatg                                                 <  1:188795/149‑1 (MQ=255)
                                                                                                       ccGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCgatgat                                               >  2:8330/1‑148 (MQ=255)
                                                                                                       ccGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGatgatg                                              <  3:286920/149‑1 (MQ=255)
                                                                                                         gcgcCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGaa                                            <  1:8330/149‑1 (MQ=255)
                                                                                                                     tGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGcc                                  >  7:4332/1‑147 (MQ=255)
                                                                                                                       gggACCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGa                                <  7:244229/147‑1 (MQ=255)
                                                                                                                        ggCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGt                             <  7:9514/149‑1 (MQ=255)
                                                                                                                         gCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGc                                                        <  5:126335/121‑1 (MQ=255)
                                                                                                                         gCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGc                                                        >  6:126335/1‑121 (MQ=255)
                                                                                                                          cccAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCTGGGTGAACTGCCAGGCGTCTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTgg                           <  1:110929/149‑1 (MQ=255)
                                                                                                                             aaGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtacta                                                                                       >  7:181406/1‑86 (MQ=255)
                                                                                                                             aaGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtacta                                                                                       <  8:181406/86‑1 (MQ=255)
                                                                                                                               gCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCga                                                   >  6:346800/1‑120 (MQ=255)
                                                                                                                               gCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCga                                                   <  5:346800/120‑1 (MQ=255)
                                                                                                                                   tctaccGGTACCACGGTGCCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGcc                         >  3:82928/4‑142 (MQ=255)
                                                                                                                                   tctaccGGTACCACGGTGCCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGcc                         <  4:82928/139‑1 (MQ=255)
                                                                                                                                    ccaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACctg                 <  8:4332/149‑1 (MQ=255)
                                                                                                                                     taccGGTACCACGGTGCCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACctga                >  7:266832/2‑148 (MQ=255)
                                                                                                                                     taccGGTACCACGGTGCAGTTAGAACTGGTATACATGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGTCCGAGGTGGCCATTAcat                  <  8:266832/146‑3 (MQ=255)
                                                                                                                                     caccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACctgc                >  2:26288/1‑149 (MQ=255)
                                                                                                                                     caccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGATGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACctgc                <  4:187843/149‑1 (MQ=255)
                                                                                                                                     caccGGCACCGTAGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACctgc                >  5:230895/1‑149 (MQ=255)
                                                                                                                                      accGGTACCACGGTGCCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACctgct               >  4:65769/1‑149 (MQ=255)
                                                                                                                                              cGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTgg             <  6:342912/143‑1 (MQ=255)
                                                                                                                                              cGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTgg             >  5:342912/1‑143 (MQ=255)
                                                                                                                                              cGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTGGTCGAc        >  2:341533/1‑148 (MQ=255)
                                                                                                                                                   tACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACGAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTGGTCGACCAGGCg  >  3:254902/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAACGCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTGGTCGACCAGGCG  >  NC_002947/4050536‑4050831

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: ‑