Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,758,141 G→C 5.5% intergenic (+1/+20) ompQ → / ← pvdP outer membrane pyoverdine efflux protein/pyoverdine biosynthesis‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,758,1410GC 5.5% 237.0 / 3.2 72intergenic (+1/+20)ompQ/pvdPouter membrane pyoverdine efflux protein/pyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base C (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

GACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGAT  >  NC_002947/4757993‑4758288
                                                                                                                                                    |                                                                                                                                                   
gaCCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAgg                                                                                                                                                     <  3:125886/149‑1 (MQ=255)
gaCCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAgg                                                                                                                                                     <  1:334747/149‑1 (MQ=255)
          aCCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGc                                                                                                                                            >  6:100543/1‑148 (MQ=255)
            cGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTa                                                                                                                                          >  7:150606/1‑148 (MQ=255)
            cGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTAt                                                                                                                                         >  8:254843/1‑149 (MQ=255)
            cGTGCTCGACACCCAGCGCACCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTa                                                                                                                                          >  4:333030/1‑148 (MQ=255)
              tGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCt                                                                                                                                       <  1:211278/149‑1 (MQ=255)
               gCTCGACACCCAGCGCAGCCTGTATCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTg                                                                                                                                      <  7:199176/149‑1 (MQ=255)
               gCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTg                                                                                                                                      <  3:84772/149‑1 (MQ=255)
                cTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACTCGCTAGGGCGTCAGGCTATCTg                                                                                                                                      <  8:310544/148‑1 (MQ=255)
                   gACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGgtg                                                                                                                                   <  2:130218/148‑1 (MQ=255)
                   gACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGgtg                                                                                                                                   >  1:130218/1‑148 (MQ=255)
                     cacCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTc                                                                                                                                <  3:330625/149‑1 (MQ=255)
                        ccAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCAc                                                                                                                              >  4:308942/1‑148 (MQ=255)
                        ccAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACg                                                                                                                             >  6:194153/1‑149 (MQ=255)
                          aGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCa                                                                                                                                               >  4:23119/1‑129 (MQ=255)
                          aGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCa                                                                                                                                               <  3:23119/129‑1 (MQ=255)
                           gcgcAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGc                                                                                                                                                   <  3:253385/124‑1 (MQ=255)
                           gcgcAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGc                                                                                                                                                   >  4:253385/1‑124 (MQ=255)
                             gcAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGGCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGc                                                                                                                        >  7:133413/1‑149 (MQ=255)
                              cAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGCCGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGCCCGCt                                                                                                                       >  8:139677/1‑149 (MQ=255)
                                gCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                                     <  4:195210/149‑1 (MQ=255)
                                     tACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGa                                                                                                                 <  3:308942/148‑1 (MQ=255)
                                               aGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTgg                                                                                                      >  2:107936/1‑149 (MQ=255)
                                                    cagcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCAGATAGCCTGACGCCCTAGCGCGTCGGGCTATcc                                                                                                                                       <  1:19527/111‑2 (MQ=38)
                                                    cagcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCAGATAGCCTGACGCCCTAGCGCGTCGGGCTATcc                                                                                                                                       >  2:19527/1‑110 (MQ=38)
                                                       cagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATa                                                                                              <  3:78742/149‑1 (MQ=255)
                                                        agGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATAc                                                                                             <  5:100543/149‑1 (MQ=255)
                                                          gCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATAc                                                                                             >  8:307514/1‑147 (MQ=255)
                                                           cccGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAg                                                                                          <  7:254843/149‑1 (MQ=255)
                                                            ccGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAgg                                                                                         >  5:70430/1‑149 (MQ=255)
                                                              gCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAgggg                                                                                       >  2:291000/1‑149 (MQ=255)
                                                                    aGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAg                                                                                                                                              >  8:98552/1‑88 (MQ=255)
                                                                    aGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAg                                                                                                                                              <  7:98552/88‑1 (MQ=255)
                                                                       tAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCGGGCTATc                                                                                                                                        <  3:47251/91‑1 (MQ=38)
                                                                       tAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCGGGCTATc                                                                                                                                        >  4:47251/1‑91 (MQ=38)
                                                                         aGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTgg                                                                                                                                     <  6:208595/92‑1 (MQ=255)
                                                                         aGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgcggc                                                                              <  5:194153/147‑1 (MQ=255)
                                                                         aGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgcggc                                                                              <  1:107936/147‑1 (MQ=255)
                                                                           ccAACTGCAGGCCAGCGTGGCCCTTTACAAAGCGTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACtt                                                                          >  6:305128/1‑149 (MQ=255)
                                                                            cAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgc                                                                                 <  4:292733/141‑1 (MQ=255)
                                                                            cAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgc                                                                                 >  3:292733/1‑141 (MQ=255)
                                                                            cAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGGGCGGCACttt                                                                         <  3:333030/149‑1 (MQ=255)
                                                                                  cAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTAc                                                                   >  5:33830/1‑149 (MQ=255)
                                                                                         gCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCg                                                            <  7:307514/149‑1 (MQ=255)
                                                                                               cccTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAg                                                      >  8:310740/1‑149 (MQ=255)
                                                                                               cccTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAg                                                      <  8:150606/149‑1 (MQ=255)
                                                                                                 cTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAgcg                                                    >  6:14173/1‑149 (MQ=255)
                                                                                                         aGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                   >  4:304102/1‑94 (MQ=255)
                                                                                                         aGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                   <  3:304102/94‑1 (MQ=255)
                                                                                                         aGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgt                                            <  6:70430/149‑1 (MQ=255)
                                                                                                         aGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgt                                            <  3:147699/149‑1 (MQ=255)
                                                                                                           ccTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgtcg                                          >  6:30036/1‑149 (MQ=255)
                                                                                                                      ggctggcAAGCGGATAGCCAGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAg                               >  5:308031/1‑149 (MQ=255)
                                                                                                                       gctggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAg                               <  1:291000/148‑1 (MQ=255)
                                                                                                                        ctggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgc                             >  8:224680/1‑149 (MQ=255)
                                                                                                                         tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTaa                                    <  5:145462/141‑1 (MQ=255)
                                                                                                                         tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTaa                                    >  6:145462/1‑141 (MQ=255)
                                                                                                                         tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgcg                            >  3:162029/1‑149 (MQ=255)
                                                                                                                         tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgcg                            >  7:104645/1‑149 (MQ=255)
                                                                                                                           gcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGctc                                                 >  2:92246/1‑126 (MQ=255)
                                                                                                                           gcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGctc                                                 <  1:92246/126‑1 (MQ=255)
                                                                                                                             aaGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAg                        <  6:33830/149‑1 (MQ=255)
                                                                                                                               gCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGc                      <  8:104645/149‑1 (MQ=255)
                                                                                                                                cGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCa                                                                 <  4:322978/105‑1 (MQ=255)
                                                                                                                                cGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCa                                                                 >  3:322978/1‑105 (MQ=255)
                                                                                                                                  gATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGctc                                                 <  7:271171/119‑1 (MQ=255)
                                                                                                                                  gATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGctc                                                 >  8:271171/1‑119 (MQ=255)
                                                                                                                                     aGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCtt                <  5:305128/149‑1 (MQ=255)
                                                                                                                                       cccGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgcggc                                                                              <  2:133954/85‑1 (MQ=255)
                                                                                                                                       cccGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgcggc                                                                              >  1:133954/1‑85 (MQ=255)
                                                                                                                                             gcgcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgtcg                                          <  5:369941/115‑1 (MQ=255)
                                                                                                                                             gcgcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgtcg                                          >  6:369941/1‑115 (MQ=255)
                                                                                                                                               gcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTcc      >  6:174080/1‑149 (MQ=255)
                                                                                                                                                   gggCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGAt  <  5:14173/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGAT  >  NC_002947/4757993‑4758288

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: