Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,758,730 T→G 7.3% S352R (AGC→CGC)  pvdP ← pyoverdine biosynthesis‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,758,7300TG 7.3% 195.3 / 3.1 65S352R (AGC→CGC) pvdPpyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base G (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCGCGCCATCGGCACCGGCATGCCGTTGCTCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTAACGCGCGAAA  >  NC_002947/4758587‑4758878
                                                                                                                                               |                                                                                                                                                    
cGTGCGCGCCATCGGCACCGGCATGCCGTTGCTCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAAc                                                                                                                                                 <  3:132896/149‑1 (MQ=255)
  tGCGCGCCATCGGCACCGGCATGCCGTTGCTCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAAcgc                                                                                                                                               <  7:92210/149‑1 (MQ=255)
           cgagCAACGGCATGCCGGTGCTCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCca                                                                                                                                       >  5:22533/4‑148 (MQ=255)
                  cGGCATGCCGTTGCTCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCg                                                                                                                                >  5:96194/1‑148 (MQ=255)
                  cGGCATGCCGTTGCTCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCAGTAACGCGACTCCCACACCTCg                                                                                                                                >  6:340361/1‑148 (MQ=255)
                        gCCGTTGCTCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTcc                                                                                                                                         <  5:1007/133‑1 (MQ=255)
                        gCCGTTGCTCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTcc                                                                                                                                         >  6:1007/1‑133 (MQ=255)
                                tCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCTAACGTCAGCCGGGACAGGTACTGCGGGTCGCTGTAACGCCACTCCCACACCTCGAAATTCCCATGAAAg                                                                                                                 >  8:105622/1‑149 (MQ=255)
                                tCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAg                                                                                                                 >  4:158901/1‑149 (MQ=255)
                                        cGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGgc                                                                                                         >  8:240631/1‑149 (MQ=255)
                                          tGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGgcgc                                                                                                       >  6:97485/1‑149 (MQ=255)
                                           gggACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCt                                                                                                      >  3:66496/1‑149 (MQ=255)
                                              aCTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCtttt                                                                                                   <  6:396996/149‑1 (MQ=255)
                                               cTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTg                                                                                                  <  5:97158/149‑1 (MQ=255)
                                                  aGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGata                                                                                               >  5:408066/1‑149 (MQ=255)
                                                   ggCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGatat                                                                                              >  6:406888/1‑149 (MQ=255)
                                                    gCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATc                                                                                             <  6:22533/149‑1 (MQ=255)
                                                      ccAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATcgc                                                                                           >  8:67796/1‑149 (MQ=255)
                                                       cAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATcgcg                                                                                          <  6:408066/149‑1 (MQ=255)
                                                       cAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGAGTCGCGTTACAGCGACcccc                                                                                                                                        <  6:418385/103‑1 (MQ=255)
                                                       cAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGAGTCGCGTTACAGCGACcccc                                                                                                                                        >  5:418385/1‑103 (MQ=255)
                                                                  cAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTgg                                                                                <  3:158901/148‑1 (MQ=255)
                                                                   agacAGGCGTGCAAGCACCTTTCAAGCTCCGAGCCGAAATGCCCCCACGTCAGCCGGGACAGGTACTGGTGGTCGCGGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGATTTTGATATCGCGCACCCACTGGgt                                                                              <  4:59/146‑1 (MQ=255)
                                                                   aGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGgt                                                                              <  2:236439/149‑1 (MQ=255)
                                                                     ccAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCAAGAAAGGTTTCgg                                                                                                          >  8:202931/1‑119 (MQ=255)
                                                                     ccAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCAAGAAAGGTTTCgg                                                                                                          <  7:202931/119‑1 (MQ=255)
                                                                       aGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTAc                                                                          <  1:326259/149‑1 (MQ=255)
                                                                       aGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTAc                                                                          <  7:67796/149‑1 (MQ=255)
                                                                       aGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTAc                                                                          >  7:365773/1‑149 (MQ=255)
                                                                                      gTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTaa                                                           <  6:96194/149‑1 (MQ=255)
                                                                                       ttCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACCCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAg                                                          >  6:186925/1‑149 (MQ=255)
                                                                                          aaGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGcc                                                        >  4:380937/1‑148 (MQ=255)
                                                                                          aaGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCa                                                       <  2:72632/149‑1 (MQ=255)
                                                                                                     cGAACTTCCCTAAATTCATCCTGGACAGGTACTGGGGTTCGCTGTAACGCGAATCCCACACCTCGATATTGCCATGAAAGGTTTCTTCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTtggt                                            <  7:105622/149‑1 (MQ=255)
                                                                                                            ccccAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAgg                                     >  6:394854/1‑149 (MQ=255)
                                                                                                            ccccAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAg                                      <  5:97485/148‑1 (MQ=255)
                                                                                                            ccccAACGTCAGCCGGGACAGGTACAGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAgg                                     <  2:43343/149‑1 (MQ=255)
                                                                                                             cccAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGtca                                                                    >  1:365722/1‑117 (MQ=255)
                                                                                                             cccAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGtca                                                                    <  2:365722/117‑1 (MQ=255)
                                                                                                             cccAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGCTTGGTGGCAAGGc                                    >  3:211201/1‑149 (MQ=255)
                                                                                                              ccAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCg                                   >  4:139812/1‑149 (MQ=255)
                                                                                                               cAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGt                                  >  6:61061/1‑149 (MQ=255)
                                                                                                                    tCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCg                             >  1:202285/1‑149 (MQ=255)
                                                                                                                      aGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGt                            <  3:380937/148‑1 (MQ=255)
                                                                                                                      aGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAACGCCTTGccac                            <  5:406888/148‑3 (MQ=255)
                                                                                                                      aGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACGCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGt                            <  3:18321/148‑1 (MQ=255)
                                                                                                                       gCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCa                          >  4:44790/1‑149 (MQ=255)
                                                                                                                          gggACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGa                       <  5:394854/149‑1 (MQ=255)
                                                                                                                              cAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTAt                   <  4:211201/149‑1 (MQ=255)
                                                                                                                               aGGTACTGGGGGTCGCTGTAACTCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCTGCGCTTTTGATATCTCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCATGGTTGTGGCAAGGCGTTGCCGTCATAATTATc                  <  6:24283/149‑1 (MQ=255)
                                                                                                                               aGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCAt                         <  7:223701/142‑1 (MQ=255)
                                                                                                                               aGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCAt                         >  8:223701/1‑142 (MQ=255)
                                                                                                                                ggTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCg                 <  7:240631/149‑1 (MQ=255)
                                                                                                                                     tGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCAATGAAAGGTTTCGGCGCTTTTGATATCGCGCACAAACTGGGTGTACTCGTCATCATCGTAAGCCAACCATGTTGGTGGCAAGGAGTTGCCGTCATGATTATCGAAGt             >  5:208117/1‑148 (MQ=255)
                                                                                                                                     tGGGAGTCGCGTTACAGCGACCCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGt             >  3:349882/1‑148 (MQ=255)
                                                                                                                                     tGGGAGTCGCGTTACAGCGAC‑CCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTaa           >  7:249081/1‑149 (MQ=255)
                                                                                                                                      gggggTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTAGCCGTCATGATTATCGAAGTaa           >  6:333173/1‑149 (MQ=255)
                                                                                                                                       ggggTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTa            <  3:44790/147‑1 (MQ=255)
                                                                                                                                        gggTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTAAcg         >  6:187177/1‑149 (MQ=255)
                                                                                                                                           tCGCTGTAGCGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTAAcgcgc      >  5:210363/1‑149 (MQ=255)
                                                                                                                                           tCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTAAcgcgc      <  5:187177/149‑1 (MQ=255)
                                                                                                                                           tCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTAAcgcgc      <  2:202285/149‑1 (MQ=255)
                                                                                                                                             gCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTAACGCGCGa    <  4:66496/149‑1 (MQ=255)
                                                                                                                                              cTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGt                                  <  7:221228/118‑1 (MQ=255)
                                                                                                                                              cTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGt                                  >  8:221228/1‑118 (MQ=255)
                                                                                                                                               tGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTAACGCGCGatc  >  3:373558/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CGTGCGCGCCATCGGCACCGGCATGCCGTTGCTCGGGTCACGTGGGACTGAGGCCCAGCGCATGTGCAGCCAGTCGTGCAAGCCCAGTTCAAGCTCCGAGCCGAACTGCCCCAACGTCAGCCGGGACAGGTACTGGGGGTCGCTGTAACGCGACTCCCACACCTCGAAATTGCCATGAAAGGTTTCGGCGCTTTTGATATCGCGCACCCACTGGGTGTACTCGTCATCATCGTAAGCCAACCAGGTTGGTGGCAAGGCGTTGCCGTCATGATTATCGAAGTAACGCGCGAAA  >  NC_002947/4758587‑4758878

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