Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,760,108 G→C 7.3% 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,760,1080GC 7.3% 527.0 / 22.7 162intergenic (+1/+20)ompQ/pvdPouter membrane pyoverdine efflux protein/pyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base G (70/79);  new base C (9/3);  total (79/83)
Fisher's exact test for biased strand distribution p-value = 7.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

ACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGATCC  >  NC_002947/4759961‑4760257
                                                                                                                                                   |                                                                                                                                                     
aCCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGc                                                                                                                                                     <  1:4401126/150‑1 (MQ=255)
   ttGCTGTCCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCGAACTGCAGGCCAGCGTGGCCCTTTACAAAGCGTTGGGCGGTGGGTGGCAAGCGGATAGCCCGTCGCGCTAGGGCGtg                                                                                                                                                  >  2:4968854/1‑149 (MQ=255)
    tGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCa                                                                                                                                                 <  2:4460414/150‑1 (MQ=255)
    tGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCa                                                                                                                                                 <  1:957360/150‑1 (MQ=255)
      cTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCggg                                                                                                                                               <  1:2051792/150‑1 (MQ=255)
      cTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCggg                                                                                                                                               <  1:2038316/150‑1 (MQ=255)
      cTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAgg                                                                                                                                               <  1:603798/150‑1 (MQ=255)
      cTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCCCTAGCGCGTCggg                                                                                                                                               <  2:4164321/150‑1 (MQ=255)
       tGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGc                                                                                                                                              >  2:2369227/1‑150 (MQ=255)
           cGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATc                                                                                                                                          >  1:4993507/1‑150 (MQ=255)
           cGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATc                                                                                                                                          >  1:4183454/1‑150 (MQ=255)
            gTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCt                                                                                                                                         <  1:2666272/150‑1 (MQ=255)
                tCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGgtg                                                                                                                                     <  2:2088560/150‑1 (MQ=255)
                       ccAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGa                                                                                                                              >  2:1202484/1‑150 (MQ=255)
                              aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTATGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                                       <  1:2523996/150‑1 (MQ=255)
                              aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCACCTGCAGGCCAGCGTGGCCCTTTACAACGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                                       <  2:2637265/150‑1 (MQ=255)
                              aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                                       <  1:3940422/150‑1 (MQ=255)
                              aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                                       <  2:3041182/150‑1 (MQ=255)
                              aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                                       <  2:3041205/150‑1 (MQ=255)
                              aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                                       <  1:1360918/150‑1 (MQ=255)
                                ccTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAACGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCt                                                                                                                     <  2:3332524/150‑1 (MQ=255)
                                ccTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCt                                                                                                                     >  1:250069/1‑150 (MQ=255)
                                ccTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCt                                                                                                                     <  2:3332417/150‑1 (MQ=255)
                                 cTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTg                                                                                                                    >  1:4696277/1‑150 (MQ=255)
                                       cAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTc                                                                                                              >  1:2905266/1‑150 (MQ=255)
                                       cAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTc                                                                                                              <  1:544756/150‑1 (MQ=255)
                                       cAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTc                                                                                                              <  2:1491911/150‑1 (MQ=255)
                                         ggCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAg                                                                                                            <  2:2476721/150‑1 (MQ=255)
                                            ccAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTg                                                                                                         >  1:3131855/1‑150 (MQ=255)
                                             cAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTgg                                                                                                        >  1:749006/1‑150 (MQ=255)
                                              aGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGc                                                                                                       <  2:640857/150‑1 (MQ=255)
                                              aGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCCGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGc                                                                                                       <  1:984555/150‑1 (MQ=255)
                                                gACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                     <  1:1495817/150‑1 (MQ=255)
                                                gACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                     <  2:896317/150‑1 (MQ=255)
                                                gACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                                     <  1:4817308/150‑1 (MQ=255)
                                                   ctgcagGCGGGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCa                                                                                                  >  1:3042727/3‑150 (MQ=255)
                                                   ctgcagGCGGGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCa                                                                                                  >  1:770077/3‑150 (MQ=255)
                                                     gcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATa                                                                                                >  1:4032920/1‑150 (MQ=255)
                                                     gcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATa                                                                                                >  1:354672/1‑150 (MQ=255)
                                                        ggCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCa                                                                                             <  1:2608258/150‑1 (MQ=255)
                                                         gCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAg                                                                                            >  1:383465/1‑150 (MQ=255)
                                                         gCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAg                                                                                            <  2:4286532/150‑1 (MQ=255)
                                                           ccGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAggg                                                                                          <  2:3910122/150‑1 (MQ=255)
                                                             gCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGggc                                                                                        >  2:1273246/1‑150 (MQ=255)
                                                                          ccAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                           <  2:2234755/150‑1 (MQ=255)
                                                                           cAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTg                                                                          <  2:3170029/150‑1 (MQ=255)
                                                                            aaCTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGc                                                                         <  1:1835305/150‑1 (MQ=255)
                                                                               tGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTa                                                                      >  2:2726893/1‑150 (MQ=255)
                                                                                  aGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCa                                                                   >  1:1805453/1‑150 (MQ=255)
                                                                                    gCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATTGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCCGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGc                                                                 <  2:511084/150‑1 (MQ=255)
                                                                                    gCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGc                                                                 >  1:1932502/1‑150 (MQ=255)
                                                                                             gCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAg                                                        <  2:1592285/150‑1 (MQ=255)
                                                                                                cTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAgcgc                                                     >  2:2470513/1‑150 (MQ=255)
                                                                                                 tttACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGct                                                    <  2:1385824/150‑1 (MQ=255)
                                                                                                  ttACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGctc                                                   <  2:4696277/150‑1 (MQ=255)
                                                                                                    aCAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGctcct                                                 >  1:4038687/1‑150 (MQ=255)
                                                                                                      aaaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcg                                               >  2:3663477/1‑150 (MQ=255)
                                                                                                      aaaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcg                                               >  2:1259802/1‑150 (MQ=255)
                                                                                                      aaaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcg                                               >  1:5096623/1‑150 (MQ=255)
                                                                                                      aaaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcg                                               >  1:1792629/1‑150 (MQ=255)
                                                                                                        aGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgtc                                             <  2:1897657/150‑1 (MQ=255)
                                                                                                         gCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgtcg                                            >  2:3082923/1‑150 (MQ=255)
                                                                                                          ccTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGc                                           >  1:1783921/1‑150 (MQ=255)
                                                                                                           cTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCt                                          >  2:2015743/1‑150 (MQ=255)
                                                                                                           cTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCt                                          <  1:367936/150‑1 (MQ=255)
                                                                                                            ttGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTg                                         >  2:4222611/1‑150 (MQ=255)
                                                                                                            ttGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTg                                         >  2:4224304/1‑150 (MQ=255)
                                                                                                                    tggctggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAg                                 >  1:1603915/1‑150 (MQ=255)
                                                                                                                     ggctggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAgg                                <  2:1805453/150‑1 (MQ=255)
                                                                                                                      gctggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgc                               <  2:472824/150‑1 (MQ=255)
                                                                                                                        tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgcgc                             >  2:1994648/1‑150 (MQ=255)
                                                                                                                        tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgcgc                             >  1:1860175/1‑150 (MQ=255)
                                                                                                                             accagATAGCCTGACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCCCCTCGTCGCTGTAAATCAGGCGCAAGGc                        >  2:937566/5‑150 (MQ=39)
                                                                                                                             aGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATAGCAGGGGCGGCGCGGCACTTTTCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCATGCGAAAGGc                        <  2:5096623/150‑1 (MQ=255)
                                                                                                                              gCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGcc                       <  1:1338503/150‑1 (MQ=255)
                                                                                                                              gCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGcc                       >  1:2395591/1‑150 (MQ=255)
                                                                                                                               cagATAGCCTGACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCa                      >  1:3811522/3‑150 (MQ=255)
                                                                                                                               cGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCa                      <  1:2058503/150‑1 (MQ=255)
                                                                                                                                ggATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAg                     >  2:1265130/1‑150 (MQ=255)
                                                                                                                                agATAGCCTGACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAg                     >  2:314088/2‑150 (MQ=255)
                                                                                                                                 gATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGc                    >  1:3164879/1‑150 (MQ=255)
                                                                                                                                     gCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCa                <  2:2114175/150‑1 (MQ=255)
                                                                                                                                     gCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCa                <  2:872535/150‑1 (MQ=255)
                                                                                                                                      cctgACGCCCTAGCGCGTCGGGCTATCCGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCAt               >  2:4566025/4‑150 (MQ=255)
                                                                                                                                      cccGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCAt               >  2:3182949/1‑150 (MQ=255)
                                                                                                                                       ccGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATg              <  2:1113637/150‑1 (MQ=255)
                                                                                                                                       ccGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATg              <  1:2688373/150‑1 (MQ=255)
                                                                                                                                            gcgcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTc         <  2:4478194/150‑1 (MQ=255)
                                                                                                                                            gcccTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTc         >  2:390267/4‑150 (MQ=255)
                                                                                                                                             cgcgAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTcc        <  2:4183444/150‑1 (MQ=255)
                                                                                                                                             cgcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTcc        <  1:3582423/150‑1 (MQ=255)
                                                                                                                                             cgcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTAACAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTcc        <  1:4969212/150‑1 (MQ=255)
                                                                                                                                              gcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCg       >  1:2597161/1‑150 (MQ=255)
                                                                                                                                              gcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCg       >  1:2597158/1‑150 (MQ=255)
                                                                                                                                              gcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCg       <  2:4993507/150‑1 (MQ=255)
                                                                                                                                                tAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGa     <  1:2937344/150‑1 (MQ=255)
                                                                                                                                                tAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGa     >  2:1780286/1‑150 (MQ=255)
                                                                                                                                                tAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGa     >  1:3583597/1‑150 (MQ=255)
                                                                                                                                                  gggCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGATc   <  2:277837/150‑1 (MQ=255)
                                                                                                                                                  gggCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGATc   <  1:1202484/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCCGGATCC  >  NC_002947/4759961‑4760257

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: