Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,758,424 C→G 9.1% V454L (GTG→CTG)  pvdP ← pyoverdine biosynthesis‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,758,4240CG 9.1% 193.7 / 9.8 66V454L (GTG→CTG) pvdPpyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base G (3/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

ATGGTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGC  >  NC_002947/4758276‑4758569
                                                                                                                                                     |                                                                                                                                                 
aTGGTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAc                                                                                                                                                   <  4:162029/149‑1 (MQ=255)
aTGGTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAc                                                                                                                                                   <  1:96688/149‑1 (MQ=255)
   gTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGTGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGAc                                                                                                                                                >  2:280796/1‑149 (MQ=255)
    tttCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACt                                                                                                                                               >  3:318462/1‑149 (MQ=255)
            tCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAgg                                                                                                                                       <  5:20984/149‑1 (MQ=255)
               aTCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                    >  4:306722/1‑149 (MQ=255)
                tctcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                    >  5:98531/1‑148 (MQ=255)
                  tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGaa                                                                                                                                 >  1:61858/1‑149 (MQ=255)
                   ccACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAAt                                                                                                                                >  5:211800/1‑149 (MQ=255)
                    cACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGaa                                                                                                                                 <  5:197141/147‑1 (MQ=255)
                    cACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                               <  6:211800/149‑1 (MQ=255)
                     aCGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGTGCCCCGGCAGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                               <  1:280796/149‑1 (MQ=255)
                       gCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGtt                                                                                                                            <  7:267377/149‑1 (MQ=255)
                         tACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAGTCGTGCGCCAgga                                                                                                                                      <  5:299871/137‑2 (MQ=255)
                         tACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAGTCGTGCGCCAgga                                                                                                                                      >  6:299871/1‑136 (MQ=255)
                           cGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCa                                                                                                                          <  4:318462/147‑1 (MQ=255)
                           cGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCa                                                                                                                          <  1:302618/147‑1 (MQ=255)
                           cGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATg                                                                                                                        >  7:106751/1‑149 (MQ=255)
                           cGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATg                                                                                                                        >  7:101545/1‑149 (MQ=255)
                            gCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAc                                                                                                                                                   >  5:132619/1‑121 (MQ=255)
                            gCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAc                                                                                                                                                   <  6:132619/121‑1 (MQ=255)
                               cTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgcg                                                                                                                    >  1:287717/1‑149 (MQ=255)
                                  ggCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCgg                                                                                                                 <  7:27571/149‑1 (MQ=255)
                                     gccgCCCCGCGACGGTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAAc                                                                                                              <  8:101545/146‑1 (MQ=255)
                                     gcAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAAc                                                                                                              <  8:106751/149‑1 (MQ=255)
                                     gcAGCCCCGCCAAGTTACTTCAGCCGTGCGTGTCCGGCCCCAGCCATGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAAc                                                                                                              <  4:244139/149‑1 (MQ=255)
                                          cccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGGGCCAGGGTGGAATCGTTCATGGGCGCGGaa                                                                                                               <  1:30456/143‑1 (MQ=255)
                                          cccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGaa                                                                                                               >  2:30456/1‑143 (MQ=255)
                                            cGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCt                                                                                                       <  8:9320/149‑1 (MQ=255)
                                                       tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcg                                                                                            >  5:328145/1‑149 (MQ=255)
                                                         cAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTc                                                                                                         <  4:41745/134‑1 (MQ=255)
                                                         cAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTc                                                                                                         >  3:41745/1‑134 (MQ=255)
                                                                    gtCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCa                                                                                                                          <  2:342905/106‑1 (MQ=255)
                                                                    gtCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCa                                                                                                                          >  1:342905/1‑106 (MQ=255)
                                                                      ccGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGcc                                                                                                                                          >  4:2068/1‑88 (MQ=255)
                                                                      ccGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGcc                                                                                                                                          <  3:2068/88‑1 (MQ=255)
                                                                         agcccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAg                                                                          <  7:185411/147‑1 (MQ=255)
                                                                              aGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAgcgc                                                                       >  3:160433/1‑147 (MQ=255)
                                                                               gCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAg                                                                    <  1:129274/149‑1 (MQ=255)
                                                                                  aaGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGaa                                                                  >  1:357486/1‑148 (MQ=255)
                                                                                     gATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAaca                                                                >  2:44125/1‑147 (MQ=255)
                                                                                     gATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAaca                                                                >  2:244104/1‑147 (MQ=255)
                                                                                          gCGACCTCCACCCAGCGCCCCGGCGCGATCCATGGCACGCCATTGACCTCCAGGCGCACGACTGCGCCAGGGTGGAATCGTGCATGGGCGCGGAACCAGTCCTCGACACCGTCGTCGATACAGCCATGGAAGCGGCATAACACAggggt                                                         >  8:122436/1‑147 (MQ=255)
                                                                                              ccTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                           >  1:190892/1‑111 (MQ=255)
                                                                                              ccTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                           <  2:190892/111‑1 (MQ=255)
                                                                                               cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGt                                                    >  8:222834/1‑149 (MQ=255)
                                                                                               cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGt                                                    <  1:44125/149‑1 (MQ=255)
                                                                                                  cACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTgg                                                  <  2:61858/148‑1 (MQ=255)
                                                                                                   aCCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGaa                                               <  6:328145/149‑1 (MQ=255)
                                                                                                          gcCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAAc                                         <  3:225985/149‑1 (MQ=255)
                                                                                                            cccGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACgg                                       >  6:82892/1‑149 (MQ=255)
                                                                                                             ccGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACggg                                      >  2:72950/1‑149 (MQ=255)
                                                                                                              cGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAgg                                                            <  6:107189/126‑1 (MQ=255)
                                                                                                              cGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAgg                                                            >  5:107189/1‑126 (MQ=255)
                                                                                                                     gAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTcccc                                 >  1:282080/1‑147 (MQ=255)
                                                                                                                     gAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCga                               <  1:244104/149‑1 (MQ=255)
                                                                                                                             gCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCg                       <  7:222834/149‑1 (MQ=255)
                                                                                                                                        aCCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGATGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgcgc             >  7:114329/1‑148 (MQ=255)
                                                                                                                                           tcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAgg                                                            <  2:298756/98‑1 (MQ=255)
                                                                                                                                           tcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAgg                                                            >  1:298756/1‑98 (MQ=255)
                                                                                                                                           tcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgc                 <  4:270220/141‑1 (MQ=255)
                                                                                                                                           tcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgc                 >  3:270220/1‑141 (MQ=255)
                                                                                                                                            cctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACgg                                       <  6:70930/118‑1 (MQ=255)
                                                                                                                                            cctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACgg                                       <  6:169664/118‑1 (MQ=255)
                                                                                                                                            cctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACgg                                       >  5:70930/1‑118 (MQ=255)
                                                                                                                                            cctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACgg                                       >  5:169664/1‑118 (MQ=255)
                                                                                                                                                 cgcACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAgcg   <  6:10224/149‑1 (MQ=255)
                                                                                                                                                  gcACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAgcgc  >  1:15272/1‑149 (MQ=255)
                                                                                                                                                    aCGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgc               <  2:172672/134‑1 (MQ=255)
                                                                                                                                                    aCGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgc               >  1:172672/1‑134 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ATGGTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGC  >  NC_002947/4758276‑4758569

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: