Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,758,192 A→G 8.2% intergenic (+7/+14) 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,1920AG 8.2% 317.3 / 12.7 109intergenic (+7/+14)ompQ/pvdPouter membrane pyoverdine efflux protein/pyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base A (50/50);  new base G (4/5);  total (54/55)
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.85e-01

Only 100 of 110 total aligned reads displayed.

CTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGATCCATCTCC  >  NC_002947/4758045‑4758341
                                                                                                                                                   |                                                                                                                                                     
ctaaCCGTGCTCGACACCCCGCGCAGCCCGTACCAGGGCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCTCACGGCAGTCCAGCGTGGCCCTTTACAAACCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAgg                                                                                                                                                     <  1:520199/147‑1 (MQ=255)
 tGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGTTTGCCCTCTTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTATGGCGTCAggg                                                                                                                                                    <  2:1230141/150‑2 (MQ=255)
     cGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATc                                                                                                                                                >  1:1751799/1‑150 (MQ=255)
     cGTGCTCGACACCCAGCGCAGCCTGTACCACGCCCAGGGCCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAACGCCTTGGGCGGTAGCTGGCAAGCAGATAGCCTGACGCCCTAGCGCGTCGGGCTATc                                                                                                                                                <  1:939623/150‑1 (MQ=255)
        gCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTgg                                                                                                                                             >  1:2151902/1‑150 (MQ=255)
         cTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTGTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGgt                                                                                                                                            <  2:788256/150‑1 (MQ=255)
             acacCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTc                                                                                                                                        <  1:1317147/150‑1 (MQ=255)
             acacCCAGCGCAGCCTGTACCAGGCCCAGGACAAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTc                                                                                                                                        >  2:300719/1‑150 (MQ=255)
                cccAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACg                                                                                                                                     >  1:2082637/1‑150 (MQ=255)
                 ccAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGTGCGTCAGGCTATCTGGTGTTCACGa                                                                                                                                    >  1:883206/1‑150 (MQ=255)
                   aGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGAcc                                                                                                                                  >  1:378244/1‑150 (MQ=255)
                        aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                                             <  1:1031084/150‑1 (MQ=255)
                          ccTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCt                                                                                                                           >  2:319224/1‑150 (MQ=255)
                            tGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGa                                                                                                                         <  1:884415/150‑1 (MQ=255)
                            tGTACCAGGCCCAGGACCAGCAGACCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGa                                                                                                                         <  1:893201/150‑1 (MQ=255)
                             gTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGAc                                                                                                                        >  2:869801/1‑150 (MQ=255)
                              tACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGAcc                                                                                                                       <  1:1527676/150‑1 (MQ=255)
                                 cAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTc                                                                                                                    >  1:1505965/1‑150 (MQ=255)
                                 cAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTc                                                                                                                    >  2:759669/1‑150 (MQ=255)
                                  aGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGTCTATCTGGTGTTCACGACCGCTGCCTGACCTTCa                                                                                                                   <  1:2199116/150‑1 (MQ=255)
                                   ggCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAg                                                                                                                  >  2:687064/1‑150 (MQ=255)
                                   ggCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAg                                                                                                                  <  2:137707/150‑1 (MQ=255)
                                      ccAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGTTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTg                                                                                                               >  2:2226143/1‑150 (MQ=255)
                                      ccAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTg                                                                                                               <  1:997996/150‑1 (MQ=255)
                                      ccAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTg                                                                                                               >  2:2185382/1‑150 (MQ=255)
                                          gACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                           <  2:2151902/150‑1 (MQ=255)
                                          gACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                           >  2:1603662/1‑150 (MQ=255)
                                          gACCACCAGGCCCGGCTGTAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                           <  2:954411/150‑1 (MQ=255)
                                             cagcagGCCCGCCTGCAGCTAAGCCCACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCa                                                                                                        >  1:62779/1‑150 (MQ=255)
                                             cagcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCa                                                                                                        >  1:995246/1‑150 (MQ=255)
                                              agcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGcc                                                                                                                            <  2:2122554/129‑1 (MQ=255)
                                              agcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGcc                                                                                                                            >  1:2122554/1‑129 (MQ=255)
                                               gcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATa                                                                                                      >  2:1307650/1‑150 (MQ=255)
                                                  ggCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCa                                                                                                   <  1:1001551/150‑1 (MQ=255)
                                                        ccTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGc                                                                                                                                                    <  2:1258481/95‑1 (MQ=255)
                                                        ccTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGc                                                                                                                                                    >  1:1258481/1‑95 (MQ=255)
                                                         cTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGggcgg                                                                                            >  2:2177992/1‑150 (MQ=255)
                                                            cAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGATAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgc                                                                                         <  1:1332966/150‑1 (MQ=255)
                                                             aGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGTGGGTAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgcg                                                                                        >  2:1957187/1‑150 (MQ=255)
                                                              gCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgcgg                                                                                       >  2:102762/1‑150 (MQ=255)
                                                                tAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCGTTCGGTGGTGCCACGCAAGCCGATAACCCGACTCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCGGGGGCGGCGCGGCa                                                                                     >  1:1213363/1‑150 (MQ=255)
                                                                    ccAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                                 <  1:2185382/150‑1 (MQ=255)
                                                                    ccAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                                 <  1:1603662/150‑1 (MQ=255)
                                                                    ccAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                                 >  2:1711872/1‑150 (MQ=255)
                                                                    ccAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCAAttt                                                                                 <  1:319224/150‑1 (MQ=255)
                                                                                cAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                           <  1:858960/112‑1 (MQ=255)
                                                                                cAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                           >  2:858960/1‑112 (MQ=255)
                                                                                cAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGccc                                                                     <  1:300719/150‑1 (MQ=255)
                                                                                    gTGGCCCTTTACGAAGCCTTGGGCGGGGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCa                                                                 <  2:361941/150‑1 (MQ=255)
                                                                                      ggCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTTTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCCGGGGCGGCGCGCCACTACCCCATGCAGCCCGTCaaa                                                               <  2:1737477/150‑1 (MQ=255)
                                                                                      ggCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCaaa                                                               <  2:2082637/150‑1 (MQ=255)
                                                                                         ccTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAgcg                                                            >  1:2041189/1‑150 (MQ=255)
                                                                                          cTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCAc                                                                                                                                      <  2:873092/75‑1 (MQ=255)
                                                                                          cTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCAc                                                                                                                                      >  1:873092/1‑75 (MQ=255)
                                                                                          cTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAgcgc                                                           >  2:650288/1‑150 (MQ=255)
                                                                                              aCAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGAc                                                                                                                        <  2:1779769/85‑1 (MQ=255)
                                                                                              aCAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGAc                                                                                                                        >  1:1779769/1‑85 (MQ=255)
                                                                                               cAAAGCCTTGGGCGGTGGCTGGGAAGCGCATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGctcctc                                                      <  2:1340260/150‑1 (MQ=255)
                                                                                                aaaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTATGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGTGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcg                                                     >  2:843473/1‑150 (MQ=255)
                                                                                                aaaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGTCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcg                                                     >  1:2290785/1‑150 (MQ=255)
                                                                                                aaaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcg                                                     >  1:1793371/1‑150 (MQ=255)
                                                                                                 aaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgt                                                    >  1:112208/1‑150 (MQ=255)
                                                                                                   gCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgtcg                                                  <  1:1957187/150‑1 (MQ=255)
                                                                                                      ttGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTg                                               <  2:1751799/150‑1 (MQ=255)
                                                                                                             gTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCGGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCa                                        <  2:1213363/150‑1 (MQ=255)
                                                                                                                  tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgcgc                                   >  1:1285612/1‑150 (MQ=255)
                                                                                                                  tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgcgc                                   >  1:401433/1‑150 (MQ=255)
                                                                                                                    gcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgc                                     >  2:1858422/1‑146 (MQ=255)
                                                                                                                    gcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgc                                     <  1:1858422/146‑1 (MQ=255)
                                                                                                                     cAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAg                                <  1:770825/150‑1 (MQ=255)
                                                                                                                      aaGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAgg                               >  1:1467735/1‑150 (MQ=255)
                                                                                                                       aGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGTTCCTCGTCGCTGTAAATCAGGCGCAAGGc                              <  1:2177992/150‑1 (MQ=255)
                                                                                                                          ggATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAg                           >  2:1792642/1‑150 (MQ=255)
                                                                                                                           gATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGc                          <  2:112208/150‑1 (MQ=255)
                                                                                                                           gATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGc                          <  2:378244/150‑1 (MQ=255)
                                                                                                                            aTAGCCTGACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCt                         >  1:1982024/1‑150 (MQ=255)
                                                                                                                            aTAGCCTGACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCt                         >  2:52686/1‑150 (MQ=255)
                                                                                                                            aTAGCCTGACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCt                         >  2:52387/1‑150 (MQ=255)
                                                                                                                               gaccGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCa                      >  2:385051/3‑150 (MQ=255)
                                                                                                                                cctgACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCaaa                                                               >  2:1271977/4‑108 (MQ=39)
                                                                                                                                cctgACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCaaa                                                               <  1:1271977/105‑1 (MQ=39)
                                                                                                                                cctgACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCAt                     <  1:52387/147‑1 (MQ=255)
                                                                                                                                cctgACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCAt                     <  1:52686/147‑1 (MQ=255)
                                                                                                                                cccGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCAt                     >  2:872539/1‑150 (MQ=255)
                                                                                                                                 ctgACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATg                    <  2:1982024/148‑1 (MQ=255)
                                                                                                                                       cgcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTTTAAATCAGGCGCAAGGCCAGCTTCATGGTTTcc              <  1:186742/150‑1 (MQ=255)
                                                                                                                                       cgcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTcc              <  1:650288/150‑1 (MQ=255)
                                                                                                                                        gcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACtt                                                                                  <  1:2299124/81‑1 (MQ=255)
                                                                                                                                        gcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACtt                                                                                  >  2:2299124/1‑81 (MQ=255)
                                                                                                                                        gcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTCCCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCCTGGTTTCCg             >  2:436519/1‑150 (MQ=255)
                                                                                                                                          tatggCGTCATGCGATCTGTTGTTAACGACAGCTGCCTGACCTTATGGCGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGa           <  2:2136141/147‑1 (MQ=255)
                                                                                                                                          tAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgc                                                                                         >  2:1836555/1‑72 (MQ=255)
                                                                                                                                          tAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgc                                                                                         <  1:1836555/72‑1 (MQ=255)
                                                                                                                                          tAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGa           <  1:872539/150‑1 (MQ=255)
                                                                                                                                          tAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGa           <  1:1276488/150‑1 (MQ=255)
                                                                                                                                          tAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGAGGCACTTTGCCTACCAGCCCGTCAAATCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGa           <  2:572302/150‑1 (MQ=255)
                                                                                                                                           atggCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGTGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGAt          <  1:843473/148‑1 (MQ=255)
                                                                                                                                              gCGTCAGGCTATCTGGTGTTCACGACAGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGATCCa       >  2:396221/1‑150 (MQ=255)
                                                                                                                                               cGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGATCCAt      <  1:2203667/150‑1 (MQ=255)
                                                                                                                                                   aGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCTCGTCAAAGAGCTCTGCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGATCCATCTcc  <  1:286298/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGATCCATCTCC  >  NC_002947/4758045‑4758341

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