Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,772,813 A→G 5.5% intergenic (+33/+36) fpvA → / ← PP_4218 TonB‑dependent outer membrane ferripyoverdine receptor FpvA/lipase/esterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,772,8130AG 5.5% 867.6 / 24.1 272intergenic (+33/+36)fpvA/PP_4218TonB‑dependent outer membrane ferripyoverdine receptor FpvA/lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base A (136/120);  new base G (11/4);  total (147/125)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

Only 100 of 276 total aligned reads displayed.

AACAACGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCT  >  NC_002947/4772664‑4772961
                                                                                                                                                     |                                                                                                                                                    
aacaacGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCa                                                                                                                                                      >  2:1863972/1‑150 (MQ=255)
    acGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTg                                                                                                                                                  >  1:4981960/1‑150 (MQ=255)
         ttCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtt                                                                                                                                             <  1:5371317/150‑1 (MQ=255)
            gACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGttttt                                                                                                                                          >  1:4888164/1‑150 (MQ=255)
             aCAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtttttt                                                                                                                                         >  1:4473171/1‑150 (MQ=255)
               aagaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                       <  1:717187/150‑1 (MQ=255)
               aagaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                       >  1:1590867/1‑150 (MQ=255)
                  aagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATgg                                                                                                                                    <  1:2709140/150‑1 (MQ=255)
                    gTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGAt                                                                                                                                  >  2:5435587/1‑150 (MQ=255)
                       ctacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGcc                                                                                                                               >  2:2200483/1‑150 (MQ=255)
                       ctacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGcc                                                                                                                               >  2:3082355/1‑150 (MQ=255)
                       ctacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGcc                                                                                                                               >  2:1298353/1‑150 (MQ=255)
                           cAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGaa                                                                                                                           <  2:4909707/150‑1 (MQ=255)
                            aGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGaaa                                                                                                                          <  1:575372/150‑1 (MQ=255)
                             gggCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAg                                                                                                                         <  2:3409437/150‑1 (MQ=255)
                                 aTTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTc                                                                                                                     <  2:5140254/150‑1 (MQ=255)
                                 aTTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTc                                                                                                                     >  2:4602631/1‑150 (MQ=255)
                                    gCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCt                                                                                                                  >  2:380048/1‑150 (MQ=255)
                                       aGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTAc                                                                                                               <  2:1964577/150‑1 (MQ=255)
                                        gCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCGCCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACa                                                                                                              <  1:4559671/150‑1 (MQ=255)
                                        gCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACa                                                                                                              <  1:5435587/150‑1 (MQ=255)
                                         caacaaCTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAg                                                                                                             >  1:2156964/1‑150 (MQ=255)
                                           acaacTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAt                                                                                                           >  2:3652709/1‑150 (MQ=255)
                                               cTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctg                                                                                                       <  1:3652709/150‑1 (MQ=255)
                                                tCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgc                                                                                                      >  1:872673/1‑150 (MQ=255)
                                                tCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgc                                                                                                      >  1:96167/1‑150 (MQ=255)
                                                  cTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgctg                                                                                                    >  2:2902111/1‑150 (MQ=255)
                                                   tGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgctgc                                                                                                   >  1:5205951/1‑150 (MQ=255)
                                                   tGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgctgc                                                                                                   >  1:3994488/1‑150 (MQ=255)
                                                    ggAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCa                                                                                                  <  1:976843/150‑1 (MQ=255)
                                                    ggAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCa                                                                                                  <  1:2368059/150‑1 (MQ=255)
                                                     gAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGTCTGAAAGGCTCCCTTACAGATGCTGCTGCAg                                                                                                 <  2:3691639/150‑1 (MQ=255)
                                                      aGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGt                                                                                                <  1:4493430/150‑1 (MQ=255)
                                                       gCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTg                                                                                               >  1:5009113/1‑150 (MQ=255)
                                                          cGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCt                                                                                            >  1:2327518/1‑150 (MQ=255)
                                                           gTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTc                                                                                           <  2:523671/150‑1 (MQ=255)
                                                             acgacgTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGa                                                                                         <  2:2614450/150‑1 (MQ=255)
                                                               gacgTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGc                                                                                                                                                 >  1:3356662/1‑92 (MQ=255)
                                                                 cgTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGaccac                                                                                     <  2:807827/150‑1 (MQ=255)
                                                                     tATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATc                                                                                 >  1:3841227/1‑150 (MQ=255)
                                                                     tATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATc                                                                                 >  2:5004822/1‑150 (MQ=255)
                                                                         gTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCa                                                                             <  1:1522533/150‑1 (MQ=255)
                                                                         gTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCa                                                                             >  2:4802968/1‑150 (MQ=255)
                                                                         gTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCa                                                                             <  2:2788906/150‑1 (MQ=255)
                                                                          tGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAg                                                                            >  1:349666/1‑150 (MQ=255)
                                                                            aCCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGc                                                                          >  1:2699838/1‑150 (MQ=255)
                                                                              ccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCt                                                                        >  1:1295334/1‑150 (MQ=255)
                                                                              ccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCt                                                                        >  2:3648534/1‑150 (MQ=255)
                                                                              ccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCt                                                                        >  1:4595514/1‑150 (MQ=255)
                                                                               cACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCtt                                                                       >  1:2560887/1‑150 (MQ=255)
                                                                               cACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCtt                                                                       <  2:4348628/150‑1 (MQ=255)
                                                                                aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTc                                                                      >  2:871576/1‑150 (MQ=255)
                                                                                  gCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCAt                                                                    >  2:1131443/1‑150 (MQ=255)
                                                                                   cAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATc                                                                   <  2:5384573/150‑1 (MQ=255)
                                                                                        tCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcct                                                              >  1:3087428/1‑150 (MQ=255)
                                                                                        tCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcct                                                              >  1:1368716/1‑150 (MQ=255)
                                                                                          aCCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctcc                                                            >  2:4137678/1‑150 (MQ=255)
                                                                                           ccATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctcct                                                           <  2:3260827/150‑1 (MQ=255)
                                                                                                  cGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCgatgat                                                    <  1:2505825/150‑1 (MQ=255)
                                                                                                  cGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAAACGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCgatgat                                                    >  1:777130/1‑150 (MQ=255)
                                                                                                   gCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATc                                                   <  1:1047852/150‑1 (MQ=255)
                                                                                                         tACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAccc                                             <  1:1131443/150‑1 (MQ=255)
                                                                                                         tACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAccc                                             <  1:2902111/150‑1 (MQ=255)
                                                                                                         tACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAccc                                             <  2:4473110/150‑1 (MQ=255)
                                                                                                           cAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCAt                                           <  2:3929712/150‑1 (MQ=255)
                                                                                                             gCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATAc                                         <  2:2363632/150‑1 (MQ=255)
                                                                                                             gCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATAc                                         <  2:933192/150‑1 (MQ=255)
                                                                                                                 cTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCAt                                     >  1:166592/1‑150 (MQ=255)
                                                                                                                 cTGATCGCCTGACGTTGAACGAAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCAt                                     <  2:3350281/150‑1 (MQ=255)
                                                                                                                  tGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATg                                    >  2:3618207/1‑150 (MQ=255)
                                                                                                                   gATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGa                                   <  2:5096299/150‑1 (MQ=255)
                                                                                                                    aTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGaa                                  <  1:2825700/150‑1 (MQ=255)
                                                                                                                      cGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGt                                >  2:2473264/1‑150 (MQ=255)
                                                                                                                        ccTGACGTTGAACGCAAAAATCCGCACGCAGGTGCGGTTTTTTCATGGATGCCTGAATGGCTCGCTTACAGATGGTGGTGCAGTGCCTCGAGCAGCATCGGCAGGCCTTCTTGGGCCTCGTCGATGATCGATGCCATAGGCTTGAAGTct                              >  2:5566481/1‑149 (MQ=255)
                                                                                                                         cTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCAt                             <  1:4268624/150‑1 (MQ=255)
                                                                                                                         cTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTACATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCAt                             <  2:4981981/150‑1 (MQ=255)
                                                                                                                          tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATg                            >  1:1720349/1‑150 (MQ=255)
                                                                                                                           gACGTTGAACGCAAAAAACCGCACCAAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATgg                           <  1:3612927/150‑1 (MQ=255)
                                                                                                                             cGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGt                         <  2:2116812/150‑1 (MQ=255)
                                                                                                                               ttGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCa                       >  2:4573228/1‑150 (MQ=255)
                                                                                                                                 gAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACt                     >  2:347720/1‑150 (MQ=255)
                                                                                                                                   catgaaaaaaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTcc                   >  2:3583523/5‑150 (MQ=255)
                                                                                                                                    atgaaaaaaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTccc                  >  1:22155/4‑150 (MQ=255)
                                                                                                                                     gCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCg                 <  2:1680197/150‑1 (MQ=255)
                                                                                                                                     gCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCg                 <  2:4238449/150‑1 (MQ=255)
                                                                                                                                      gaaaaaaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCgg                >  1:3251797/2‑150 (MQ=255)
                                                                                                                                      cAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACGGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCgg                <  2:5291021/150‑1 (MQ=255)
                                                                                                                                       aaaaaaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCggg               >  1:3947720/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCggg               >  1:2078966/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCggg               >  1:2582678/1‑150 (MQ=255)
                                                                                                                                         aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGAt             >  1:4808162/1‑150 (MQ=255)
                                                                                                                                         aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGAt             >  1:1022351/1‑150 (MQ=255)
                                                                                                                                            aCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCAc          <  1:1863972/150‑1 (MQ=255)
                                                                                                                                             ccGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCAcc         <  2:3762610/150‑1 (MQ=255)
                                                                                                                                             ccGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCAcc         >  2:344567/1‑150 (MQ=255)
                                                                                                                                              cGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCg        >  2:4314202/1‑150 (MQ=255)
                                                                                                                                              cGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCg        <  1:5485922/150‑1 (MQ=255)
                                                                                                                                                   ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGc   >  2:61788/1‑150 (MQ=255)
                                                                                                                                                   ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGc   <  2:1334717/150‑1 (MQ=255)
                                                                                                                                                     aGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGc       <  2:2055963/144‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACAACGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCT  >  NC_002947/4772664‑4772961

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